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This third edition, like the first two, is an introduction to the basic methods, theory, 
and applications of differential equations. A knowledge of elementary calculus is 
presupposed. 

The detailed style of presentation that characterized the previous editions of the text 
has been retained. Many sections have been taken verbatim from the second edition, 
while others have been rewritten or rearranged with the sole intention of making them 
clearer and smoother As in the earlier editions, the text contains many thoroughly 
worked out examples. Also, a number of new exercises have been added, and assorted 
exercise sets rearranged to make them more useful in teaching and learning. 

The book is divided into two main parts. The first part (Chapters 1 through 9) deals 
with the material usually found in a one-semester introductory course in ordinary 
differential equations. This part is also available separately as Introduction to Ordinary 
Differential Equations , Third Edition (John Wiley & Sons, New' York, 1980). The 
second part of the present text (Chapters 10 through 14) introduces the reader to 
certain specialized and more’ advanced methods and provides an introduction to 
fundamental theory. The table of contents indicates just what topics are treated. 

The following additions and modifications are specifically noted. 

1. Material emphasizing the second-order linear equation has been inserted at 
appropriate places in Section 4.1. 

2. New illustrative examples, including an especially detailed introductory one, have 
been written to clarify the Method of Undetermined Coefficients in Section 4.3, 
and a useful table has also been supplied. 

3. Matrix multiplication and inversion have been added to the introductory 
material on linear algebra in Section 7.5. 

4. Additional applications now appear in the text in Sections 3.3 and 1.2. 

5. Section 7.6 is a completely new section on the application of matrix algebra to the 
solution of linear systems with constant coefficients in the special case of two 
equations in two unknown functions. The theory that occupied this section in the 
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second edition now appears in Chapter 1 1 (see note 9 following). We believe that 
this change represents a major improvement for both introductory and inter- 
mediate courses. 

6. Section 7.7 extends the matrix method of Section 7.6 to the case of linear systems 
with constant coefficients involving n equations in n unknown functions. Several 
detailed examples illustrate the method for the case n = 3. 

7. Both revised and new material on the Laplace Transform of step functions, 
translated functions, and periodic functions now appears in Section 9.1. 

8. The basic existence theory for systems and higher-order equations, formerly 
located at the beginning of Chapter 11, has now been placed at the end of 
Chapter 10. This minor change has resulted in better overall organization. 

9. Chapter 11, the Theory of Linear Differential Equations, has been changed 
considerably. Sections 1 1.1 through 1 1.4 present the fundamental theory of linear 
systems. Much of this material was found in Section 7,6 in the second edition, and 
some additional results are also included here. Sections 1L5 through 11.7 now 
present the basic theory of the single nth-order equation, making considerable use 
of the material of the preceding sections. Section 1 1 .8 introduces second-order 
self-adjoint equations and proceeds through the fundamentals of classical Sturm 
Theory. We believe that the linear theory is now presented more coherently than 
in the previous edition. 

10. An appendix presents, without proof, the fundamentals of second and third order 
determinants. 

The book can be used as a text in several different types of courses. The more or less 
traditional one-semester introductory course could be based on Chapter 1 through 
Section 7.4 of Chapter 7 if elementary applications are to be included. An alternative 
one-semester version omitting applications but including numerical methods and 
Laplace transforms could be based on Chapters 1, 2, 4, 6, 7, 8, and 9. An introductory 
course designed to lead quickly to the methods of partial differential equations could 
be based on Chapters 1, 2 (in part), 4, 6, 12, and 14. 

The book can also be used as a text in various intermediate courses for juniors and 
seniors who have already had a one-semesj^r introduction to the subject. An interme- 
diate course emphasizing further methods could be based on Chapters 8, 9, 12, 13, 
and 14. An intermediate course designed as an introduction to fundamental theory 
could be based on Chapters 10 through 14. We also note that Chapters 13 and 14 can 
be interchanged advantageously. 

I am grateful to several anonymous reviewers who made useful comments and 
suggestions. I thank my colleagues William Bonnice and Robert O. Kimball for helpful 
advice. I also thank my son, Shepley L. Ross, II, graduate student in mathematics, 
University of Rochester, Rochester, New York, for his careful reviewing and helpful 
suggestions. 

I am grateful to Solange Abbott for her excellent typing. I am pleased to record my 
appreciation to Editor Gary Ostedt and the Wiley staff for their constant helpfulness 
and cooperation. 

As on several previous occasions, the most thanks goes to my wife who offered 
encouragement, understanding, patience, and help in many different ways. Thanks, 
Gin. 
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CHAPTER ONE 

Differential Equations and Their Solutions 



The subject of differential equations constitutes a large and very important branch of 
modern mathematics. From the early days of the calculus the subject has been an area 
of great theoretical research and practical applications, and it continues to be so in our 
day. This much stated, several questions naturally arise. Just what is a differential 
equation and what does it signify? Where and how do differential equations originate 
and of what use are they? Confronted with a differential equation, what does one do 
with it, how does one do it, and what are the results of such activity? These questions 
indicate three major aspects of the subject: theory, method, and application. The 
purpose of this chapter is to introduce the reader to the basic aspects of the subject and 
at the same time give a brief survey of the three aspects just mentioned. In the course of 
the chapter, we shall find answers to the general questions raised above, answers that 
will become more and more meaningful as we proceed with the study of differential 
equations in the following chapters. 



1-1 CLASSIFICATION OF DIFFERENTIAL EQUATIONS; THEIR ORIGIN AND 
APPLICATION 

A. Differential Equations and Their Classification 
DEFINITION 

An equation involving derivatives of one or more dependent variables with respect to one or 
more independent variables is called a differential equation.* 



* In connection with this basic definition, we do not include in the class of differential equations those 
equations that are actually derivative identities. For example, we exclude such expressions as 



d 

— (e~) - ae*\ 
dx 



— (uv) - u 
dx 



dv du 

- — - + v — v 

dx dx 



and so forth. 



3 





4 DIFFERENTIAL EQUATIONS AND THEIR SOLUTIONS 



► Example 1.1 



For examples of differential equations we list the following: 

d 2 y 



dx 2 
d 4 x 

di 4 



( dy \ 2 

+ xyf -p- 1 - 0, 



_ d 2 x 

+ 5 2 + 3x ; 

di 2 



sin t , 



(M) 

( 1 . 2 ) 



d 2 u 

dx 2 



dv dv 

— + — * y, 

ds dt 

d 2 u d 2 u 
+ dy 2 + dz 2 



0 . 



(1.3) 

0.4) 



From the brieflist of differential equations in Example 1.1 it is clear that the various 
variables and derivatives involved in a differential equation can occur in a variety of 
ways. Clearly some kind of classification must be made. To begin with, we classify 
differential equations according to whether there is one or more than one independent 
variable involved. 



DEFINITION 

A differential equation involving ordinary derivatives of one or more dependent variables 
with respect to a single independent variable is called an ordinary differential equation. 



► Example 1.2 

Equations (1.1) and (1.2) arfc ordinary differential equations. In Equation (1.1) the 
variable x is the single independent variable, and y is a dependent variable. In Equation 
(1,2) the independent variable is t, whereas x is dependent. 

DEFINITION 

A differential equation involving partial derivatives of one or more dependent variables 
with respect to more than one independent variable is called a partial differential 
equation. 



► Example 1 .3 

Equations (1.3) and (1.4) are partial differential equations. In Equation (1.3) the 
variables s and t are independent variables and u is a dependent variable. In Equation 
(1.4) there are three independent variables: x, y, and z ; in this equation u is dependent. 

We further classify differential equations, both ordinary and partial, according to the 
order of the highest derivative appearing in the equation. For this purpose we give the 
following definition. 
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DEFINITION 

The order of the highest ordered derivative involved in a differential equation is called the 
order of the differential equation. 



► Example 1 .4 

The ordinary differentia] equation (1,1) is of the second order, since the highest 
derivative involved is a second derivative. Equation (1.2) is an ordinary differential 
equation of the fourth order. The partial differential equations (1.3) and (1.4) are of the 
first and second orders, respectively. 

Proceeding with our study of ordinary differential equations, we now introduce the 
important concept of linearity applied to such equations. This concept will enable us to 
classify these equations still further. 



DEFINITION 



A linear ordinary differential equation of order n, in the dependent variable y and the 
independent variable x , is an equation that is in , or can be expressed in, the form 



, d n y d n i y 

«.W S ; + «,W S ^T + 



dy 

+ a^^x)— + a n (x)y = Hx), 



where a 0 is not identically zero . 



Observe (1) that the dependent variable yand its various derivatives occur to the first 
degree only, (2) that no products of y and/or any of its derivatives are present, and (3) 
that no transcendental functions of y and/or its derivatives occur. 



► Example 1 .5 



The following ordinary differential equations are both linear. In each case y is the 
dependent variable. Observe that y and its various derivatives occur to the first degree 
only and that no products of y and/or any of its derivatives are present. 



&♦*£♦*-* 



(1.5) 



d*y , d 3 y , dy 

i? +x d? +x £= x ' 



( 1 . 6 ) 



DEFINITION 

A nonlinear ordinary differential equation is an ordinary differential equation that is not 
linear . 
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► Example 1 .6 

The following ordinary differential equations are all nonlinear: 

^ + 5 (S) +6y = 0, (1 - 8) 

d 2 y dy 

d? + 5 y i +6y =°- ™ 

Equation (1.7) is nonlinear because the dependent variable y appears to the second 
degree in the term by 1 . Equation (1.8) owes its nonlinearity to the presence of the term 
5 (dy/dx)*> which involves the third power of the first derivative, Finally, Equation (1.9) 
is nonlinear because of the term 5 y{dyfdx) y which involves the product of the 
dependent variable and its first derivative. 

Linear ordinary differential equations are further classified according to the nature 
of the coefficients of the dependent variables and their derivatives. For example, 
Equation (1.5) is said to be linear with constant coefficients » while Equation ( 1 .6) is linear 
with variable coefficients . 



B. Origin and Application of Differential Equations 

Having classified differential equations in various ways, let us now consider briefly 
where, and how, such equations actually originate. In this way we shall obtain some 
indication of the great variety of subjects to which the theory and methods of 
differential equations may be applied. 

Differential equations occur in connection with numerous problems that are 
encountered in the various branches of science and engineering. We indicate a few such 
problems in the following list, which could easily be extended to fill many pages. 

1. "The problem of determining the motion of a projectile, rocket, satellite, or planet. 

2. The problem of determining the charge or current in an electric circuit. 

3. The problem of the conduction of heat in a rod or in a slab. 

4. The problem of determining the vibrations of a wire or a membrane. 

5. The study of the rate of decomposition of a radioactive substance or the rate of 
growth of a population. 

6. The study of the reactions of chemicals, 

7. The problem of the determination of curves that have certain geometrical 
properties. 

The mathematical formulation of such problems give rise to differential equations. 
But just how does this occur? In the situations under consideration in each of the above 
problems the objects involved obey certain scientific laws. These laws involve various 
rates of change of one or more quantities with respect to other quantities. Let us re- 
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call that such rates of change are expressed mathematically by derivatives. In the 
mathematical formulation of each of the above situations, the various rates of change 
are thus expressed by various derivatives and the scientific laws themselves become 
mathematical equations involving derivatives, that is, differential equations. 

In this process of mathematical formulation, certain simplifying assumptions 
generally have to be made in order that the resulting differential equations be tractable. 
For example, if the actual situation in a certain aspect of the problem is of a relatively 
complicated nature, we are often forced to modify this by assuming instead an 
approximate situation that is of a comparatively simple nature. Indeed, certain 
relatively unimportant aspects of the problem must often be entirely eliminated. The 
result of such changes from the actual nature of things means that the resulting 
differential equation is actually that of an idealized situation. Nonetheless, the 
information obtained from such an equation is of the greatest value to the scientist. 

A natural question now is the following: How does one obtain useful information 
from a differential equation? The answer is essentially that if it is possible to do so, one 
solves the differential equation to obtain a solution; if this is not possible, one uses the 
theory of differential equations to obtain information about the solution. To 
understand the meaning of this answer, we must discuss what is meant by a solution of 
a differential equation; this is done in the next section. 



Exercises 

Classify each of the following differential equations as ordinary or partial differential 
equations; state the order of each equation; and determine whether the equation under 
consideration is linear or nonlinear. 



dy y 

— + x z y ~ xe x . 
dx 


2 . 


d*y d 2 y c dy 

S s + 4 C'~ 5 * +3, ’* S ‘"*' 


d 2 u d 2 u 
dx 2 + dy 2 


4. 


x 2 dy + y 2 dx = 0. 


g + ,(£)\*-a 


6. 


d 4 u d 2 u d 2 u 

dx 1 dy 1 + dx 2 + dy 2 + U 


d 2 y 


8. 


d 2 y 

dx 2 + x sin ,v = 0. 


d 6 x (d*x\{d 3 x\ 


10. 


fdrV jd 2 r 

\ds J ds 2 + 



1.2 SOLUTIONS 

A. Nature of Solutions 

We now consider the concept of a solution of the nth-order ordinary differential 
equation. 
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DIFFERENTIAL EQUATIONS AND THEIR SOLUTIONS 



DEFINITION 



Consider the nth-order ordinary differential equation 



F 




dy_ 

'dx 



d*y\ 



= 0 , 



( 1 . 10 ) 



dy 

where F is a real function of its (n -f 2) arguments x f y,— 



d n y 
’ d x" 



L Let f be a real function defined for all x in a real interval I and having an nth 
derivative ( and hence also all lower ordered derivatives ) for all x e L The function f is 
called an explicit solution of the differential equation (1.10) on I if it fulfills the following 
two requirements : 

FO./W, / ,n) M] (A) 

is defined for all x c I, and 

F[x> f(x), f'(x ), . . . , f {n) (xf] = 0 (B) 

for all xe L That is, the substitution of f(x) and its various derivations for y and its 
corresponding derivatives, respectively , in (1.10) reduces (1.10) to an identity on / . 

2 . A relation g(x, y) - 0 is called an implicit solution of(L)O) if this relation defines 
at least one real function f of the variable x on an interval l such that this function is an 
explicit solution of (1.10) on this interval . 

3. Both explicit solutions and implicit solutions will usually be called simply solutions. 

Roughly speaking, then, we may say that a solution of the differential equation (1.10) 
is a relation — explicit or implicit — between x and y , not containing derivatives, which 
identically satisfies (1.10). 



► Example 1.7 




The function / defined for all real x by 




f(x) = 2 sin x + 3 cos x 


o.ii) 


is an explicit solution of the differential equation 




d 2 y 


(1.12) 



for all real x . First note that / is defined and has a second derivative for all real x. Next 
observe that 

f\x) - 2 cos x — 3 sin x y 
f"(x) = -2 sin x - 3 cos x. 

- Upon substituting f"(x) for d 2 y/dx 2 and /(x)for yin the differential equation (1.1 2), it 
reduces to the identity 

(—2 sin x — 3 cos x) + (2 sin x + 3 cos x) — 0, 
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Thus the relation (1.15) formally satisfies the differential equation (1.1 4). Can we 
conclude from this alone that (1.15) is an implicit solution of (1,14)? The answer to this 
question is “no,” for we have no assurance from this that the relation (1.15) defines any 
function that is an explicit solution of (1.14) on any real interval L All that we have 
shown is that ( 1 . 1 5) is a relation between x and y that, upon implicit differentiation and 
substitution, formally reduces the differential equation (1.14) to a formal identity. It is 
catled a formal solution; it has the appearance of a solution; but that is all that we know 
about it at this stage of our investigation. 

Let us investigate a little further. Solving (1.15) for y, we find that 

±/-25-“^ 2 . 

Since this expression yields nonreal values of y for all real values of x, we conclude 
that the relation (1.15) does not define any real function on any interval. Thus the 
relation (LI 5) is not truly an implicit solution but merely a forma! solution of the 
differential equation (1.14). 

In applying the methods of the following chapters we shall often obtain relations that 
we can readily verify are at least formal solutions. Our main objective will be to gain 
familiarity with the methods themselves and we shall often be content to refer to the 
relations so obtained as “$olutions,”although we have no assurance that these relations 
are actually true implicit solutions. If a critical examination of the situation is required, 
one must undertake to determine whether or not these formal solutions so obtained are 
actually true implicit solutions which define explicit solutions. 

In order to gain further insight into the significance of differential equations and 
their solutions, we now examine the simple equation of the following example. 

► Example 1 .9 

Consider the first-order differential equation 

£ = 2, (MO) 

The function f 0 defined for all real x by f 0 (x) — x 2 is a solution of this equation. So also 
are the functions /j,/ 2 , and / 3 defined for all real x by /,(x) - x 2 + 1 ,/ 2 (x) = x 2 + 2, 
and / 3 (x) = x 2 + 3, respectively. In fact, for each real number c, the function f c defined 
for all real x by 

f{x)^x 2 + c (1.17) 

is a solution of the differential equation (LI 6). In other words, the formula (1.17) defines 
an infinite family of functions, one for each real constant c, and every function of this 
family is a solution of (1.16). We call the constant c in (1.17) an arbitrary constant or 
parameter and refer to the family of functions defined by (LI 7) as a one-parameter 
faimily of solutions of the differential equation (1.16). We write this one-parameter 
family of solutions as 

y — x 2 + c. (1 .18) 

Although it is clear that every function of the family defined by (LIB) is a solution of 
(1.16), we have not shown that the family of functions defined by (1.18) includes all of 
the solutions of (1,16). However, we point out (without proof) that this is indeed the 
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case here; that is, every solution of (1.16) is actually of the form (1.18) for some 
appropriate real number c. 

Note. We must not conclude from the last sentence of Example 1 .9 that every first- 
order ordinary differential equation has a so-called one-parameter family of solutions 
which contains all solutions of the differential equation, for this is by no means the case. 
Indeed, some first-order differential equations have no solution at all (see Exercise 7(a) 
at the end of this section), while others have a one-parameter family of solutions plus 
one or more “extra” solutions which appear to be “different” from all those of the 
family (see Exercise 7(b) at the end of this section). 

The differential equation of Example 1.9 enables us to obtain a better understanding 
of the analytic significance of differential equations. Briefly stated, the differential 
equation of that example defines functions , namely, its solutions. We shall see that this 
is the case with many other differential equations of both first and higher order. Thus 
we may say that a differential equation is merely an expression involving derivatives 
which may serve as a means of defining a certain set of functions: its solutions. Indeed, 
many of the now familiar functions originally appeared in the form of differential 
equations that define them. 

We now consider the geometric significance of differential equations and their 
solutions. We first recall that a real function F may be represented geometrically by a 
curve y = F(x) in the xy plane and that the value of the derivative of F at x, E'(x), may 
be interpreted as the slope of the curve y = F(x) at x. Thus the general first-order 
differential equation 

0 - 19 ) 

where / is a real function, may be interpreted geometrically as defining a slope /(x, y) at 
every point (x, y) at which the function / is defined. Now assume that the differential 
equation (1.19) has a so-called one-parameter family of solutions that, can be written in 
the form 

( 1 . 20 ) 

where c is the arbitrary constant or parameter of the family. The one-parameter family 
of functions defined by (1.20) is represented geometrically by a so-called one-parameter 
family of curves in the xy plane, the slopes of which are given by the differential 
equation (1.19). These curves, the graphs of the solutions of the differential equation 
(1.19), are called the integral curves of the differential equation (1.19). 



► Example 1.10 



Consider again the first-order differential equation 




0-16) 



of Example 1 .9. This differential equation may be interpreted as defining the slope 2x at 
the point with coordinates (x, y) for every real x. Now, we observed in Example 1 .9 that 
the differential equation (1.16) has a one-parameter family of solutions of the form 

y = x 2 4- c, (1.18) 
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where c is the arbitrary constant or parameter of the family. The one-parameter family 
of functions defined by (1.1 8) is represented geometrically by a one-parameter family of 
curves in the xy plane, namely, the family of parabolas with Equation (1.18). The slope 
of each of these parabolas is given by the differential equation (1. 1 6) of the family. Thus 
we see that the family of parabolas (1.18) defined by differential equation (1.16) is that 
family of parabolas, each of which has slope 2x at the point (. x , y) for every real x, and 
all of which have the y axis as axis. These parabolas are the integral curves of the 
differential equation (1.16). See Figure 1.1. 



B. Methods of Solution 

When we say that we shall solve a differential equation we mean that we shall find one 
or more of its solutions. How is this done and what does it really mean? The greater 
part of this text is concerned with various methods of solving differential equations. 
The method to be employed depends upon the type of differential equation under 
consideration, and we shall not enter into the details of specific methods here. 

But suppose we solve a differential equation, using one or another of the various 
methods. Does this necessarily mean that we have found an explicit solution / 
expressed in the so-called closed form of a finite sum of known elementary functions? 
That is, roughly speaking, when we have solved a differential equation, does this 
necessarily mean that we have found a “formula” for the solution? The answer is “no.” 
Comparatively few differential equations have solutions so expressible; in fact, a 
closed-form solution is really a luxury in differential equations. In Chapters 2 and 4 we 
shall consider certain types of differential equations that do have such closed-form 
solutions and study the exact methods available for finding these desirable solutions. 
But, as we have just noted, such equations are actually in the minority and we must 
consider what it means to “solve” equations for which exact methods are unavailable. 
Such equations are solved approximately by various methods, some of which are 
considered in Chapters 6 and 8. Among such methods are series methods, numerical 



l 
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methods, and graphical methods. What do such approximate methods actually yield? 
The answer to this depends upon the method under consideration. 

Series methods yield solutions in the form of infinite series; numerical methods give 
approximate values of the solution functions corresponding to selected values of the 
independent variables; and graphical methods produce approximately the graphs of 
solutions (the integral curves). These methods are not so desirable as exact methods 
because of the amount of work involved in them and because the results obtained from 
them are only approximate; but if exact methods are not applicable, one has no choice 
but to turn to approximate methods. Modern science and engineering problems 
continue to give rise to differential equations to which exact methods do not apply, and 
approximate methods are becoming increasingly more important. 



Exercises 

1. Show 1 that each of the functions defined in Column I is a solution of the 
corresponding differential equation in Column II on every interval a < x < b of 
the x axis. 

1 11 

(a) f(x) = x + 3e~* lbc + y = X + l 

(b) /(*) = 2e ix -5e**- 0-7^+12> = O 

(c) /(x) = e* + 2x 2 + 6x + 7 ~ - 3^ + 2y = 4x 2 

w) m-rh ,l+x!) S + 4 x l +2),=0 

2. (a) Show that x 3 -f 3 xy 2 = 1 is an implicit solution of the differential equation 

2xy{dy/dx) + x 2 + y 2 =*= 0 on the interval 0 < x < 1. 

(b) Show that 5x 2 y 2 — 2x 3 y 2 = 1 is an implicit solution of the differential 
equation x{dy/dx) + y = x 3 y 3 on the interval 0 < x < f . 

3. (a) Show that every function / defined by 

fix) = (x 3 -1- c)e 3x > 

where c is an arbitrary constant, is a solution of the differential equation 

ly = lx 2 e' 3x . 
ax 

(b) Show that every function / defined by 

f(x) = 2 + ce-^\ 

where c is an arbitrary constant, is a solution of the differential equation 

dy 
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4. (a) Show that every function / defined by f(x) =* c x e* x + c 2 e~ 2x 9 wherecj andc 2 
are arbitrary constants, is a solution of the differential equation 



d 2 y 2 dy 
dx 2 dx 



- 8y *= Q. 



(b) Show that every function # defined by g(x) = c x e 2 * + c 2 xe lx + c 3 e~ 2x , where 
c lf c 2 , and c 3 are arbitrary constants, is a solution of the differential equation 



<Py 0 d^y 

dx 3 dx 2 



- 4 ^ + 8)> = 0. 
dx 



5. (a) For certain values of the constant m the function / defined by f(x) = e mx is a 
solution of the differential equation 



d 3 y 
d ? 





12y = 0. 



Determine all such values of m. 

(b) For certain values of the constant n the function g defined by g(x) = x" is a 
solution of the differential equation 



d 3 y 



,fy 



dy 



x 3 -4 + 2x 2 -4 - lOx - %y = 0. 



dx 



dx 



dx 



Determine all such values of n. 

6. (a) Show that the function / defined by /(x) - (2x 2 + 2e 3x + 3)e~ 2x satisfies the 
differential equation 

dy _ - 

~ + 2y = 6e* 4* 4xe 2x 
\ dx 

and also the condition /( 0) = 5. 

(b) Show that the function / defined by /(x) = 3e 2x - 2xe 2x - cos 2x satisfies 
the differential equation 

d 2 y dy 

-pT - 4 -p- + 4 y - - 8 sin 2x 
dx" 1 dx 



and also the conditions that /(0) = 2 and /'(0) = 4. 



7. (a) Show that the first-order differential equation 



dy 

dx 



+ 1 * 1 + 1 =0 



has no (real) solutions. 

(b) Show that the first-order differential equation 



fdyV 

/ ~ 4 >^ 0 
\dx J 



has a one-parameter family of solutions of the form f{x) = (x + c) 2 , where c is 
an arbitrary constant, plus the “extra” solution #(x) = 0 that is not a member 
of this family f(x) = (x + c) 2 for any choice of the constant c. 
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1.3 INITIAL-VALUE PROBLEMS, BOUNDARY-VALUE PROBLEMS, AND 
EXISTENCE OF SOLUTIONS 

A. Initial-Value Problems and Boundary-Value Problems 

We shall begin this section by considering the rather simple problem of the following 
example. 



► Example 1.11 



Problem. Find a solution / of the differential equation 




such that at x = 1 this solution / has the value 4. 



0 21 ) 



Explanation. First let us be certain that we thoroughly understand this problem. 
We seek a real function / which fulfills the two following requirements: 



1. The function / must satisfy the differential equation (1.21). That is, the function / 
must be such that f'(x ) = 2x for ail real x in a real interval I. 

2. The function / must have the value 4 at x = 1. That is, the function/ must be such 
that /(l) = 4. 



Notation. The stated problem may be expressed in the following somewhat 
abbreviated notation: 




y{ l) = 4. 



In this notation we may regard y as representing the desired solution. Then the 
differential equation itself obviously represents requirement 1, and the statement 
y(l) = 4 stands for requirement 2. More specifically, the notation y(l) — 4 states 
that the desired solution y must have the value 4 at x — I; that is, y = 4 at x = 1. 



Solution. We observed in Example 1.9 that the differential equation (1.21) has a 
one-parameter family of solutions which we write as 

y = (1-22) 

where c is an arbitrary constant, and that each of these solutions satisfies requirement 1 . 
Let us now attempt to determine the constant c so that (1.22) satisfies requirement 2, 
that is, y = 4 at x = 1. Substituting x ^ 1, y = 4 into (1.22), we obtain 4 = 1 + c, and 
hence c — 3. Now substituting the value c = 3 thus determined back into (1.22), we 
obtain ^ 



y = x 2 + 3, 
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which is indeed a solution of the differential equation (1.21), which has the value 4 at 
x = 1. In other words, the function / defined by 

fix) -x 2 + 3, 

satisfies both of the requirements set forth in the problem. 

Comment on Requirement 2 and Its Notation . In a problem of this type, require- 
ment 2 is regarded as a supplementary condition that the solution of the differential 
equation must also satisfy. The abbreviated notation y(l) = 4, which we used to 
express this condition, is in some way undesirable, but it has the advantages of being 
both customary and convenient. 

In the application of both first- and higher-order differential equations the problems 
most frequently encountered are similar to the above introductory problem in that they 
involve both a differential equation and one or more supplementary conditions which 
the solution of the given differential equation must satisfy. If all of the associated 
supplementary conditions relate to one x value, the problem is called an initial-value 
problem (or one-point boundary-value problem). If the conditions relate to two dif- 
ferent x values, the problem is called a two-point boundary-value problem (or simply a 
boundary-value problem). We shall illustrate these concepts with examples and then 
consider one such type of problem in detail. Concerning notation, we generally employ 
abbreviated notations for the supplementary conditions that are similar to the ab- 
breviated notation introduced in Example 1.1 1. 

► Example 1.12 



d 2 y 

dx 



izi + y = °> 



rtU-3, 

/(l) = -4. 

This problem consists in finding a solution of the differential equation 

d 2 y 

d? +y = 0 ' 



which assumes the value 3 at x = 1 and whose first derivative assumes the value -4 at 
x = 1. Both of these conditions relate to one x value, namely, x = 1. Thus this is an 
initial-value problem. We shall see later that this problem has a unique solution. 



► Example 1 .1 3 




y(0) 



o, 
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In this problem we again seek a solution of the same differential equation, but this time 
the solution must assume the value 1 at x = 0 and the value 5 at x = tz/2. That is, the 
conditions relate to the two different x values, 0 and rc/2. This is a (two-point) 
boundary-value problem. This problem also has a unique solution; but the boundary- 
value problem 



J 2 y 

dx 1 



+ 



y-o, 



y( 0) - I. . y(n) = 5, 

has no solution at all! This simple fact may lead one to the correct conclusion that 
boundary-value problems are not to be taken lightly! 



We now turn to a more detailed consideration of the initial- value problem fora first- 
order differential equation. 



DEFINITION 



Consider the first-order differential equation 



^ = /(x>y)> (1-23) 

where f is a continuous function of x and y in some domain * D of the xy plane; and let 
(x 0 ,y 0 ) be a point of D . The initial-value problem associated with { 1.23 ) is to find 
a solution <f> of the differential equation {1 .23), defined on some real interval containing x 0 , 
and satisfying the initial condition 



4>(x 0 ) = Jo- 

In the customary abbreviated notation , this initial-value problem may be written 



dy 



y{* o) = y 0 - 

To solve this problem, we must find a function <f) that not only satisfies the differential 
equation (1.23) but that also satisfies the initial condition that it has the value y 0 when x 
has value x 0 . The geometric interpretation of the initial condition, and hence of the 
entire initial-valine problem, is easily understood. The graph of the desired solution 4> 
must pass through the point with coordinates (x 0 , y 0 ). That is, interpreted geometri- 
cally, the initial-value problem is to find an integral curve of the differential equation 
(1.23) that passes through the point (x 0 , y 0 ). 

The method of actually finding the desired solution 0 depends upon the nature of the 
differential equation of the problem, that is, upon the form of /(x, y). Certain special 
types of differential equations have a one-parameter family of solutions whose 
equation may be found exactly by following definite procedures (see Chapter 2). If the 
differential equation of the problem is of some such special type, one first obtains the 
equation of its one-parameter family of solutions and then applies the initial condition 



* A domain is an open, connected set. For those unfamiliar with such concepts, D may be regarded as the 
interior of some simple closed curve in the plane. 
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to this equation in an attempt to obtain a “particular” solution <f> that satisfies the entire 
initial-value problem. We shall explain this situation more precisely in the next 
paragraph. Before doing so, however, we point out that in general one cannot find 
“the equation of a one-parameter family of solutions of the differential equation; 
approximate methods must then be used (see Chapter 8). 

Now suppose one can determine the equation 

9(x,y,c) = 0 (1.24) 

of a one-parameter family of solutions of the differential equation of the problem. 
Then, since the initial condition requires that y ~ y 0 at x - x 0 , we let x = x 0 and 
y - y 0 in (1.24) and thereby obtain 

3(xoiy 0 >c) ~ 0. 

Solving this for c, in general we obtain a particular value of c which we denote here by 
c 0 ■ We now replace the arbitrary constant c by the particular constant c 0 in ( 1 .24), thus 
obtaining the particular solution 

g(x t y,c 0 ) = 0. 

The particular explicit solution satisfying the two conditions (differential equation and 
initial condition) of the problem is then determined from this, if possible. 

We have already solved one initial-value problem in Example 1.11. We now give 
another example in order to illustrate the concepts and procedures more thoroughly. 



► Example 1.14 



Solve the initial-value problem 



dy x 
dx = 



(1.25) 



y( 3 ) = 4 , 



(1.26) 



given that the differential equation (1.25) has a one-parameter family of solutions which 
may be written in the form 

x 2 + | 2 = c 2 (1.27) 

The condition ( 1 .26) means that we seek the solution of ( 1 .25) such that y - 4 at x - 3. 
Thus the pair of values (3,4) must satisfy the relation (1.27). Substituting x - 3 and 
y = 4 into (1.27), we find 

9+ 16 = c 2 or c 1 = 25. 

Now substituting this value of c 2 into (1.27), we have 

x 2 + y 2 ~ 25. 

Solving this for y, we obtain 

±y/25 - x 2 . 

Obviously the positive sign must be chosen to give y the value + 4 at x = 3. Thus the 
function / defined by 



f{x) = 725 - x 5 , 



— 5 < x < 5, 
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is the solujiofHrf Jhe problem. In the usual abbreviated notation, we write this solution 
as y = y/25 — x\ 



B. Existence of Solutions 



In Example 1.14 we were able to find a solution of the initial-value problem under 
consideration. But do all initial-value and boundary-value problems have solutions? 
We have already answered this question in the negative, for we have pointed out that 
the boundary-value problem 



d 2 y 

dx 



2 + y - o, 



y( o) - 1, 

y(7t) = 5, 

mentioned at the end of Example 1.13, has no solution! Thus arises the question of 
existence of solutions: given an initial-value or boundary : value problem, does it 
actually have a solution? Let us consider the question for the initial-value problem 
defined on page 1 7. Here we can give a definite answer. Every initial-value problem that 
satisfies the definition on page 17 has at least one solution. 

But now another question is suggested, the question of uniqueness. Does such a 
problem ever have more than one solution? Let us consider the initial-value problem 



y( 0) - o. 



One may verify that the functions f 1 and f 2 defined, respectively, by 



/i (x) = 0 for all real x; 

and 

/ 2 (x) = (fx) 3 ' 2 , x > 0; f 2 (x) - 0, x < 0; 

are both solutions of this initial-value problem! In fact, this problem has infinitely many 
solutions! The answer to the uniqueness question is clear: the initial-value problem, as 
stated, need not have a unique solution. In order tb ensure uniqueness, some additional 
requirement must certainly be imposed. We shall see what this is in Theorem 1.1, which 
we shall now state. 



THEOREM 1.1. BASIC EXISTENCE AND UNIQUENESS THEOREM 



Hypothesis. Consider the differential equation 



dx 



- /(*. y). 



(1.28) 
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where 

/. The function f is a continuous function of x and y in some domain D of the xy plane , 
and 

2 . The partial derivative df/dy is also a continuous function of x and y in D; and let 
(x 0 , yo) he a point in D, 

Conclusion. There exists a unique solution 0 of the differential equation (/.2S), 
defined on some interval |x — x 0 | < h, where h is sufficiently small , that satisfies the 
condition 



<t>(x 0 )=yo • (1.29) 

Explanatory Remarks . This basic theorem is the first theorem from the theory of 
differential equations which we have encountered. We shall therefore attempt to 
explain its meaning in detail. 



1. It is an existence and uniqueness theorem. This means that it is a theorem which 
tells us that under certain conditions (stated in the hypothesis) something exists (the 
solution described in the conclusion) and is unique (there is only one such solution). It 
gives no hint whatsoever concerning how to find this solution but merely tells us that 
the problem has a solution. 

2. The hypothesis tells us what conditions are required of the quantities involved. It 
deals with two objects: the differential equation (1.28) and the point (x 0 , y 0 ). As far as 
the differential equation^ 1.28) is concerned, the hypothesis requires that both the 
function f and the function df/dy (obtained by differentiating / (x, y) partially with 
respect to y) must be continuous in some domain D of the xy plane. As far as the point 
(x 0 , y 0 ) is concerned, it must be a point in this same domain D, where / and df/dy are so 
well behaved (that is, continuous). 

3. The conclusion tells us of what we can be assured when the stated hypothesis is 
satisfied. It tells us that we are assured that there exists one and only one solution 0 of 
the differential equation, which is defined on some interval | x — x 0 1 <, h centered about 
x 0 and which assumes the value y 0 when x takes on the value x 0 . That is, it tells us that, 
under the given hypothesis on f(x, y), the initial-value problem 



dx 



= /(*> y)> 



y(x o) = y 0 > 

has a unique solution that is valid in some interval about the initial point x 0 . 

4. The proof of this theorem is omitted. It is proved under somewhat less restrictive 
hypotheses in Chapter 10. 

5. The value of an existence theorem may be worth a bit of attention. What good is 
it, one might ask, if it does not tell us how to obtain the solution? The answer to this 
question is quite simple: an existence theorem will assure us that there is a solution to^ 
look for! It would be rather pointless to spend time, energy, and even money in trying to 
find a solution when there was actually no solution to be found! ; .s for the value of the 
uniqueness, it would be equally pointless to waste time and energy finding one 
particular solution only to learn later that there were others and that the one found was 
not the one wanted! . 
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We have included this rather lengthy discussion in the hope that the student, who has 
probably never before encountered a theorem of this type, will obtain a clearer idea of 
what this important theorem really means. We further hope that this discussion will 
help him to analyze theorems which he will encounter in the future, both in this book 
and elsewhere. We now consider two simple examples which illustrate Theorem LI. 



► Example 1 .1 5 

Consider the initial-value problem 



y( 0 = 3. 

Let us apply Theorem 1*1. We first check the hypothesis. Here /(x, y) — x 2 + y 1 

and — 2 y. Both of the functions / and df/dy are continuous in every domain 

D of the xy plane. The initial condition y(l) — 3 means that x 0 = 1 and y 0 — 3, and the 
point (1, 3) certainly lies in some such domain D. Thus all hypotheses are satisfied and 
the conclusion holds. That is, there is a unique solution (j) of the differential equation 
dy/dx = x 2 4- y 2 , defined on some interval |x — 1 1 < h about x 0 — 1, which satisfies 
that initial condition, that is, which is such that <f>(\) - 3. 



► Example 1.16 



Consider the two problems: 

y( 1 ) = 2 , 



'■ E- 



y 

j? 



dy 

dx 

Here 



_y_ 

J? 



= -^=, y( 0) - 2. 






y 



and 



d/(x, y ) . 

dy 



- 1 / 2 * 



These functions are both continuous except for x = 0 (that is, along the y axis). In 
problem 1, x 0 = 1, y 0 — 2. The square of side 1 centered about (1,2) does not contain 
the y axis, and so both / and df/dy satisfy the required hypotheses in this square. Its 
interior may thus be taken to be the domain D of Theorem 1.1; and (1,2) certainly lies 
within it. Thus the conclusion of Theorem 1.1 applies to problem 1 and we know the 
problem has a unique solution defined in some sufficiently small interval about x 0 — 1. 

Now let us turn to problem 2. Here x 0 - 0,y 0 - 2. At this point'neither / nor df/dy 
are continuous. In other words, the point (0, 2) cannot be included in a domain D where 
the required hypotheses are satisfied. Thus we can nor conclude from Theorem 1.1 that 
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problem 2 has a solution. We are not saying that it does not have one. Theorem 1.1 
simply gives no information one way or the other. 



Exercises 

1. Show that 

y = 4e 2 * + 2e^ x 

is a solution of the initial-value problem 

y{ 0 ) - 6 , 

/( 0 ) - 2 . 

Is y = 2e 2 * + 4e~ 3 * also a solution of this problem? Explain why or why not. 

2. Given that every solution of 

dy 



dx 



+ y = 2xe ” 



may be written in the form y = {x 2 + cje x , for some choice of the arbitrary 
constant c, solve the following initial-value problems: 



(.) g + 

>-( 0 ) = 2 , 

3. Given that every solution of 



dy 

( b ) ^ + y = 2xe *, 
y(- l) = e + 3. 






may be written in the form 



y = c x e* x + c 2 e '‘ 3x , 



for some choice of the arbitrary constants c x and c 2 , solve the following initjal- 
‘ value problems: 






y(0)= -2, 

m = 6. 



m - 5, 

/(0) = 6. 

4. Every solution of the differential equation 




1.3 
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may be written in the form y = c, sin x + c 2 cos x, for some choice of the arbitrary 
constants c t and c 2 - Using this information, show that boundary problems (a) and 
(b) possess solutions but that (c) does 'not. 



d^ +y = °' 

m - 0 , 

y(n! 2) =4 1. 
d 2 y 

( C ) + y = 



S T+>' = 0 . 

y(0) = i, 
y'M=-i. 



y(0) = 0, 
y(n) = I. 

5. Given that every solution of 



X A, 3x^ 

dx 3 dx 2 



dv 

+ 6x -p- - 6y *= 0 
dx 



may be written in the form y ~ CjX + c 2 x 2 + c 3 x 3 for some choice of the arbitrary 
constants c l ,c 2 , and c 3 , solve the initial-value problem consisting of the above 
differential equation plus the three conditions 

y{ 2) = o, y'( 2) = 2, y"(2) - 6. 

6. Apply Theorem 1.1 to show that each of the following initial-value problems has 
a unique solution defined on some sufficiently small inverval |x - 1| < h about 
x 0 = 1: 



(a) ~ = x 2 sin y, 

y(D- —2. 

7. Consider the initial- value problem 



dy = y z 
dx x - 2’ 

y(l) = 0. 



^ = P(x)y 2 + eWy, 

y(2) = 5, 

where P(x) and Q(x) are both third-degree polynomials in x. Has this problem a 
unique solution on some interval Jjc — 2| <h about x 0 = 2? Explain why or why 
not. 

8. On page 19 we stated that the initial-value problem 

dx y ’ 



y( 0) = 0, 



has infinitely many solutions. 
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(a) Verify that this is indeed the case by showing that 

_ £p, * < c, 

y ICfOt-c)] 3 ' 2 , X>c, 

is a solution of the stated problem for every real number c > 0. 

(b) Carefully graph the solution for which c — 0. Then, using this particular 
graph, also graph the solutions for which c — 1, c — 2, and c = 3. 



t 



= CHAPTER TWO 

First-Order Equations for Which Exact Solutions Are Obtainable 



In this chapter we consider certain basic types of first-order equations for which exact 
solutions may be obtained by definite procedures. The purpose of this chapter is to gain 
ability to recognize these various types and to apply the corresponding methods of 
solutions. Of the types considered here, the so-called exact equations considered in 
Section 2.1 are in a sense the most basic, while the separable equations of Section 2.2 are 
in a sense the “easiest.” The most important, from the point of view of applications, are 
the separable equations of Section 2.2 and the linear equations of Section 2.3. The 
remaining types are of various very special forms, and the correponding methods 
of solution involve various devices. In short, we might describe this chapter as a col- 
lection of special “methods,” “devices ” “tricks,” or “recipes,” in descending order of 
kindness! 



2,1 EXACT DIFFERENTIAL EQUATIONS AND INTEGRATING FACTORS 



A, Standard Forms of First-Order Differential Equations 



The first-order differential equations to be studied in this chapter may be expressed in 
either the derivative form 



dx 



= y) 



( 2 . 1 ) 



f 



or the differential form 



M(x, y) dx + N(x, y) dy = 0. 



( 2 . 2 ) 

25 
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An equation in one of these forms may readily be written in the other form. For 
example, the equation 

dy _ x 2 + y 2 
dx x — y 

is of the form (2J). It may be written 

(x 2 + y 2 ) dx + (y - x) dy - 0, 
which is of the form (2.2). The equation 

(sin x + y) dx 4- (x + 3 y) dy - 0, 
which is of the form (2.2), may be written in the form (2.1) as 

dy _ sin x + y 
dx~ x + 3 y ' 

In the form (2.1) it is clear from the notation itself that y is regarded as the dependent 
variable and x as the independent one; but in the form (2.2) we may actually regard 
either variable as the dependent one and the other as the independent. However, in 
this text, in all differential equations of the form (2.2) in x and y, we shall regard y 
as dependent and x as independent, unless the contrary is specifically stated. 



B. Exact Differential Equations 
DEFINITION 



Let F be a function of two real variables such that F has continuous first partial 
derivatives in a domain D. The total differential dF of the function F is defined by the 
formula 



dF(x, y) — 






dF(x> y) 
dy 



dy 



for all (x, y) e D. 



txampie 2A 



Let F be the function of two real variables defined by 

F(x, y) = xy 2 + 2x 3 y 

for all real (x, y). Then 

eF,x ’ ,> 



— 2xy + 2x 3 , 



dx ' " dy 

and the total differential dF is defined by 

dF(x, y) = (y 2 4- 6x 2 y) dx 4- ( 2xy 4- 2x 3 ) dy 



for all real (x, y). 
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DEFINITION 



The expression 



M(x, y) dx + N(x;y) dy 



(2.3) 



is called an exact differential in a domain D if there exists a function F of two real vari- 
ables such that this expression equals the total differential dF(x y y) for all ( x , y) e D. 
That is , expression (2.3) is an exact differential in D if there exists a function F such 
that 



fix 



= M(x, y) and 



dF(x,y) 

dy 



= N(x,y) 



for all (x, y) e D. 

If M(x, y) dx + N(x, y) dy is an exact differential, then the differential equation 
M(x , y) dx 4- N(x, y) dy = 0 
is called an exact differential equation. 



► Example 2.2 
The differential equation 

y 2 dx + 2xy dy - 0 (2.4) 

is an exact differential equation, since the expression y 2 dx + 2xy dy is an exact 
differential. Indeed, it is the total differential of the function F defined for all (x, y) 
by F(x, y) = xy 2 , since the coefficient of dx is dF{x, y)/(dx) = y 2 and that of dy is 
dF(x, y)/(dy) = 2xy. On the other hand, the more simple appearing equation 

y dx + 2x dy = 0, (2.5) 

obtained from (2.4) by dividing through by y, is not exact. 

In Example 2.2 we stated without hesitation that the differential equation (2.4) is 
exact but the differential equation (2.5) is not. In the case of Equation (2.4), we verified 
our assertion by actually exhibiting the function F of which the expression y 2 dx + 

2xy dy is the total differential. But in the case of Equation (2.5), we did not back up 
our statement by showing that there is no fucntion F such that y dx + 2x dy is its total ft 
differential. It is clear that we need a simple test to determine whether or not a given 
differential equation is exact. This is given by the following theorem. 



THEOREM 2.1 

Consider the differential equation 

M(x, y) dx + N(x, y) dy - 0, (2.6) 

where M and N have continuous first partial derivatives at all points (x, y) in a rectangular 
domain D . 
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/. If the differential equation (2,6) is exact in £), then 

dM(x> y) dN(x, y) 



for all (x, y) e D. 

2. Conversely , if 



dy 



dx 



(2.7) 



dM(x, y) dN(x , y) 



dy dx 

for all (x, y) e Z>, then the differential equation (2.6) is exact in D. 



Proof. Part I. If the differential equation (2.6) is exact in D, then M dx + N dy is an 
exact differential in D. By definition of an exact differential, there exists a function F 
such that 

dF(x, & j dF(x ' $ \ 

= M(x, y) and — — — = N(x , y) 



dx 

for all (x, y) e D. Then 

d z F(x,y) dM(x,y) 



dy dx 



dy 



and 



dy 



d 2 F(x,y) dN(x, y) 



dx dy 



dx 



for all (x, y) e D . But, using the continuity of the first partial derivatives of M and /V, we 
have 



and therefore 



d 2 F(x, y) __ d 2 F(x, y) 
dy dx dx dy 

dM(x, y) dN(x, y) 



dy 



dx 



for all (x, y) e D. 

Part 2. This being the converse of Part 1, we start with the hypothesis that 

dM(x,y) dN(x, y) 



dy , 



dx 



for all (x, y) e Z), and set out to show that M dx + N dy - 0 is exact in D. This means 
that we must prove that there exists a function F such that 



dF(x,y) 

dx 



= M(x, y) 



and 



dy 



( 2 . 8 ) 



(2.9) 



for all (x, y) e /). We can certainly find some F(x, y) satisfying either (2.8) or (2.9), but 
what about both? Let us assume that F satisfies (2.8) and proceed. Then 



F(x, >>) = | M(x, 



y)dx + 



(2.10) 
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where j* M(x, y) dx indicates a partial integration with respect to x, holding y constant, 
and 4> is an arbitrary function of y only. This <j){y) is'needed in (2.10) so that F(x t y) 
given by (2 J 0) will represent all solutions of (2.8). It corresponds to a constant of inte- 
gration in the “one- variable” case. Differentiating (2.10) partially with respect to y, 
we obtain 



.V) s - f h/r / 



Now if (2.9) is to be satisfied, we must have 






y) dx + 



y) + 



and hence 



^ = N(x,y)~^M(x,y)dx. 

Since (f> is a function of y only, the derivative d(j>/dy must also be independent of x That 
is, in order for (2.1 1) to hold, 



*’ ^ dy\ M ^ X ’ 



must be independent of x. 
We shall show that 



2[N(x, y )-||M(x.y)5x] = °. 



We at once have 



I ["<*• »> - iy l m - y> Sx ] - - 'ir - - eh \ m - 11 dx ■ 

If (2.8) and (2.9) are to be satisfied, then using the hypothesis (2.7), we must have 



; M(x, 

ft 



, a d 2 F(x,y) d 2 F(x, y) d 2 [ a 

>x " -mT " ~Wsr ■ 37 * J "<*• y) ix - 



Thus we obtain 



• [' V(X * y) ~fy \ M(X * * dX ] = -sVTx \ M(x * y) 5x 



and hence 



|[N(x,y ) -|jM(x,y ) dx] = 



dN(x, y) dM{x,y) 



dM(x, y) __ dN(x, y) 
dy dx 



But by hypothesis (2.7), 
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; 



for all (x, y) e D, Thus 

for ail (x, y) e D , and so (2.12) is independent of x. Thus we may write 

<f>(y) = [ [-V(x, y) - | ax] dy. 

Substituting this into Equation (2.10), we have 

F(x, y) = | M(x, y) Sx + J [^{x, y) - j ax] dy. (2. 1 3) 

This F(x, y) thus satisfies both (2.8) and (2.9) for all (x, y) e D> and so M dx 4- N dy — 0 
is exact in D. Q.E.D. 



Students well versed in the terminology of higher mathematics will recognize that 
Theorem 2.1 may be stated in the following words: A necessary and sufficient condition 
that Equation (2.6) be exact in D is that condition (2.7) hold for all (x, y) e D. For 
students not so well versed, let us emphasize that condition (2.7), 

dM(x,y) _ 8N( Xi y) 
dy dx 

is the criterion for exactness. If (2.7) holds, then (2.6) is exact; if (2.7) does not hold, then 
(2.6) is not exact. 



► Example 2.3 

We apply the exactness criterion (2.7) to Equations (2.4) and (2.5), introduced in 
Example 2.2. For the equation 

y 2 dx + 2 xy-dy = 0 (2.4) 

we have 

M{x, y) = y\ Nix, y) = 2xy, 

dMix, y) _ dNix, y) 

dy y ~ dx 

for all (x, y). Thus Equation (2.4) is exact in every rectangular domain D . On the other 
hand, for the equation 

y dx + 2x dy = 0, (2.5) 

we have 

M(x, y) = y> N(x> y) = 2x, 
dM(x,y) dN (x, y) 

for all (x, y\ Thus Equation (2.5) is not exact in any rectangular domain D. 






j 
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► Example 2.4 

Consider the differential equation 

(2x sin y + y 2 e x ) dx + ( x 2 cos y -f 3y % e x ) dy 0. 



Here 



M(x, y) = 2x sin y + y 3 e x > 
N(x, y) = x 2 cos y + 3y 2 e*, 



ay 



ax 



in every rectangular domain D. Thus this differential equation is exact in every such 
domain. 



These examples illustrate the use of the test given by (2.7) for determining whether or 
not an equation of the form M(x, y) dx + N(x, y) dy = 0 is exact. It should be observed 
that the equation must be in the standard form M(x> y) dx + N(x, y) dy = 0 in order to 
use the exactness test (2.7), Note this carefully: an equation may be encountered in the 
nonstandard form M(x,y)dx = N{x, y)dy f and in this form the test (2.7) does not 
apply. 



C. The Solution of Exact Differential Equations 

Now that we have a test with which to determine exactness, let us proceed to solve exact 
differential equations. If the equation A/(x, y)dx + N(x, y) dy = 0 is exact in a 
rectangular domain D, then there exists a function F such that 

SF ^' ^ = M(x, y) and — = N(x, y) for all (x, y) e D. 

ox dy 

Then the equation may be written 

^ ^ dy -0 or simply dF(x , y) = 0, 

The relation ^(x, y) = cis obviously a solution of this, where c is an arbitrary constant. 
We summarize this observation in the following theorem. 



THEOREM 2-2 

Suppose the differential equation M(x, y) dx 4- N(x, y) dy = 0 satisfies the differentia- 
bility requirements of Theorem 2 A and is exact in a rectangular domain D. Then a one- 
parameter family of solutions of this differential equation is given byF(x, y) = c r where F 
is a function such that 

= and = N(x, y) for all (x, y) £ D. 

ox dy 



and c is an arbitrary constant . 
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Referring to Theorem 2.1, we observe that F(x, y) is given by formula (2.13). 
However, in solving exact differential equations it is neither necessary nor desirable to 
use this formula. Instead one obtains F(x \ y) either by proceeding as in the proof of 
Theorem 2.1, Part 2, or by the so-called ‘‘method of grouping,” which will be explained 
in the following examples. 



► Example 2.5 

Solve the equation 

(3x 2 + 4xy) dx -l- (2x 2 4- 2y) dy = 0. 

Our first duty is to determine whether or not the equation is exact. Here 
M(x, y ) ~ 3x 2 + 4xy, N(x, y) = 2x 2 + 2y, 

SM{x, y) = 4x ^ dNjx, y) = 4x 



dy 



dx 



for all real (x, y), and so the equation is exact in every rectangular domain D, Thus we 
must find F such that 

^ (x ’ , ' ) -M(x,»).3x= + 4x, and - N(x, ,) - 2x- + ly. 



dx 



dy 



From the first of these, 

F(x> y) = J M(x, y) dx 4- 4>(y) = J (3x 2 + 4xy) dx + < p(y ) 
— x 3 + 2x 2 y 4- (j){y ). 



Then 



But we must have 



, My) 



dy 

y ) 

dy 



Thus 



or 



■ - 2x + ■ . 

dy 

« N(x, y) = 2x 2 4- 2y. 

d<t>(y) 



2x 2 4* 2y - 2x 2 4- 



dy 



d<t>(y) 

dy 



= 2 y. * 



Thus <p(y) - y 2 4- c 0 , where c 0 is an arbitrary constant, and so 



F(x, y) = x 3 4- 2x 2 y + y 2 4- c 0 . 
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Hence a one-parameter family of solution is F(x , y) — c l9 or 

x 3 + 2x 2 y 4- y 2 4- c 0 = Cj. 

Combining the constants c 0 and c x we may write this solution as 

x 3 4 - 2x 2 y H- y 2 *= c, 

where c = Cj — c 0 is an arbitrary constant. The student will observe that there is no loss 
in generality by taking c 0 = 0 and writing <j>(y ) = y z . We now consider an alternative 
procedure. 



Method of Grouping. We shall now solve the differential equation of this exam- 
ple by grouping the terms in such a way that its left member appears as the sum of 
certain exact differentials. We write the differential equation 

(3x 2 4- 4xy) dx 4- (2x 2 + 2y) dy = 0 

in the form 



3x 2 dx + (4xy dx + 2x 2 dy) 4* 2 y dy = 0. 

We now recognize this as 

d(x 3 ) + d(2x 2 y) 4* d(y 2 ) = d(c . ), 
where c is an arbitrary constant* or 

d(x 3 + 2x 2 y + y 2 ) — d(c). 

From this we have at once 

x 3 4- 2x 2 y 4 y 2 = c. 

Clearly this procedure is much quicker, but it requires a good ‘‘working knowledge" of 
differentials and a certain amount of ingenuity to determine just how the terms should 
be grouped. The standard method may require more “work" and take longer, but it is 
perfectly straightforward. It is recommended for those who like to follow a pattern and 
for those who have a tendency to jump at conculsions. 

Just to make certain that we have both procedures well in hand, we shall consider an 
initial-value problem involving an exact differential equation. 



► Example 2.6 



Solve the initial-value problem 

(2x cos y 4 - 3x 2 y) dx 4- (x 3 — x 2 sin y — y) dy ~ 0, 

y(0) = 2. 

We first observe that the equation is exact in every rectangular domain D, since 

3M(x,y) ^ . , , 2 5N'x,y) 

— 2x sin y + 3x = — 

oy dx 



for all real (x* y). 
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Standard Method. We must find F such that 

3F(x,y) 



and 



Then 



dx 



W(x,y) 



by 



M(x , y) = 2x cos y + 3 x 2 y 



N(x, y) = x 3 — x 2 sin y — y. 



, y) = | M(x, 
- 1 < 



y) dx + 0 (y) 



(2x cos 7 + 3x 2 y) dx + <£(y) 

= x 2 cos y + x 3 y + <£{y), 

3F(x,y) 2 , 3 , d<p(y) 

dy ' dy 



But also 



dF(x y y) ... 3 2 . 

— - — - = N(x, y) — x — x 2 sin y - y 

and so 

#()>) 
dy y 

and hence 

y 2 

4>(y) = —j + c o- 

Thus 

y ^ 

T(x, y) = x 2 cos y + x 3 y - — + c 0 . 

Hence a one-parameter family of solutions is F(x y y) = c { , which may be expressed as 

y ^ 

x 2 cos y + x 3 y — — = c. 

Applying the initial condition y = 2 when x = 0, we find c — — 2 Thus the solution of 
the given initial-value problem is 

y 2 

x 2 cos y + x 3 y — — ~ — 2. 



Method of Grouping. We group the terms as follows: 

(2v* cos y dx — x 2 sin y dy) + (3x 2 y dx t x 3 dy) -- y dy - 0, 
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Thus we have 

d{x 2 cos y) + d(x 3 y) — = <i(c); 

and so 



x 2 cos y + x 3 y — 

2 

is a one-parameter family of solutions of the differential equation. Of course the initial 
condition y{0) *= 2 again yields the particular solution already obtained. 



D. Integrating Factors 



Given the differential equation 



M (x, y) dx + N{x, y) dy - 0, 
if 

3M(x, y) _ dN( X) y) 
dy dx 

then the equation is exact and we can obtain a one-parameter family of solutions by 
one of the procedures explained above. But if 



dM{x,y) dN(x>y) 
dy dx 

then the equation is not exact and the above procedures do not apply. What shall we do 
in such a case? Perhaps we can multiply the nonexact equation by some expression that 
will transform it into an essentially equivalent exact equation. If so, we can proceed to 
solve the resulting exact equation by one of the above procedures. Let us consider 
again the equation 



y dx + 2 xdy = 0, (2.5) 

which was introduced in Example 2.2. In that example we observed that this equation is 
not exact. However, if we multiply Equation (2.5) by y, it is transformed into the 
essentially equivalent equation 

y 2 dx + 2xy dy = 0, (2.4) 

which is exact (see Example 2.2). Since this resulting exact equation (2.4) is integrable, 
we call y an integrating factor of Equation (2.5). In general, we have the following 
definition: 



DEFINITION 

If the differential equation 

M(x y) dx + N (x, y) dy = 0 (2. 14) 
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| 



i 



is not exact in a domain D but the differential equation 

/i(x, y)M(x , y) dx 4- fi(x 9 y)N(x , y) dy = 0 (2. 1 5) 

is exact in D, then /i(x, y) is called an integrating factor of the differential equation (2. 14 ). 



► Example 2.7 



Consider the differential equation 

(3y + 4xy 2 ) dx + (2x + 3 x 2 y) dy - 0. 

This equation is of the form (2.14), where 

M(x, y) = 3y + 4 xy 2 , N(x , y) = 2x + 3 x 2 y, 



dM(x, y) 



= 3 + 8xy, 



3N(x,y) 



2 + 6xy. 



x y) flV(x, y) 

dy dx 

except for (x, y) such that 2xy -1-1=0, Equation (2.16) is not exact in any rectangular 
domain D. 

Let /j(x, y) - x 2 y. Then the corresponding differential equation of the form (2.15) is 
(3x 2 y 2 + 4x 3 y 3 ) dx -1- (2x 3 y + 3x 4 y 2 ) dy = 0. 

This equation is exact in every rectangular domain Z), since 1 



d[/i(x, y)M(x, y)] 
dy 



* 2 .^32 3 Mx, y)N(x, y)] 

6x 2 y + 12 x*y 2 = 



for all real (x, y). Hence ji(x, y) — x 2 y is an integrating factor of Equation (2.16). 



Multiplication of a nonexact differential equation by an integrating factor thus 
transforms the nonexact equation into an exact one. We have referred to this result- 
ing exact equation as “essentially equivalent 1 * to the Original. This so-called essentially 
equivalent exact equation has the same one-parameter family of solutions as the non- 
exact original. However, the multiplication of the original equation by the integrating 
factor may result in either (1) the loss of (one or more) solutions of the original, or (2) the 
gain of (one or more) functions which are solutions of the “new” equation but not of the 
original, or (3) both of these phenomena. Hence, whenever we transform a nonexact 
equation into an exact one by multiplication by an integrating factor, we should check 
carefully to determine whether any solutions may have been lost or gained. We shall 
illustrate an important special case of these phenomena when we consider separable 
equations in Section 2.2. See also Exercise 22 at the end of this section. 

The question now arises: How is an integrating factor found? We shall not attempt to 
answer this question at this time. Instead we shall proceed to a study of the important 
class of separable equations in Section 2.2 and linear equations in Section 2.3. We shall 
see that separable equations always possess integrating factors that are perfectly 



2.1 
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obvious, while linear equations always have integrating factors of a certain special 
form. We shall return to the question raised above in Section 2.4. Our object here has 
been merely to introduce the concept of an integrating factor. 



Exercises 



In Exercises 1-10 determine whether or not each of the given equations is exact; solve 
those that are exact. 

1. (3x 4 2y) dx 4 (2x 4 y) dy - 0. 

2. ( y 2 4 3) dx 4 (2 xy -4 ) dy = 0. 

3. (2xy 4 1) dx 4 (x 2 4 Ay) dy = 0. 

4. (3x 2 y + 2) dx ~ (x 3 4- y) dy = 0. 

5. (6xy 4 2 y 2 — 5) dx 4 (3x 2 4 4 xy — 6) dy*= 0. 

6. {6 2 4 l)cos r dr 4 26 sin r d8 = 0. 

7. (y sec 2 x 4 sec x tan x) dx 4 (tan x 4 2y) dy = 0. 



8 - (p + dx + (jr + yj d y = o. 



9. 



10 . 



2s -l 



— — dx 4 (3x ,/2 y 1/2 - 1) dy - 0. 



Solve the initial-value problems in Exercises 11-16, 

1 1 . (2xy - 3) dx 4 (x 2 4 Ay) dy = 0, y( 1) = 2. 

12. (3x 2 y 2 - y 3 4 2x) dx 4 (2x 3 y — 3xy 2 4 1) dy = 0, y( — 2) = 1. 

1 3. (2y sin x cos x 4 y 2 sin x) dx 4 (sin 2 x — 2y cos x) dy = 0, y(0) = 3. 

14. ( ye x 4 2e x 4 y 2 ) dx 4 ( e x 4 2xy) dy = 0, y(0) = 6. 



15 , 



16. 






1 4 8xy 2 ' 3 J 2x 4/3 y 2/3 — x 1/3 J „ 
* dx + ^ dy = 0, 



x 2/3 ^ 1/3 



u */ 3 



yd) = 8 - 



17. In each of the following equations determine the constant 4 such that the 
equation is exact, and solve the resulting exact equation: 

(a) (x 2 4 3xy) dx 4 (Ax 2 4 4y) dy = 0. 



(b) {h + p) dx + (^') dy = °- 
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18. In each of the following equations determine the constant A such that the 
equation is exact, and solve the resulting exact equation: 

(a) (Ax 2 y 4- 2y 2 ) dx 4- (x 3 4- 4xy) dy = 0. 

< b) 

19. In each of the following equations determine the most general function N{x, y) 
such that the equation is exact: 

(a) (x 3 4- xy 2 ) dx 4- N(x , y) dy = 0. 

(b) (x “ 2 y~ 2 4- xy ~ 3 ) dx + N(x , y) dy - 0. 

20. In each of the following equations determine the most general function M(x, y) 
such that the equation is exact: 

(a) A/(x, y) dx 4- (2x 2 y* + x*y) dy — 0. 

(b) Af(x, y) dx + (2ye x + y 2 e 2x ) dy - 0. 

21. Consider the differential equation 

(4x -I- 3 y 2 ) dx + 2 xy dy « 0. 

(a) Show that this equation is not exact. 

(b) Find an integrating fa dor of the form x fl , where n is a positive integer. 

(c) Multiply the given equation through by the integrating factor found in (b) 
and solve the resulting exact equation. 

22. Consider the differential equation 

(y 2 4- 2xy) dx - x 2 dy - 0. 

(a) Show that this equation is not exact. 

(b) Multiply the given equation through by y", where n is an integer, and then 
determine n so that y n is an integrating factor of the given equation. 

(c) Multiply the given equation through by the integrating factor found in (b) 
and solve the resulting exact equation. 

(d) Show that y — 0 is a solution of the original nonexact- equation but is not a 
solution of the essentially equivalent exact equation found in step (c). 

(e) Graph several integral curves of the original equation, including all those 
whose equations are (or can be written) in some “special” form. 

23. Consider a differential equation of the form 
[y + xf(x 2 + y 2 )] dx 4- tyf(x 2 + y 2 ) - xj dy = 0. 

(a) Show that an equation of this form is not exact. 

(b) Show that l/(x 2 4- y 2 ) is an integrating factor of an equation of this form. 
Use the result of Exercise 23(b) to solve the equation 

[y 4- x(x 2 4- y 2 ) 2 ] dx + [y(x 2 + y 2 ) 2 - x] dy = 0. 



24 . 



2.2 SEPARABIf EQUATIONS AND EQUATIONS REDUCIBLE TO THIS FORM 



39 


















2.2 SEPARABLE EQUATIONS AND EQUATIONS REDUCIBLE TO THIS FORM 
A. Separable Equations 

DEFINITION 
An equation of the form 



F(x)G(y)dx + f(x)g(y)dy = 0 (2,17) 

is called an equation with variables separable or simply a separable equation. 

For example, the equation (x — 4)y A dx — x 3 (y 2 — 3)dy = 0 is a separable equa- 
tion. 

In general the separable equation (2.17) is not exact, but it possesses an obvious^ 
integrating factor, namely l/f(x)G(y). For if we multiply Equation (2.17) by this 
expression, we separate the variables, reducing (2.17) to the essentially equivalent 
equation 



F(x) 

fix) 



dx + 



9(y) 

G(y) 



dy — 0. 



(2.18) 



This equation is exact, since 

d p*X*)~| = o _ a [g(y)T 

L/wJ dx [g(>>) j 

Denoting F(x)//(x) by Af(x) and g(y)/G{y) by N{y\ Equation (2.18) takes the form 
M(x) dx + N(y) dy = 0. Since M is a function of x only and N is a function of y only, 
we see at once that a one-parameter family of solutions is 

| M (x) dx + | N(y) dy = c, (2.19) 

where c is the arbitrary constant. Thus the problem of finding such a family of solutions 
of the separable equation (2J7) has reduced to that of performing the integrations 
indicated in Equation (2.19). It is in this sense that separable equations are the simplest 
first-order differential equations. 

Since we obtained the separated exact equation (2.18) from the nonexact equation 
(2.17) by multiplying (2.17) by the integrating factor l//(x)G(y), solutions may have 
been lost or gained in this process. We now consider this more carefully. In formally 
multiplying by the integrating factor 1 //(x)G(y), we actually divided by /(x)G(y). We 
did this under the tacit assumption that neither /(x) nor G(_y) is zero; and, under this 
assumption, we proceeded to obtain the one-parameter family of solutions given by 
(2,19). Now, we should investigate the possible loss or gain of solutions that may have 
occurred in this formal process. In particular, regarding y as the dependent variable as 
usual, we consider the situation that occurs if G(y) is zero. Writing the original dif- 
ferential equation (2.17) in the derivative form 



f{x)g(y)^ + F (x)G(y) = 0, 
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we immediately note the following: If y 0 is any real number such that G(y 0 ) — 0, then 
y = y 0 is a (constant) solution of the original differential equation; and this solution 
may (or may not) have been lost in the formal separation process. 

In finding a one-parameter family of solutions a separable equation, we shall 
always make the assumption that any factors by which we divide in the formal 
separation process are not zero. Then we must find the solutions y = y Q of the 
equation G(y) = 0 and determine whether any of these are solutions of the original 
equation which were lost in the formal separation process. 



► Example 2.8 

Solve the equation 

(x - 4)y 4 dx - x 3 (y 2 - 3) dy = 0. 

The equation is separable; separating the variables by dividing by x 3 y 4 , we obtain 

(x-4 )dx (y 2 -l)dy n 

3 A ^ " 



or 

0 x~ 2 — 4x~ 3 ) dx - (y ~ 2 — 3 y^ 4 ) dy - 0. 

Integrating, we have the one-parameter family of solutions 

_I A 1__L_ 

where c is the arbitrary constant. 

In dividing by x 3 y 4 in the separation process, we assumed that x 3 ^ 0 and y 4 ^ 0. 
We now consider the solution y = 0 of y 4 = 0. It is not a member of the one-parameter 
family of solutions which we obtained. However, writing the original differential 
equation of the problem in the derivative form 

dy (jc - 4)j> 4 
dx x 3 (y 2 - 3)’ 

it is obvious that y = 0 is a solution of the original equation. We conclude that it is a 
solution which was lost in the separation process. 



► Example 2.9 

Solve the initial-value problem that consists of the differential equation 

x sin y dx + (x 2 + 1) cos y dy = 0 (2.20) 

and the initial condition 




( 2 . 21 ) 



We first obtain a one-parameter family of solutions of the differential equation (2.20). 



2.2 
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Separating the variables by dividing by (x 2 + 1) sin y, we obtain 



x 

X 2 + ! 



dx + 



cos y 
sin y 



dy = 0. 



Thus 



x dx r cc 

x 2 -f I J si 



cos y 
sin y 



dy - c 0 , 



where c 0 is an arbitrary constant. Recall that 

+ C and |w| = !| 

Then, carrying out the integrations, we find 



u > 0, 
u < 0* 



j In (x 2 + 1) + In j sin y] = c 0 . (2.22) 

We could leave the family of solutions in this form, but we can put it in a. neater form in 
the following way. Since each term of the left member of this equation involves the 
logarithm of a function, it would seem reasonable that something might be accom- 
plished by writing the arbitrary constant c 0 in the form In \c Y \. This we do, obtaining 



} ln(x 2 + l) + ln | sin y) = In |cj. 
Multiplying by 2, we have 

In (x 2 + 1) + 2 In | sin y| = 2 In \c y \. 

Since 

2 In | sin y| — In (sin y) 2 , 

and 

2 In \c\ | - In c\ = In c, 



where 

c = c\ > 0, 

we now have 

In (x 2 + 1) + In sin 2 y — In c. 
Since In A + In B - In AB , this equation may be written 

In (x 2 + l)sin 2 y = In c. 

From this we have at once 

(x 2 + l)sin 2 y — c. 



(2.23) 



Clearly (2.23) is of a neater form than (2.22). 

In dividing by (x 2 + l)sin y in the separation process, we assumed that sin y ^ 0. 
Now consider the solutions of sin y - 0. These are given by y - nn (n = 0, ± 1, 
±2,...). Writing the original differential equation (2.20) in the derivative'form, it is 
clear that each of these solutions y = mz(n — 0, ±1, ±2,. . .) of sin y ~ 0 is a constant 
solution of the original differential equation. Now, each of these constant solutions 
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y =t nn is a member of the one-parameter family (2.23) of solutions of (2.20) for c = 0. 

Thus none of these solutions was lost in the separation process. 

We now apply the initial condition (2.21) to the family of solutions (2.23). We have 

(l 2 + l)sin 2 | = c 

and so c — 2. Therefore the solution of the initial-value problem under consideration is 

(x 2 + l)sin 2 y = 2. 

B. Homogeneous Equations 

We now consider a class of differential equations that can be reduced to separable 
equations by a change of variables. 

DEFINITION 

The first-order differential equation M(x, y) dx + N(x, y) dy = 0 is said to be homo- 
geneous if, when written in the derivative form ( dy/dx ) = f{x, y), there exists a function g 
such that /(x, y) can be expressed in the form g(y/x). 

► Example 2. tO 

The differential equation (x 2 - 3y 2 ) dx + 2xy dy = 0 is homogeneous. To see this, we / 

first write this equation in the derivative form ] 

dy _ 3 y 2 - x 2 
dx 2 xy 

Now observing that 

3y 2 - x 2 __ 3y x _ 3 /y\ 1 / 1 \ 

2xy 2x 2y 2 \x ) 2 \ y/x J 

we see that the differential equation under consideration may be written as 

1/JJ \ 

dx 2 \xj 2 \yfx)' 

in which the right member is of the form g(y/x) for a certain function g . 

► Example 2.11 

The equation 

(y + y/x 2 H- y 2 ) dx — xdy = 0 
is homogeneous. When written in the form 

<iy _ y + x/x 2 + y 2 

Wv V ! 
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the right member may be expressed as 




depending on the sign of x. This is obviously of the form g(y/x ). 

Before proceeding to the actual solution of homogeneous equations we shall 
consider slightly different procedure for recognizing such equations. A function F is 
called homogeneous of degree n if F(fx, ty) = t n F(x , y). This means that if tx and ty are 
substituted for x and y, respectively, in F(x, y), and if t n is then factored out, the other 
factor that remains is the original expression F(x, y) itself. For example, the function F 
given by F(x, y) = x 2 + y 2 is homogeneous of degree 2, since 

F(fx, ty) = (tx) 2 + (ty) 2 = t 2 (x 2 + y 2 ) = t 2 F(x, y). 

Now suppose the functions M and N in the differential equation M(x, y) dx + 
N(x, y)dy = 0 are both homogeneous of the same degree n. Then since M(tx, ty) — 
t n M(x, y), if we let t = 1/x, we have 

"(*■*■ *•’’)" (*)' ^ 

Clearly this may be written more simply as 

and from this we at once obtain 

Likewise, we find 

Now writing the differential equation Af(x, y) dx -1- N(x, y) dy = 0 in the form 

dy _ M (x, y) 
dx N(x, y) ’ 

we find 




Clearly the expression on the right is of the form g(y/x) y and so the equation A/(x, y) 
dx + N(x, y) dy - 0 is homogeneous in the sense of the original definition of homo- 
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geneity. Thus we conclude that if M and N in M(x, y) dx + N(x , y) dy — 0 are both 
homogeneous functions of the same degree n , then the differential equation is a 
homogeneous differential equation. 

Let us now look back at Examples 2.10 and 2.11 in this light. In Example 2.10, 
M(x, y) =*= x 2 — 3y 2 and iV(x, y) = 2 xy. Both M and JV are homogeneous of degree 2. 
Thus we know at once that the equation ( x 2 — 3y 2 ) dx 4- 2 xy dy = 0 is a homogeneous 
equation. In Example 2.1 1, Af(x, y) = y 4 - > jx 1 + y 2 and N(x, y) = -x. Clearly N is 
homogeneous of degree 1. Since 

M(tx, ty) = ty + J{tx) 2 4- (fy) 2 = t(y + ,/x 2 + y 2 ) = t 1 M(x, y), 
we see that M is also homogeneous of degree 1. Thus we conclude that the equation 
(y + Jx 1 + y 2 ) dx — x dy — 0 

is indeed homogeneous. 

We now show that every homogeneous equation can be reduced to a separable 
equation by proving the following theorem. 



THEOREM 2.3 



M(x , y) dx 4- N(x, y) dy = 0 (2.24) 

is a homogeneous equation , f/ie change of variables y — vx transforms (224) into a 

separable equation in the variables v and x. 

Proof. Since Af(x, y) dx 4- N(x y y) dy — 0 is homogeneous, it may be written in the 
form 



Let y = yx. Then 



and (2.24) becomes 



dy dv 

-2- = v + x — 
dx dx 



dv 

v 4- x — = g(v) 



O - g(v)2 dx + x dv = 0. 



This equation is separable. Separating the variables we obtain 



‘—JT + — = 0 . 

» - g(t>) x 



Q.E.D. 
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Thus to solve a homogeneous differential equation of the form (2.24), we let y = vx 
and transform the homogeneous equation into a separable equation of the form (2.25). 
From this, we have 



f dv [dx 

+ J x ~ Ci 



where c is an arbitrary constant Letting F{v) denote 

dv 



1= 



» - g{v) 

and returning to the original dependent variable y, the solution takes the form 



^) + u,w-c. 



► Example 2,12 

Solve the equation 



( x 2 — 3 y 2 ) dx + 2 xy dy = 0. 



We have already observed that this equation is homogeneous. Writing it in the form 



dy = _ * 

dx 2 y lx 



and letting y = vx, we obtain 



or 



or, finally, 



dv 1 3l; 

’ + X T*-- 5 + T- 



dv 1 V 

X d^~ ~S + 2’ 



dv v 2 — ! 






dx 2v 

This equation is separable. Separating the variables, we obtain 

2v dv dx 



Integrating, we find 
and hence 



t> 2 — 1 x 



In |t; 2 — 1| a* In |x| + In | c j , 



|v 2 - 11 = lex), 
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where c is an arbitrary constant. The reader should observe that no solutions were lost 
in the separation process. Now, replacing v by y/x we obtain the solutions in the 
form 




= k*l 



or 



\y 2 — x 2 \ = |cx|x 2 . 



If y > x >: 0, then this may be expressed somewhat more simply as 



y 2 — x 2 — cx 3 . 



► Example 2.1 3 

Solve the initial-value problem 

( y + Jx 1 + y 2 ) dx - x dy = 0, 

yW = 0 . 

We have seen that the differential equation is homogeneous. As before, we write it in the 
form 



dy __ y + J~x 2 + y 2 
dx x 

Since the initial x value is 1, we consider x > 0 and take x — v / x^and obtain 



We let y = vx and obtain 



+ /i + i' ] . 
dx X V \X ; 



dv 



UV f J 

+ X — = U + ./l + V 2 

si v v 



or 



Separating variables, we find 



dx 

dv r t 

dv dx 



Using tables, we perform the required integrations to obtain 



or 



In \v + yjv 2 4* 1| = In |x| -f In |c|, 



v H- yj v 1 + 1 = cx. 
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Now replacing v by y/x , we obtain the general solution of the differential equation in 
the form 



' + *j + l =c* 
x v x* 



y + Jx 2 4 y 1 — cx 2 . 



The initial condition requires that y — 0 when x = \. This gives c - hand hence 



from which it follows that 



y + sfx 1 + y 2 - x 2 . 



y = \{x 2 - l). 



Exercises 

Solve each of the differential equations in Exercises 1—14. 

1. 4 xy dx 4 {x 2 4- 1) dy - 0. 

2. (xy 4 2x 4* y 4- 2) dx 4 (x 2 4- 2x) dy = 0. 

3. 2r(s 2 4- 1) dr 4- (r 4 4- l)ds = 0. 

4. esc y dx 4- sec x dy — 0. 

5. tan 6 dr 4- 2r d6 = 0. 

6. (e v 4 l)cos u du 4- e v (s\n u 4 1) dv - 0. 

7. (x 4 4)(y 2 4 1) dx 4* y(x 2 4- 3x 4- 2) dy = 0. 

8. (x 4 y) dx - x dy = 0. 

9. (2xy 4 3y 2 ) dx - (2 xy 4 x 2 ) dy = 0. 

10. u 3 du 4 (u 3 — m> 2 ) du = 0. 



11. |^x tan - 4 yj dx - x dy = 0. 

12. (2s 2 4- 1st 4- t 2 ) ds 4- (s 2 4- 2sr - t 2 ) dt = 0. 

13. (x 3 + y 2 ^/x 2 4* y 2 ) dx - xy Jx 2 4 y 2 dy *= 0. 

14. (^xTy 4 v/x- y) dx 4 (^/x - y - >/* + y) dy « 0. 
Solve the initial-value problems in Exercises 15-20. 

15. (y 4 2) dx 4 y(x 4 4) dy = 0, y( — 3) = — 1. 



16. 8 cos 2 y dx 4 esc 2 x dy - 0, 



17. (3x 4 8)(y 2 4 4) dx - 4y(x 2 4 5x 4 6) dy = 0, y(l) = 2. 
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18. (x 2 + 3y 2 ) dx - 2xy dy = 0, y(2) = 6. 

19. (2x-5y)dx + (4x- y)dy = 0, y( 1) - 4. 

20. (3x 2 + 9 xy + 5y 2 ) dx - (6x 2 + 4xy) dy ~ 0, y(2) = -6. 

21. (a) Show that the homogeneous equation 

(Ax 4- By) dx + (Cx + Dy) dy - 0 

is exact if and only if B = C. 

(b) Show that the homogeneous equation 

(Ax 2 4- Bxy + Cy 2 ) dx + (Dx 2 + Exy -|- Fy 2 ) dy - 0 

is exact if and only if B — 2D and E = 2 C. 

22. Solve each of the following by two methods (see Exercise 21(a)): 

(a) (x + 2y) dx + (2x — y)dy = 0. 

(b) (3x - y) dx ~ (x + y) dy - 0. 

23. Solve each of the following by two methods (see Exercise 21(b)): 

(a) <x 2 + 2y 2 ) dx + (4xy - y 2 ) dy = 0, 

(b) (2x 2 + 2xy + y 2 ) dx + (x 2 + 2xy) dy = 0. 

24. (a) Prove that if M dx + N dy = 0 is a homogeneous equation, then the change 

of variables x — uy transforms this equation into a separable equation in the 
variables u and x. 

(b) Use the result of (a) to solve the equation of Example 2.12 of the text. 

(c) Use the result of (a) to solve the equation of Example 2.13 of thte text. 

25. Suppose the equation M dx + N dy = 0 is homogeneous. Show that the trans- 
formation x = r cos 0, y = r sin 6 reduces this equation to a separable equation in 
the variables r and 0 . 

26. (a) Use the method of Exercise 25 to solve Exercise 8. 

(b) Use the method of Exercise 25 to solve Exercise 9. 

27. Suppose the equation 

M dx + N dy = 0 (A) 

is homogeneous. 

(a) Show that Equation (A) is invariant under the transformation 

x = fcf, y = tof, (B) 

where Ac is a constant. 

(b) Show that the general solution of Equation (A) can be written in the form 

(c) 

where c is an arbitrary constant. 
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(c) Use the result of (b) to show that the solution (C) is also invariant under the 
transformation (B). 

(d) Interpret geometrically the results proved in (a) and (c). 



2.3 LINEAR EQUATIONS AND BERNOULLI EQUATIONS 
A. Linear Equations 

In Chapter 1 we gave the definition of the linear ordinary differential equation of order 
n ; we now consider the linear ordinary differential equation of the first order. 



DEFINITION 

A first-order ordinary differential equation is linear in the dependent variable y and the 
independent variable x if it is, or can be, written in the form 



~ + P(x)y = G(jc). 



For example, the equation 



+ (x + \)y = x 3 

is a first-order linear differential equation, for it can be written as 

I + (' + 

which is of the form (2.26) with P(x) — 1 + (1/x) and Q(x) = x*. 
Let us write Equation (2.26) in the form 

[ P(x)y - Q(x)] dx + dy - 0. 

Equation (2.27) is of the form 

M(x, y)dx + N (x, y) dy = 0, 



M(x,y) = P(x)y-g(x) and N(x,y)= 1. 



and "M-OL 
dy dx 

Equation (2.27) is not exact unless P(x) = 0, in which case Equation (2.26) degenerates 
into a simple separable equation. However, Equation (2.27) possesses an integrating 
factor that depends on x only and may easily be found. Let us proceed to find it. Let us 
multiply Equation (2.27) by p(x), obtaining 

\_p(x)P(x)y - M*)£2(x)] dx + M*) dy = 0. (2.28) 
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By definition, fi(x) is an integrating factor of Equation (2.28) if and only if Equation 
(2.28) is exact; that is, if and only if 

J- Mx)p(x)y - ^ Mx)l 

This condition reduces to 

M(x)PW = ^ [Mx)]. (2.29) 



In (2.29), P is a known function of the independent variable x, but \i is an unknown 
function of x that we are trying to determine. Thus we write (2.29) as the differential 
equation 



du 

^ x)= f x > 



in the dependent variable // and the independent variable x, where P is a known 
function of x. This differential equation is separable; separating the variables, we have 

du 

— = P(x) dx. 

P 

Integrating, we obtain the particular solution 

tn|/j| =|p(x)iJx 
or 

li « e* Pix) d \ (2.30) 

where it is clear that fi > 0. Thus the linear equation (2.26) possesses an integrating 
factor of the form (2.30). Multiplying (2.26) by (2.30) gives 

e iP( x )dx^_ e sr{ X )dx p{ x )y — Q(x)e iP(x)dx y 

which is precisely 

— [e J/,(x)dx y] rs Q(x)e iPix>dx . { 

Integrating this we obtain the solution of Equation (2.26) in the form 

e \p{x)dxy j* e l p W d *Q(x) dx + c, 

where c is an arbitrary constant. 

Summarizing this discussion, we have the following theorem: 



THEOREM 2.4 

The linear differential equation 



dy 
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has an integrating factor of the form 

€ SPix)dx 

A one-parameter family of solutions of this equation is 

y e !P(x)dx = J e fP(x)d*Q^ fa + 

that is , 

y = e~ l J e iF{x)d *Q(x) dx 4- cj. 

Furthermore , it can be shown that this one-parameter family of solutions of the linear 
equation (2.26) includes all solutions of (2.26). 

We consider several examples. 



► Example 2.14 



Here 



dx 






<T 2 *. 



(2.31) 



P(*) = 



2x + 1 



and hence an integrating factor is 



expj^J P(x) dxj «= expj^ J ^ = exp(2x + In |*|) 

= exp(2x) exp(ln |x|) = x exp(2x).* 

Multiplying Equation (2.31) through by this integrating factor, we obtain 

xe 2x + <? 2x (2x + 1 )y = x 



or 



dx 

Integrating, we obtain the solutions 



(x e 2x y) = x. 



xe ix y - — + c 



or 



y = \xe 2x + — e 2x , 



where c is an arbitrary constant. 



The expressions e* and exp x are identical. 
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► Example 2*1 5 

Solve the initial- value problem that consists of the differential equation 



and the initial condition 



(x 2 + 1) ^ + Axy = x 



3 * 2 ) - 1 . 



The differential equation (2.32) is not in the form (2.26). We therefore divide by x 2 + 1 
to obtain 

dy , Ax x 

— + — 7 y = -= 7. (2.34) 

dx X 2 -I- 1 x 2 + 1 

Equation (2.34) is in the standard form (2.26), where 



x 2 + r 



An integrating factor is 



cxpf" J P(x) dx 1 = expQ —~) = exp[ln(x 2 + l) 2 ] = (x 2 + l) 2 
Multiplying Equation (2.34) through by this integrating factor, we have 
(x 2 + !) 2 + 4x(x 2 4- l)y = x(x 2 + 1) 



£ K* a + i) 2 y3 = x 3 + x. 

We now integrate to obtain a one-parameter family of solutions of Equation (2.23) in 
the form 

(x 2 + l) 2 y — + ~2 + c. 

Applying the initial condition (2.33), we have 

25 = 6 + c. 

Thus c - 19 and the solution of the initial-value problem under consideration is 

x 2 

(x 2 + \) 2 y = — + — + 19. 



► Example 2.16 

Consider the differential equation 

y 1 dx + (3 xy - 1) dy — 0. 
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Solving for dy/dx> this becomes 

dy = ,v 2 

dx 1-3 xy 9 

which is clearly not linear in y, Also, Equation (2.35) is not exact, separable, or 
homogeneous. It appears to be of a type that we have not yet encountered; but let us 
look a little closer. In Section 21, we pointed out that in the differential form of a first- 
order differential equation the roles of x and y are interchangeable, in the sense that 
either variable may be regarded as the dependent variable and the other as the 
independent variable. Considering differential equation (2.35) with this in mind, let us 
now regard x as the dependent variable and y as the independent variable. With this 
interpretation, we now write (2.35) in the derivative form 

dx 1 — 3 xy 
dy ~ y 2 



or 



dx 3 1 

dy y y 



Now observe that Equation (2,36) is of the form 



dx 

dy 



+ ?(y)x = Q(y) 



(2.36) 



and so is linear in x. Thus the theory developed in this section may be applied to 
Equation, (2.36) merely by interchanging the roles played by x and y. Thus an inte- 
grating factor is 

exp|^| P(y) dyj = exp^j* ~ dy^ = exp(ln | y) 3 ) = y 3 . 



Multiplying (2.36) by y 3 we obtain 

. dx , , 

> 3f +3 ’ x ~f 



£ Cy,xl - y - 

Integrating, we find the solutions in the form 




x = 



1 c 

2 y + ?’ 



where c is an arbitrary constant. 
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B. Bernoulli Equations 

We now consider a rather special type of equation that can be reduced to a linear 
equation by an appropriate transformation. This is the so-called Bernoulli equation. 

DEFINITION 



An equation of the form 



£ + P{x)y = Q(x)y* (2.37) 

is called a Bernoulli differential equation. 

We observe that if n = 0 or 1, then the Bernoulli equation (2.37) is actually a linear 
equation and is therefore readily solvable as such. However, in the general case in which 
n ^ 0 or l, this simple situation does not hold and we must proceed in a different 
manner. We now state and prove Theorem 2.5, which gives a method of solution in the 
general case. 

THEOREM 2.5 

Suppose n ^Oor 1. Then the transformation v — y l reduces the Bernoulli equation 

^ + P(x)y = Q(x)y” (2.37) 

to a linear equation in v. 

Proof. We first multiply Equation (2.37) by y~\ thereby expressing it in the 
equivalent form 

y- n ~ + p (x)y'- n = Q(x). (2.38) 



If we let v = y 1 "'*, then 



dv x dy 



and Equation (238) transforms into 



1 — n dx 



+ P(x)» = Q(x) 



or, equivalently, 



+ (1 - n)P(x)v - (1 - n)Q(x). 



Letting 



-n)P{*) 
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and 



this may be written 



which is linear in u. 



CM*) = (1 - n)Q(x), 

dv 



dx 



+ AM® = GiM* 



Q.£.D. 



► Example 2.1 7 



dy i 



(2.39) 



This is a Bernoulli differential equation, where n = 3. We first multiply the equation 
through by y -3 , thereby expressing it in the equivalent form 

dy 

y -r- + y =* x. 

dx 

If we let v — y 1 "" - y“ 2 , then dv/dx = — 2y' 3 (dy/dx) and the preceding differential 
equation transforms into the linear equation 

1 dv 

-iTx + ” x - 



Writing this linear equation in the standard form 

2» = — 2x, 
dx 

we see that an integrating factor for this equation is 

e s P{x)dx _ e ~f 2 dx __ e - 

Multiplying (2.40) by e ~ 2x , we find 

. dv 



(2.40) 



„~2x 



dx 



or 



Integrating, we find 



-r~{e 2x v) = — 2xe 2 *. 
dx 

e~ 2x v — je“ 2 *(2x + 1) + c, 
u = x -f- £ -f ce 2 *. 



where c is an arbitrary constant. But 



v = ■ 
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Thus we obtain the solutions of (2.39) in the form 

-4- = * 4- { + ce 2x . 

y 

Note . Consider the equation 

~~ + mm = ew, (2.41) 

where / is a known function of y. Letting v — f(y ), we have 

dv dv dy df(y) dy 

dx dy dx dy dx 9 

and Equation (2.41) becomes 

^ + P(x)v = Q(x), 

which is linear in v. We now observe that the Bernoulli differential equation (2.37) is a 
special case of Equation (2.41). Writing (2.37) in the form 

y-"^ + P(x)y l -” = Q(x) 

and then multiplying through by (1 - n), we have 

^~n)y-"^ c + P l (x)y i -^Q l (x), 

where P x (x) = (i - n)F(x) and Q x {x) — (1 - n)Q(x). This is of the form (2.41), where 
f(y) = y 1 ^; letting v =* y 1 “ n , it becomes 

£ + Pt(x)» = Qi(x), 

which is linear in v. For other special cases of (2.41), see Exercise 37. 



Exercises 



Solve the given differential equations in Exercises 1-18. 



1 . + 

dx x 

3. ^ + 3y = 3x 2 e~ 3x . 

dx 

- dx x 1 

5. — + ^ 

dt V t 2 

dy 2x + 1 

x 1 L + —— r y~x-l. 

dx x + 1 



2. x 4 ^ + 2x 3 y=l. 
dx 

4. ~ + 4xy — 8x. 
dx 

6. (u 2 4- i)~~ 4 4uu - 3u. 
du 



1 . 
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8 . 



(x 2 + x-2)^ + 3(x+ l)y = x- 1. 
dx 



9. 

10 . 

11 . 

12 . 

13. 

14. 

15. 
17. 



x dy + (xy ± y - \) dx = 0. 
y dx 4- ( xy 2 + x — y) dy = 0. 
dr 

— + r tan 6 = cos 6. 
uu 



cos 6 dr + (r sin 0 — cos 4 0) dS = 0. 

(cos 2 x — y cos x) dx — (1 + sin x) dy = 0. 
(y sin 2x — cos x) dx + (1 + sin 2 x) dy = 0. 

‘h_y_ = _y^ 16 

dx x x ' 

dy + (4y-8y~ 3 )xdx = 0. 18. 



dx t + 1 ( + 1 

-3 — | — - — x — . 

dt 2 t xt 



Solve the initial-value problems in Exercises 19-30. 

19. x~ — 2y — 2x 4 > y(2) = 8. 

ax 

d v 

20. — - + 3x 2 y - x 2 , y(0) - 2. 

/IV 



21 . 

22 . 

23. 

24. 

25. 

26. 



e x ly - 3(e x + l) 2 ] dx + (e x + 1) dy = 0, 
2x(y + 1) dx - (x 2 + 1) dy — 0, y(l) = 



+ r tan 6 = cos 2 0. 
du 






dx 

— x = sin It, 

dt 



x(0) - 0. 



x~+y = (xy) 3 ' 2 , y{l) = 4. 



j>(0) - 4. 
-5. 



27. ~ + y ~ f(x), where /(x) = | 2, ® ^ * < y(0) = 0. 

2S. j; + y = f(x), where /(*)-{* = 6 - 

29 tx + y = f{x) ' Where m = \ e e \ ll X 2 <X ^ = '• 

30. (x + 2) + y = f(x), where f(x) = j 2 *’ ^ * < 2 ’ y(0) = 4. 
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31. Consider the equation a(dyjdx) + by = ke where a, b, and k are positive 
constants and A is a nonnegative constant. 

(a) Solve this equation. 

(b) Show that if X = 0 every solution approaches kjb as x -+ go, but if X > 0 
every solution approaches 0 as x -» oo. 

32. Consider the differential equation 

£+«*>- ° 

(a) Show that if / and g are two solutions of this equation and c i and c 2 are 
arbitrary constants, then c t / -f c 2 g is also a solution of this equation. 

(b) Extending the result of (a), show that if /i , f% are n solutions of this 

equation and c lf c 2 ,. ..,c„ are n arbitrary constants, then 

Y c kfk 

I *=1 

is also a solution of this equation.. 

33. Consider the differential equation 

~ + /’W3' = 0, (A) 



where P is continuous on a real interval /. 

(a) Show that the function / such that f(x) — 0 for all x e 1 is a solution of this 
equation. 

(b) Show that if / is a solution of (A) such that /(x 0 ) = 0 for some x 0 e 1 , then 
/(x) = Ofor all x e /. 

(c) Show that if / and g are two solutions of (A) such that / (x 0 ) = ^(x 0 ) for some 
x 0 e /, then f(x) — g(x) for all x e /. 

34. (a) Prove that if / and g are two different solutions of 

~ + my = G(x), (A) 

then / — g is a solution of the equation 

E + '■<*»-»• 

(b) Thus show that if / and g are two different solutions of Equation (A) and c is 
an arbitrary constant, then 

c(f ~g) + f 

is a one-parameter family of solutions of (A). 



35. (a) Let f x be a solution of 



+ P(x).K = Qy(x) 
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and f 2 be a solution of 



— + P(x)y = Q 2 {x\ 



where P, Q u and Q 2 are all defined on the same real interval /. Prove that 
f\ + fi is a solution of 

$ + P(x)y = Q 1 (x) + e 2 (x) 



(b) Use the result of (a) to solve the equation 

dy 

— -f v 2 sin x + 5 sin 2x. 
dx 

36. (a) Extend the result of Exercise 35(a) to cover the case of the equation 

% + PM* = t &<*), 

ax k=i 

where P, Q k {k = 1, 2, . . . , n) are all defined on the same real interval /. 
(b) Use the result obtained in (a) to solve the equation 

dy A 

/ + y - t sin kx . 

dx k= i 

37. Solve each of the following equations of the form (2.41): 

dy 1 

(a) cos y — — + — sin y = 1. 

dx x 

( b ) (y + *) % + x( ~y 2 + 2 y ) = *• 



38. The equation 



^ = A(x)y' l + B(x)y + C(x) 



is called Riccati's equation . 

(a) Show that if A(x) = 0 for all x> then Equation (A) is a linear equation, 
whereas if C(x) = 0 for all x > then Equation (A) is a Bernoulli equation. 

(b) Show that if / is any solution of Equation (A), then the transformation 

reduces (A) to a linear equation in v. 

In each of Exercises 39-41, use the result of Exercise 38(b) and the given solution to 
find a one-parameter family of solutions of the given Riccati equation: 



39. — = (1 - x)y 2 + (2x - l)y - x; given solution f{x) - 1. 
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y 2 4 xy 4 1; given solution f(x) ~ x. 

8xy 2 4 4x(4x 4 l)y — (8x 3 4- 4x 2 — 1); given solution f(x) = x. 

Exercises: Miscellaneous Review 

Solve each of the differential equations in Exercises 1-14. Several can be solved by at 
least two different methods. 

1. 6x 2 y dx — (x 3 + 1) dy = 0. 

2. (3x 2 y 2 — x) dy 4- (2xy 3 — y) dx = 0. 

3. (y — 1) dx 4- x(x 4-1 )dy - 0. 

4. (x 2 - 2y) dx — x dy = 0. 

5. (3x - 5 y) dx 4 (x 4- y) dy - 0. 

6. e 2x y 2 dx 4- {e 2 *y — 2y) dy = 0. 

7. (8x 3 y — 12x 3 ) dx 4 (x 4 4 1) dy = 0. 

8. (2x 2 4 xy 4 y 2 ) dx 4- 2x 2 dy ~ 0. 

dy _ 4x 3 y 2 - 3x 2 y 
dx x 3 - 2x 4 y 

10. (x 4 1 ) 4* xy — 

ax 

dy = 2x - 7y 
dx 3 y - 8x ‘ 

V 

12. x 2 ^ + xy = xy 3 . 

13. (x 3 + 1) ^- + 6x 2 y = 6x 2 . 

dx 

i4 2x2 + yl 

dx 2 xy — x 2 ' 

Solve the initial value problems in Exercises 15-24. 

15. (x 2 4- y 1 ) dx - 2xy dy = 0, y(l) = 2. 

16. 2(y 2 4 4) dx 4- (1 - x 2 )y dy = 0, y(3) = 0. 

17. (e 2x y 2 — 2x) dx + e 2x y dy = 0, y(0) — 2. 

18. (3x 2 4- 2xy 2 ) dx 4- (2x 2 y 4 6y 2 ) dy = 0, y(l) = 2. 

19. 4xy -” = y 2 4l, y(2) = 1. 



40. 

dx 



41. 
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20. 








21. 








22. 


^ + y = fix), where f(x) = 


fl, 0 < x < 2, 
jo, X > 2, 


y(0) = o. 


23. 


(x + 2) ^ + y = /(x), where 
dx 


Vf A 
o * 

H 

CM 

II 

TT 


x<,2, 

2 y(0) = 4. 


24. 









2.4 SPECIAL INTEGRATING FACTORS AND TRANSFORMATIONS 

We have thus far encountered five distinct types of first-order equations for which 
solutions may be obtained by exact methods, namely, exact, separable, homogeneous, 
linear, and Bernoulli equations. In the case of exact equations, we follow a definite 
procedure to directly obtain solutions. For the other four types definite procedures for 
solution are also available, but in these cases the procedures are actually not quite so 
direct In the cases of both separable and linear equations we actually multiply by 
appropriate integrating factors that reduce the given equations to equations that are of 
the more basic exact type. For both homogeneous and Bernoulli equations we make 
appropriate transformations that reduce such equations to equations that are of the 
more basic separable and linear types, respectively. 

This suggests two general plans of attack to be used in solving a differential equation 
that is not of one of the five r types mentioned. Either (1) we might multiply the given 
equation by an appropriate integrating factor and directly reduce it to an exact 
equation, or (2) we might make an appropriate transformation that will reduce the 
given equation to an equation of some more basic type (say, one of the five types 
already studied). Unfortunately no general directions can be given for finding an 
appropriate integrating factor or transformation in all cases. However, there is a 
variety of special types of equations that either possess special types of integrating 
factors or to which special transformations may be applied. We shall consider a few of 
these in this section. Since these types are relatively unimportant, in most cases we shall 
simply state the relevant theorem and leave the proof to the exercises. 



A. Finding Integrating factors 

The so-called separable equations considered in Section 2.2 always possess integrating 
factors that may be determined by immediate inspection. While it is true that some 
nonseparable equations also possess integrating factors that may be determined “by 
inspection,” such equations are rarely encountered except in differential equations texts 
on pages devoted to an exposition of this dubious “method.” Even then a considerable 
amount of knowledge and skill are often required. 
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Let us attempt to attack the problem more systematically. Suppose the equation 
M(x, y) dx 4- N(x , y) dy = 0 (2.42) 

is not exact and that ^(x, y) is an integrating factor of it. Then the equation 

/i(x, y)M(x, y) dx + j*(*> y)A r (x, y) dy = 0 (2.43) 

is exact. Now using the criterion (2.7) for exactness, Equation (2.43) is exact if and only if 

j- y y) M ( x > y)] = ^ CM*. yW(x> y)J- 

This condition reduces to 



N(x,y) 



Qfijx, y) 
dx 



-M{x,y) 



fy(x, y ) 

dy 



\dM(x,y) dN(x,y)~\ 

L~5 5^— ^ 



Here M and N are known functions of x and y, but pi is an unknown function of x and y 
that we are trying to determine. Thus we write the preceding condition in the form 



< i44 > 

Hence fx is an integrating factor of the differential equation (2.42) if and only if it is a 
solution of the differential equation (2.44). Equation (2.44) is a partial differential 
equation for the general integrating factor pt, and we are in no position to attempt to 
solve such an equation. Let us instead attempt to determine integrating factors of 
certain special types. But what special types might we consider? Let us recall that the 
linear differential equation 

■j- + P(*>* = QW 



always possesses the integrating factor e /FU)d *, which depends only upon x. Perhaps 
other equations also have integrating factors that depend only upon x. We therefore 
multiply Equation (2.42) by pt(x), where pi depends upon x alone. We obtain 



/i(x)M(x, y) dx + MW# (x* y) dy = 0. 

This is exact if and only if 

lfi(x)M{x, y)] = [fi(x)N(x, y)]. 

Now M and N are known functions of both x and y, but here the intergrating factor pt 
depends only upon x. Thus the above condition reduces to 






dy 



dx 



+ N(x, y) 



Mx) 

dx 



or 



Mx) 1 |~ 8M(x, y) 

n(x) N(x, y) |_ dy 



dN(x, y)~ 



dx. 



dx 



(2-45) 
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If 



1 pM(xj) dN(x,y) 1 

M(x> y) |_ dy dx J 

involves the variable y, this equation then involves two dependent variables and we 
again have difficulties. However, if 

1 [ dM(x, y) _ dN(x y y) 1 

N(x, y) [ dy dx ] 

depends upon x only. Equation (2,45) is a separated ordinary equation in the single 
independent variable x and the single dependent variable p. In this case we may 
integrate to obtain the integrating factor 

In like manner, if 

1 pN(x,y) gM(x, y n 
M(x, y) \_ dx dy J 

depends upon y only, then we may obtain an integrating factor that depends only on y. 
We summarize these observations in the following theorem. 



THEOREM 2.6 



Consider the differential equation 

M(x, y) dx + N{x, y) dy = 0. 

If 

I f £M (x. y) 8N(x, y) 

N{x, y) j_ dy dx 

depends upon x only, then 



] 



_ n ff 1 \cM(x,y) dN(x, y)"| 1 

e * p iJ^(3r^L— ^ 



is an integrating factor of Equation {2.42), If 



1 



I" dN{x , y) _ dM{x , y) ~ 
dx dy 



depends upon y on/y, then 



f f 1 pN(x.y) dAf(x,y)l | 

*r\ y i 



is an integrating factor of Equation (2.42). 



(2.42) 

(2.46) 

(2.47) 

(2.48) 
(2:49) 



We emphasize that, given a differential equation, we have no assurance in general 
that either of these procedures will apply. It may well turn out that (2.46) involves y 
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and (2,48) involves x for the differential equation under consideration. Then we must 
seek other procedures. However, since the calculation of the expressions (2.46) and 
(2.48) is generally quite simple, it is often worthwhile to calculate them before trying 
something more complicated. 



► Example 2.18 



Consider the differential equation 

(2x 2 + y) dx + (x 2 y — x)dy = 0. (2.50) 

Let us first observe that this equation is not exact, separable, homogeneous, linear, 
or Bernoulli. Let us then see if Theorem 2.6 applies. Here A/(x, y) — 2x 2 + y, and 
N(x , y) - x 2 y - x, and the expression (2.46) becomes 



x 2 y ■ 



[1 — (2xy — 1)] = 



2( 1 - xy) 
x(xy- 1) 



This depends upon x only, and so 



exp 






exp( — 2 In | x | ) = 



is an integrating factor of Equation (2.50). Multiplying (2.50) by this integrating factor, 
we obtain the equation 



(2 + p)* + (y-j)^ = o- 



(2.51) 



The student may readily verify that Equation (2.5 1 ) is indeed exact and that the solution 
is 



2x + 



x 



c. 



More and more specialized results concerning particular types of integrating factors 
corresponding to particular types of equations are known. However, instead of going 
into such special cases we shall now proceed to investigate certain useful 
transformations. 



B. A Special Transformation 

We have already made use of transformations in reducing both homogeneous and 
Bernoulli equations to more tractable types. Another type of equation that can be 
reduced to a more basic type by means of a suitable transformation is an equation of 
the form 

(a { x +b Y y + cj dx + {a 2 x -f b 2 y + c 2 ) dy - 0. 

We state the following theorem concerning this equation. 



2,4 SPECIAL INTEGRATING FACTORS AND TRANSFORMATIONS 65 



THEOREM 2.7 

Consider the equation 

(a x x + b x y + c y ) dx 4- (a 2 x + b 2 y + c 2 ) dy = 0, (2.52) 

where a u b Xi c u a 2i b 29 and c 2 are constants . 

Case 1. // a 2 fa y ^ b 2 fb yi then the transformation 

x = x + K • 
y « y + fc, 

w/iere (/i, k) is the so/ution 0 / the system 

a x h + b x k + c y = 0, 
a 2 h + b 2 k + c 2 = 0, 
reduces Equation (2.52) to the homogeneous equation 

(a l X-¥b l Y)dX^(a 2 X^b 2 Y)dY = 0 
in the variables X and Y. 

Case 2. If a 2 /a l = = k, then the. frans/ormathw z = a y x + b x y reduces the 

equation (2.52) to a separable equation in the variables x and z. 

Examples 2.19 and 2.20 illustrate the two cases of this theorem. 



► Example 2.19 

(x — 2y+ 1) dx 4- (4* - 3 y ~ 6} dy = 0. 

Here a x - 1, b x - -2,a 2 - 4, b 2 ~ - 3, and so 

a 2 . . b 2 3a 2 

— = 4 but . 

< 2 , 2 ^ 

Therefore this is Case 1 of Theorem 2.7. We make the transformation 

x - X + K 

y=Y + K 

where (h, k) is the solution of the system 

h - 2k + 1 = 0, 

4h — 3/c — 6 — 0. 

The solution of this system is h = 3, k = 2, and so the transformation is 

x = * + 3, 

y = Y 4* 2. 
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This reduces Equation (2.53) to the homogeneous equation 
(X - 2 Y) dX + (4X - 3Y)dY = 0. 



(2.54) 



Now following the procedure in Section 2.2 we first put this homogeneous equation in 
the form 

dY 1-2 (YIX) 



and let Y ~ vX to obtain 



This reduces to 



dX 3 (YIX) - 4 

dv 1 - 2v 
V + X dX~ 3u — 4' 

(3t; — 4) dv 



dX 
X ‘ 



3v 2 - 2v - 1 

Integrating (we recommend the use of tables here), we obtain 

|3w — 3 1 



(155) 



i In |3w 2 — 2f — 1| In 



or 



or 



ln(3ir — 2v — \) 2 - In 



3v + 1 
3v — 3 



- — In | A" | -h In \Ci | 



■3v + 1 



-"#)• 



In 



(3v + l) f 



v — 1 



= In 



or, finally, 



X A \{3v 4- 1) 5 | =a c\v — 1 1, 



where c = cf. These are the solutions of the separable equation (2.55). Now replacing 
v by Y/X, we obtain the solutions of the homogeneous equation (2.54) in the form 

\3Y + X\ 5 =:c\Y -X |. 

Finally, replacing X by x — 3 and Y by y — 2 from the -original- transformation, vve 
obtain the solutions of the differential equation (2,53) in the form 

1 3(> — 2) + (x — 3)( 5 = c|y — 2 — x + 3| 

or 

|x + 3y - 9| 5 = c\y - x + 1|. 



► Example 2.20 



(x + 2y + 3) dx + (2x + 4y — 1) dy = 0. 



(2.56) 
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Here a x » 1, 6j = % a 2 =2,b 2 = 4, and a 2 fa i = b 2 /hi = 2. Therefore, this is Case 2 of 
Theorem 2.7. We therefore let 



z = x + 2y, 

and Equation (2.56) transforms into 

{z + 3)dx + (2z-l)^~j = 0 
or 

7 dx + (2z - 1) dz = 0, 
which is separable. Integrating, we have 

7x + z 2 — z = c. 

Replacing z by x 4- 2y we obtain the solution of Equation (2.56) in the form 
7x + (x 4* 2 y) 1 - (x 4- 2y) = c 
or 

x 2 4- 4 xy 4- 4y 2 + 6x — 2y - c. 



C. Other Special Types and Methods; An Important Reference 

Many other special types of first-order equations exist for which corresponding special 
methods of solution are known. We shall not go into such highly specialized types 
in this book. Instead we refer the reader to the book Differentialgleichungen: 
Losungsmethoden und Losurtgen ; by E. Kamke (Chelsea, New York, 1948). This re- 
markable volume contains discussions of a large number of special types of equations 
and their solutions. We strongly suggest that the reader consult this book whenever he 
encounters an unfamiliar type of equation. Of course one may encounter an equation 
for which no exact method of solution is known. In such a case one must resort to 
various methods of approximation. We shall consider some of these general methods 
in Chapter 8. 



Exercises 

Solve each differential equation in Exercises 1-4 by first finding an integrating factor. 
1. (5xy + 4y 2 + 1) dx 4- (x 2 + 2xy) dy - 0. 

2: (2x + tan y) dx + (x — x 2 tan y) dy = 0. 

3. [y 2 (x + 1) 4- y] dx + (2xy + 1) dy - 0. 

4. (2xy 2 4- y) dx 4- (2y 3 - x) dy = 0. 

In each of Exercises 5 and 6 find an integrating factor of ti.e form x p y q and solve. 

5. (4xy 2 + 6y) dx 4- (5x 2 y + 8x) dy - 0. 

6. (8x 2 y 3 — 2y 4 ) dx 4- (5x 3 y 2 — 8xy 3 ) dy = 0. 
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Solve each differential equation in Exercises 7-10 by making a suitable transformation. 

7. (5x + 2y + \)dx + (2x + y + l)rfy = 0. 

8. (3x — y 4- 1 ) dx — (6x — 2y — 3) dy — 0. 

9. (x — 2y — 3) dx + (2x 4- y — 1) dy = 0. 

10. (lOx - 4 y 4- 12) dx - (x 4- 5y 4- 3) dy - 0. 

Solve the initial-value problems in Exercises 11-14. 

11. (6x + 4y 4~ 1) dx 4- (4x + 2y 4- 2) dy = 0, y(|) - 3. 

12. (3x - y - 6) dx 4- (x 4- y 4- 2)'dy = 0, y(2) - -2. 

13. (2x + 3y + 1) dx + (4x + 6y 4- 1) dy - 0, y(- 2) = 2. 

14. (4x 4- 3y 4- 1) dx + (x + y + 1) dy — 0, y( 3) = — 4. 

15. Prove Theorem 2.6. 

1 6. Prove Theorem 2.7. 

17. Show that if /i(x, y) and t;(x, y) are integrating factors of 

M(x , y) dx 4- N(x, y) rfy = 0 (A) 

such that /i(x, y)/y(x, y) is not constant, then 

/z(x, y) = cv{x , y) 

is a solution of Equation (A) for every constant c. 

18. Show that if the equation 

M(x, y) dx + jV(x, y) dy = 0 (A) 

is homogeneous and M(x, y)x 4- N(x, y)y ^ 0, then l/[M(x, y)x 4- N(x, y)y] is 
an integrating factor of (A). 

19. Show that if the equation M(x, y) dx + N(x, y) dy = 0 is both homogeneous and 
exact and if M(x, y)x 4- N(x, y)y is not a constant, then the solution of this 
equation is M(x, y)x 4- N(x> y)y = c, where c is an arbitrary constant. 

20. An equation that is of the form 

y = px+ f{p), (A) 

where p s dy/dx and / is a given function, is called a Clairaut equation. Given 
such an equation, proceed as follows: 

1. Differentiate (A) with respect to x and simplify to obtain 

[X + f'(p)-] ^ = o. (B) 

Observe that (B) is a first-order differential equation in x and p. 

2. Assume x 4- f'(p) # 0, divide through by this factor, and solve the resulting 
equation to obtain 

p = c, (C) 

where c is an arbitrary constant. 
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3. Eliminate p between (A) and (C) to obtain 

y = cx + f(c). (D) 

Note that (D) is a one-parameter family of solutions of (A) and compare the form 
of differential equation (A) with the form of the family of solutions (D). 

4. Remark Assuming x 4* f\p) — 0 and then eliminating p between (A) and 
x + f{p) = 0 may lead to an “extra 1 ’ solution that is not a member of the one- 
parameter family of solutions of the form (D). Such an extra solution is usually 
called a singular solution. For a specific example, see Exercise 21. 

Consider the Clairaut equation 



y = px + p 2 , where p 



dy 

dx‘ 



(a) Find a one-parameter family of solutions of this equation. 

(b) Proceed as in the Remark of Exercise 20 and find an “extra” solution that is 
not a member of the one-parameter family found in part (a). 

(c) .Graph the integral curves corresponding to several members of the one- 
parameter family of part (a); graph the integral curve corresponding to the 
“extra” solution of part (b); and describe the geometric relationship between 
the graphs of the members of the one-parameter family and the graph of the 
“extra” solution. 



CHAPTER THREE 

Applications of First-Order Equations 



In Chapter 1 We pointed out that differential equations originate from the mathemati- 
cal formulation of a great variety of problems in science and engineering. In this 
chapter we consider problems that give rise to some of the types of first-order ordinary 
differential equations studied in Chapter 2. First, we formulate the problem mathemat- 
ically j thereby obtaining a differential equation. Then we solve the equation and 
attempt to interpret the solution in terms of the quantities involved in the original 
problem. 



3.1 ORTHOGONAL AND OBLIQUE TRAJECTORIES 
A. Orthogonal Trajectories 

DEFINITION 

Let 

F(x, y, c) — 0 (3.1) 

be a given one-parameter family of curves in the xy plane . A curve that intersects the 
curves of the family (5./) at right angles is called an orthogonal trajectory of the given 
family . 



► Example 3.1 
Consider the family of circles 



f 
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(3.2) 
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with center at the origin and radius c. Each straight line through the origin, 

y - kx, (3.3) 

is an orthogonal trajectory of the family of circles (3.2). Conversely, each circle of the 
family (3.2) is an orthogonal trajectory of the family of straight lines (3.3). The families 
(3.2) and (3.3) are orthogonal trajectories of each other. In Figure 3.1 several members 
of the family of circles (3.2), drawn solidly, and several members of the family of 
straight lines (3.3), drawn with dashes, are shown. 

The problem of finding the orthogonal trajectories of a given family of curves arises 
in many physical situations. For example, in a two-dimensional electric field the lines of 
force (flux lines) and the equipotential curves are orthogonal trajectories of each other. 
We now proceed to find the orthogonal trajectories of a family of curves 

F(x,y,c) = 0. : (3.1) 

We obtain the differential equation of the family (3.1) by first differentiating Equation 
(3.1) implicitly with respect to x and then eliminating the parameter c between the 
derived equation so obtained and the given equation (3.1) itself. We assume that the 
resulting differential equation of the family (3.1) can be expressed in the form 

Thus the curve C of the given family (3.1) which passes through the point (x, y) has 
the slope f(x , y) there. Since an orthogonal trajectory of the given family intersects 
each curve of the family at right angles, the slope of the orthogonal trajectory to C at 

(*> y) is 

1 

H*,yY 




Figure 3.1 
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Thus the differential equation of the family of orthogonal trajectories is 

dy ^ 1 

dx ~ f(x , yY 

A one-parameter family 

G(x, y, c) = 0 



(3.5) 



or 



y = F(x,c) 

of solutions of the differential equation (3.5) represents the family of orthogonal 
trajectories of the original family (3.1), except possibly for certain trajectories that are 
vertical lines. 

We summarize this procedure as follows: 

Procedure for Finding the Orthogonal Trajectories of a 
Given Family of Curves 

Step I. F rom the equation 

F(x,y,c) = 0 (3.1) 

of the given family of curves, find the differential equation 

^ = /(*. y) (3.4) 



of this family. 



Step 2. In the differential equation dy/dx = /(x, y) so found in Step 1, replace 
f(x, y) by its negative reciprocal — l//(x, y). This gives the differential equation 



dy _ I 
dx ~ f{x, y) 

of the orthogonal trajectories. 



(3-5) 



Step J. Obtain a one-parameter family 

G(x, y y c) ~ 0 or y = F(x,c) 

of solutions of the differential equation (3.5), thus obtaining the desired family of 
orthogonal trajectories (except possibly for certain trajectories that are vertical lines 
and must be determined separately). 

Caution . In Step 1, in finding the differential equation (3.4) of the given family, be 
sure to eliminate the parameter c during the process. 



► Example 3.2 

In Example 3.1 we stated that the set of orthogonal trajectories of the family of circles 

x 2 + y 2 = c 2 (3.2) 
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is the family of straight lines 

y = kx. 

Let us verify this using the procedure outlined above. 
Step 1, Differentiating the equation 



of the given family, we obtain 



x 2 + y 2 = c 2 



dy 



From this we obtain the differential equation 



of the given family (3.2). (Note that the parameter c was automatically eliminated in this 
case.) 

Step 2. We replace — xjy by its negative reciprocal y/x in the differential equation 
(3.6) to obtain the differential equation 

T - 1 

dx x 

of the orthogonal trajectories. 

Step 3 . We now solve the differential equation (3.7). Separating variables, we have 

dy _ dx 

y * 1 

integrating, we obtain 

y - kx. (3.3) 

This is a one-parameter family of solutions of the differential equation (3.7) and thus 
represents the family of orthogonal trajectories of the given family of circles (3.2) 
(except for the single trajectory that is the vertical line x — 0 and this may be deter- 
mined by inspection). 



► Example 3.3 

Find the orthogonal trajectories of the family of parabolas y — cx 2 . 
Step I. We first find the differential equation of the given family 



y = cx 2 . 



Differentiating, we obtain 



& - 2 „. 

dx 



1 
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Eliminating the parameter c between Equations (3.8) and (3.9), we obtain the differ- 
ential equation of the family (3.8) in the form 



dy^2y 

dx x ' 



(3.10) 



Step 2. We now find the differential equation of the orthogonal trajectories by 
replacing 2 y/x in (3.10) by its negative reciprocal, obtaining 



dy _ x 
~dx = ~2y 



(3.11) 



Steps. We now solve the differential equation (3.11). Separating variables, we have 



2y dy = ~x dx. 

Integrating, we obtain the one-parameter family of solutions of (3.11) in the form 



x 2 + 2y 2 = k 2 , 

where k is an arbitrary constant. This is the family of orthogonal trajectories of (3.8); it 
is clearly a family of ellipses with centers at the origin and major axes along the x axis. 
Some members of the original family of parabolas and some of the orthogonal 
trajectories (the ellipses) are shown in Figure 3.2. 



B. Oblique Trajectories 

DEFINITION 

Let 



F(x , y,c) = 0 



<3.12) 
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be a one-parameter family of curves. A curve that intersects the curves of the family (3. 12) 
at a constant angle a # 90° is called an oblique trajectory of the given family. 

Suppose the differential equation of a family is 

| = (3.13) 

Then the curve of the family (3. 13) through the point (x, y)has slope /(x, y)at(x, y)and 
hence its tangent line has angle of inclination tan _1 [/(x, y)] there. The tangent line 
of an oblique trajectory that intersects this curve at the angle a will thus have angle 
of inclination 

tan~‘[/(x,y)] + a 

at the point (x, y). Hence the slope of this oblique trajectory is given by 

) / (x, y) 4- tan a 

M 

Thus the differential equation of such a family of oblique trajectories is given by 

dy /(x, y) 4- tan a 
dx 1 — /(x, y)tan a' 

Thus to obtain a family of oblique trajectories intersecting a given family of curves at 
the constant angle a ^ 90°, we may follow the three steps in the above procedure (page 
72) for finding the orthogonal trajectories, except that we replace Step 2 by the fol- 
lowing step: 

Step 2\ In the differential equation dy/dx = /(x, y) of the given family, replace 
/(x, y) by the expression 

/(x. y) + tan a 

1 — /(*> y)tan a’ ' ’ 



► Example 3.4 

Find a family of oblique trajectories that intersect the family of straight lines y ^ cx at 
angle 45°. 

Step I. From y = cx, we find dy/dx = c. Eliminating c, we obtain the differential 
equation 

T- 1 <3.<S 

dx x 

of the given family of straight lines. 

Step 2\ We replace /(x, y) =,y/x in Equation (3.15) by 

f{x, y) 4- tan a _ y/x + i x + y . 

\ - /(x, y)tan <x 1 - y/x ~ x - y 
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(tan a = tan 45° *= 1 here). Thus the differential equation of the desired oblique 
trajectories is 



dy _ x + y 
dx x — y' 



(3.16) 



Step 3. We now solve the differential equation (3.16). Observing that it is a 
homogeneous differential equation, we let y — vx to obtain 



dv 1 + v 
v + x— . 

dx 1 — 17 - 



After simplifications this becomes 



(v — 1) dv dx 

V 2 + 1 X ' 



Integrating we obtain 



^ ln(u 2 + 1) — arctan v = —In |x| —In |c| 



or 



In c 2 x 2 (t> 2 -hi) — 2 arctan v = 0. 

Replacing v by y/x, we obtain the family of oblique trajectories in the form 



In c 2 (x 2 + y 2 ) - 2 arctan - = 0. 



Exercises 

In Exercises 1-9 find the orthogonal trajectories of each given family of curves. In each 
case sketch several members of the family and several of the orthogonal trajectories on 
the same set of axes. 

I. y = cx 3 . 

3. cx 2 + y 2 = 1. 

5. y = x — 1 + ce~ x . 

7. x 2 +y 2 — cx 3 . 



10. Find the orthogonal trajectories of the family of ellipses having center at the 
origin, a focus at the point (c, 0), and semimajor axis of length 2c. 

1 1. Find the orthogonal trajectories of the family of circles which are tangent to they 
axis at the origin. 

12. Find the value of K such that the parabolas y — c A x 2 4- K are the orthogonal 
trajectories of the family of ellipses x 2 + 2y 2 — y = c 2 . 



2. y 2 = cx. 

4 y e cx 

6. x ~ y - cx 2 . 

8. x 2 =2y- 1 +c<T 2 '. 



\ 
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13. Find the value of n such that the curves x n + y n = c { are the orthogonal 
trajectories of the family 

x 

y = 

14. A given family of curves is said to be self-orthogonal if its family of orthogonal 
trajectories is the same as the given family. Show that the family of parabolas 
y 2 — 2 cx + c z is self-orthogonal. 

15. Find a family of oblique trajectories that intersect the family of circles x 2 -\-y 2 — 
c 2 at angle 45°. 

16. Find a family of oblique trajectories that intersect the family of parabolas y 2 — cx 
at angle 60°. 

17. Find a family of oblique trajectories that intersect the family of curves x+y — 
cx 2 at’ angle a such that tan a = 2. 



3.2 PROBLEMS IN MECHANICS 
A. Introduction 

Before we apply our knowledge of differential equations to certain problems in 
mechanics, let us briefly recall certain principles of that subject. The momentum of a 
body is defined to be the product mv of its mass m and its velocity v. The velocity v and 
hence the momentum are vector quantities. We now state the following basic law of 
mechanics: 




F = kma> 



(3.18) 
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where k — \/K and a = dv/dt is the acceleration of the body. The form (3.17) is a direct 
mathematical statement of the manner in which Newton's second law is usually 
expressed in words, the mass being considered constant. However, we shall make use of 
the equivalent form (3. 1 8). The magnitude of the constant of proportionality k depends 
upon the units employed for force, mass, and acceleration. Obviously the simplest 
systems of units are those for which k — L When such a system is used (3T8) reduces to 

F — ma. (3.19) 

It is in this form that we shall use Newton's second law. Observe that Equation (3. 1 9) is 
a vector equation. 

Several systems of units for which k = 1 are in use. In this text we shall use only three: 
the British gravitational system (British), the centimeter-gram-second system (cgs), and 
the meter-kilogram-second system (mks). We summarize the various units of these 
three systems in Table 3.1. v 

Recall that the force of gravitational attraction that the earth exerts on a body is 
called the weight of the body. The weight, being a force, is expressed in force units. Thus 
in the British system the weight is measured in pounds; in the cgs system, in dynes; and 
in the mks system, in newtons. 

Let us now apply Newtons’s second law to a freely falling body (a body falling 
toward the earth in the absence of air resistance). Let the mass of the body be m and let 
w denote its weight. The only force acting on the body is its weight and so this is the 
resultant force. The acceleration is that due to gravity, denoted by g, which is 
approximately 32 ft/sec 2 in the British system, 980 cm/sec 2 in the cgs system, and 
9.8 m/sec 2 in the mks system (for points near the earth’s surface). Newton’s second law 
F — ma thus reduces to w = mg. Thus 

m - (3.20) 

, 9 

a relation that we shall frequently employ. 

Let us now consider a body B in rectilinear motion, that is, in motion along a straight 
line L. On L we choose a fixed reference point as origin 0, a fixed direction as posi- 
tive, and a unit of distance. Then the coordinate x of the position of B from the origin 0 
tells us the distance or displacement of B. (See Figure 3.3.) The instantaneous velocity of 
B is the time rate of change of x: 



TABL^ 3.1 





British System 


cgs System 


mks System 


force 


pound 


dyne 


newton 


mass 


slug 


gram 


kilogram 


distance 


foot 


centimeter 


meter 


time 


second ' 


second 


second 


acceleration 


ft/sec 2 


cm /sec 2 


m/sec 2 
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Figure 3.3 



3 



and thejnstanfaneous acceleration of 5 is the time rate of change of v: 

dv _ d 2 x 
a dt dt 1 



Note that x,v , and a are vector quantities. All forces, displacements, velocities, and 
accelerations in the positive direction on L are positive quantities; while those in the 
negative direction are negative quantities. 

If we now apply Newton’s second law F = ma to the motion of B along L, noting 
that 



dv dv dx dv 

dt dx dt dx ’ 



we may express the law in any of the following three forms: 

dv 
m — 
dt 

d 2 x 

dv 
mv — 
dx 



= F, 


(3.21) 


= F, 


(3-22) 


= F, 


(3-23) 



where F is the resultant force acting on the body. The form to use depends upon the way 
in which F is expressed. For example, if F is a function of time t only and w r e desire to 
obtain the velocity t; as a function off, we would use (3.21); whereas if F is expressed as a 
function of the displacement x and we wish to find v as a function of x, we would 
employ (3,23). 



8. Falling Body Problems 

We shall now consider some examples of a body falling through air toward the earth. In 
such a circumstance the body encounters air resistance as it falls. The amount of air 
resistance depends upon the velocity of the body, but no general law exactly expressing 
this dependence is known. In some instances the law R - kv appears to be quite 
satisfactory, while in others R = kv 2 appears to be more exact. In any case, the con- 
stant of proportionality k in turn depends on several circumstances. In the ex- 
amples that follow we shall assume certain reasonable resistance laws in each case. 
Thus we shall actually be dealing with idealized problems in which the true resis- 
tance law is approximated and in which certain comparatively negligible factors are 
disregarded. 



1 
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► Example 3.5 

A body weighing 8 ib falls from rest toward the earth from a great height. As it falls, air 
resistance act upon it, and we shall assume that this resisience (in pounds) is 
numerically equal to 2v, where v is the velocity (in feet per second). Find the velocity and 
distance fallen at time t seconds. 



Formulation. We choose the positive x axis vertically downward along the path of 
the body B and the origin at the point from which the body fell. The forces acting on the 
body are: 

1. F u its weight, 8 lb, which acts downward and hence is positive. 

2. F 2 , the air resistance, numerically equal to 2v , which acts upward and hence is the 
negative quantity — 2tx 

See Figure 3.4, where these forces are indicated. 



Newton's second law, F — ma , becomes 



dv 

or, taking g = 32 and using m — w/g = 

1 dv 

- — = 8 - 2u. 

4 dt 

Since the body was initially at rest, we have the initial condition 

v(0) = 0. 



(124) 



(3.25) 



0 



x 



F 2 - —2u 
Fi = 8 



+ 



Earth 

Figure 3.4 
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Solution. Equation (3.24) is separable. Separating variables, we have 



dv . 
8~— 2u 



= 4 dv 



Integrating we find 

— ^ In 1 8 — 2v\ = 4t + c 0 , 

which reduces to 



8 -2v = c 1 e‘ st . 

Applying the condition (3.25) we find c y — 8. Thus the velocity at time t is given by 

u - 4(1 - e“ 8 '). (3.26) 

Now to determine the distance fallen at time £, we write (3.26) in the form 

and note that x(0) - 0. Integrating the above equation, we obtain 

x = 4(t + £e“ 8 *) + c 2 . 

Since x = 0 when t — 0, we find c 2 = — j and hence the distance fallen is given by 

X = 4 (t + ir 8 ' - £). (3.27) 

Interpretation of Results. Equation (3.26) shows us that as t -» oo, the velocity v 
approaches the limiting velocity 4(ft/sec). We also observe that this limiting velocity is 
approximately attained in a very short time. Equation (3:27) states that a st-+ oo,x also 
oo. Does this imply that the body will plow through the earth and continue forever? 
Of course not; for when the body reaches the earth’s surface its motion will certainly 
cease. How then do we reconcile this obvious end to the motion with the statement of 
Equation (3.27)? It is simple; when the body reaches the earth’s surface, the differential 
equation (3.24) and hence Equation (3.27) no longer apply? 

► Example 3.6 

A skydiver equipped with parachute and other essential equipment falls from rest 
toward the earth. The total weight of the man plus the equipment is 160 lb. Before the 
parachute opens, the air resistance (in pounds) is numerically equal to where v is the 
velocity (in feet per second). The parachute opens 5 sec after the fall begins; after it 
opens, the air resistance (in pounds) is numerically equal to ju 2 , where v is the velocity 
(in feet per second). Find the velocity of the skydiver (A) before the parachute opens, 
and (B) after the parachute opens. 

Formulation. We again choose the positive x axis vertically downward with the 
origin at the point where the fall began. The statement of the problem suggests that we 
break it into two parts: (A) before the parachute opens; (B) after it opens. 

We first consider problem (A). Before the parachute opens, the forces acting upon the 
skydiver are: 

1. F u the weight, 160 lb, which acts downward and hence is positive. 
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2. F 2j the air resistance, numerically equal to which acts upward and hence is the 
negative quantity — iu. 

We use Newton’s second law F ~ ma , where F = F t + F 2> let m = w/g, and take 
g = 32. We obtain 

dv 

5 - 7 - = 160 -U 
dt 2 

Since the skydiver was initially at rest, v = 0 when t — 0. Thus, problem (A), concerned 
with the time before the parachute opens, is formulated as follows: 

5~=160-ir. (3.28) 

u(0) = 0. (3.29) 

We now turn to the formulation of problem (B). Reasoning as before, we see that 
after the parachute opens, the forces acting upon the skydiver are: 

L F^- 160, exactly as before. 

2, F z = -ft; 2 (instead of — \v). 

Thus, proceeding as above, we obtain the differential equation 

S T , *'■ 

Since the parachute opens 5 sec after the fall begins, we have v = v l when t = 5, where 
is the velocity attained when the parachute opened. Thus, problem (B), concerned 
with the time after the parachute opens, is formulated as follows: 

5 ~ = 160 — $t> 2 , (3.30) 

t>(5) = (3.31) 

Solution. We shall first consider problem (A). We find a one-parameter family of 
solution of 

5^=160-^. (3.28) 

Separating variables, we obtain 



a — 320 



= -To dt. 



Integration yields 

ln(« - 320) = — j'gt + c 0 , 
which readily simplifies to the form 

v = 320 + ce -,,, °. 
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Applying the initial condition (3.29) that v = 0 at t = 0, we find that c — —320. Hence 
the solution to problem (A) is 

t; - 320(1 - <T ,/10 ), (3.32) 

which is valid for 0 ^ t < 5. In particular, where t = 5, we obtain 

»,= 320(1 -e~ 1/2 )» 126, (3.33) 

which is the velocity when the parachute opens. 

Now let us consider problem (B). We fiist find a one-parameter family of solutions of 
the differential equation 

5~r= 160 — (3.30) 

at 



Simplifying and separating variables, we obtain 



Integration yields 



dv _ dt 
u 2 — 256 “ 



1 v — 16 
32 n d + 16 — ~8 + 



or 



In 



v — 16 
v+ 16 



-4t 



This readily simplifies to the form 



v — 16 
v ( + 16 



= ce 



- 4r 



and solving this for v we obtain 



16(ce -4 ' + 1) 
1 — ce~ M 



(3.34) 



(3.35) 



Applying the initial condition (3.31) that v = v, at t = 5, where i>, is given by (3.33) and 
is approximately 126, to (3.34), we obtain 

c = Me 20 . 



Substituting this into (3.35) wc obtain 



16(H§e 20 ~ 4t + 1) 

l-ttS * 20 ' 41 ’ 



(3.36) 



which is valid for t £ 5. 



Interpretation of Results. Let us first consider the solution of problem (A), given 
by Equation (3.32). According to this, as / -+ oo, v approaches the limiting velocity 
320 ft/sec. Thus if the parachute never opened, the velocity would have been approx- 
imately 320 ft/sec at the time when the unfortunate skydiver would have struck the 
earth! But, according to the statement of the problem, the parachute does open 5 sec 



84 APPLICATIONS OF FIRST-ORDER EQUATIONS 



after the fall begins (we tacitly and thoughfully assume 5 « T, where T is the time 
when the earth is reached!). Then, referring to the solution of problem (B), Equation 
(3.36), we see that as t -► ao, v approaches the limiting velocity 1 6 ft/sec. Thus, assuming 
that the parachute opens at a considerable distance above the earth, the velocity is 
approximately 16 ft/sec when the earth is finally reached. We thus obtain the' well- 
known fact that the velocity of impact with the open parachute is a small fraction of 
the impact velocity that would have occured if the parachute had not opened. The 
calculations in this problem are somewhat complicated, but the moral is clear: Make 
certain that the parachute opens! 



C. Frictional Forces 

If a body moves on a rough surface, it will encounter not only air resistance but also 
another resistance force due to the roughness of the surface. This additional force is 
called friction . It is shown in physics that the friction is given by /iN, where 

1. fj. is a constant of proportionality called the coefficient of friction , which depends 
upon the roughness of the given surface; and 

2. N is the normal (that is, perpendicular) force which the surface exerts on the body. 

We now apply Newton’s second law to a problem in which friction is involved. 



► Example 3.7 

An object weighing 48 lb is released from rest at the top of a plane metal slide that is 
inclined 30° to the horizontal. Air resistance (in pounds) is numerically equal to one- 
half the velocity (in feet per second), and the coefficient of friction is one-quarter. 



A. What is the velocity of the object 2 sec after it is released? 

B. If the slide is 24 ft long, what is the velocity when the object reaches the bottom? 

Formulation. The line of motion is along the slide. We choose the origin at the top 
and the positive x direction down the slide. If we temporarily neglect the friction and air 
resistance, the forces acting upon the object A are: 



1. Its weight, 48 lb, which acts vertically downward; and 

2. The normal force, N, exerted by the slide which acts in an upward direction 
perpendicular to the slide. (See Figure 3.5.) 



The components of the weight parallel and perpendicular to the slide have magnitude 

48 sin 30° = 24 



and 

48 cos 30° = 24^/3, 

respectively, The components perpendicular to the slide are in equilibrium and hence 
the normal force N has magnitude 24 x /3. 
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. Now, taking into consideration the friction and air resistance, we see that the forces 
acting on the object as it moves along the slide are the following: 

1. F l3 the component on the weight parallel to the plane, having numerical value 24. 
Since this force acts in the positive (downward) direction along the slide, we have 

F, - 24. 

2. F 2 , the frictional force, having numerical value juN = £(24^/3). Since this acts in the 
negative (upward) direction along the side, we have 

f 2 = -6^/3. 

3. F 3 , the air resistance, having numerical value Since v > 0 and this also acts in 
the negative direction, we have 

We apply Newton’s second law F — ma. Here F — F { + F 2 + F 3 = 24 — 6^/3 — 
and m = w/g = ff = |. Thus we have the differential equation 

~ - 24 - 6^3 - *«>. (3.37) 

Since the object is released from rest, the initial condition is 

u(0) = 0. (3.38) 

Solution. Equation (3.37) is separable; separating variables we have 

dv _ dt 

48 - 12 v /3 - v ~ 3 ' 

Integrating and simplifying, we find 

p = 48-12 v /3-c 1 e'" 3 . 

The condition (3.38) gives c t = 48 - 12^/3. Thus we obtain 

v = (48 - 1273X1 - e-" 3 ). (3.39) 

Question A is thus answered by letting r = 2 in Equation (3.39). We find 
c(2) = (48 — 12 v /i)(l - e~ 213 ) * !0.2(ft/sec). 
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In order to answer question B, we integrate (3.39) to obtain 

x - (48 - iiy/y&t + 3*r' /3 ) + c 2 . 

Since x(0) = 0, c 2 = —(48 — 12 x /3)(3). Thus the distance covered at time t is given by 
x = (48 - + 3c'" 3 - 3). 

Since the slide is 24 ft long, the object reaches the bottom at the time T determined from 
the transcendental equation 

24 = (48 - 12 n / 3)(7’ + 3e" T/3 - 3), 

which may be written as 

2e~ Tf3 = 47 + _ j 

13 

The value of T that satisfies this equation is approximately 2.6. Thus from Equation 
(3.39) the velocity of the object when it reaches the bottom is given approximately by 

(48 - 12 x /3)(l - e~ 0 ' 9 ) * 12.3 (ft/sec). 



Exercises 

1. A stone weighing 4 lb falls from rest toward the earth from a great height. As it 
falls it is acted upon by air resistance that is numercially equal to jv (in pounds), 
where v is the velocity (in feet per second). 

(a) Find the velocity and distance fallen at time t sec. 

(b) Find the velocity and distance fallen at the end of 5 sec. 

2. A ball weighing 6 lb is thrown vertically downward toward the earth from a height 
of 1000 ft with an initial velocity of 6 ft/sec. As it falls it is acted upon by air 
resistance that is numerically equal to f v (in pounds), wherey is the velocity (in feet 
per second). 

(a) What is the velocity and distance fallen at the end of one minute? 

(b) With what velocity does the ball strike the earth? 

3. A ball weighing | lb is thrown vertically upward from a point 6 ft above the 
surface of the earth V/ith an initial velocity of 20 ft/sec. As it rises it is acted upon 
by air resistance that is numerically equal to -^v (in pounds), where v is the 
velocity (in feet per second). How high will the ball rise? 

4. A ship which weighs 32,000 tons starts from rest under the force of a constant 
propeller thrust of 100,000 lb. The resistance in pounds is numerically equal to 
8000u, where v is in feet per second. 

(a) Find the velocity of the ship as a function of the time. 

(b) Find the limiting velocity (that is, the limit of v as t -* + oo). 

(c) Find how long it takes the ship to attain a velocity of 80% of the limiting 
velocity. 

5. A body of mass 100 g is dropped from rest toward the earth from a height of 
1000 m. As it falls, air resistance acts upon it, and this resistance (in newtons) is 
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proportional to the velocity v (in meters per second). Suppose the limiting velocity 
is 245 m/sec. 

(a) Find the velocity and distance fallen at time t secs. 

(b) Find the time at which the velocity is one-fifth of the limiting velocity. 

6. An object of mass 100 g is thrown vertically upward from a point 60 cm above the 
earth’s surface with an initial velocity of 150 cm/sec. It rises briefly and then falls 
vertically to the earth, all of which time it is acted on by air resistance that is 
numerically equal to 200i> (in dynes), where v is the velocity (in cm/sec). 

(a) Find the velocity 0.1 sec after the object is thrown. 

(b) Find the velocity 0.1 sec after the object stops rising and starts falling. 

7. Two people are riding in a motorboat and the combined weight of individuals, 
motor, boat, and equipment is 640 lb. The motor exerts a constant force of 20 lb 
on the boat in the direction of motion, while the resistance (in pounds) is 
numerically equal to one and one-half times the velocity (in feet per second). If the 
boat started from rest, find the velocity of the boat after (a) 20 sec, (b) l min. 

8. A boat weighing 150 lb with a single rider weighing 170 lb is being towed in a 
certain direction at the rate of 20 mph. At time t = 0 the tow rope is suddenly cast 
off and the rider begins to row in the same direction, exerting a force equivalent to 
a constant force of 1 2 lb in this direction. The resistance (in pounds) is numerically 
equal to twice the velocity (in feet per second). 

(a) Find the velocity of the boat 1 5 sec after the tow rope was cast off. 

(b) How many seconds after the tow rope is cast off will the velocity be one-half 

that at which the boat was being towed? 

9. A bullet weighing 1 oz is fired vertically downward from a stationary helicopter 
with a muzzle velocity of 1200 ft/sec. The air resistance (in pounds) is numerically 
equal to 16 " 5 v 2 , where v is the velocity (in feet per second). Find the velocity of the 
bullet as a function of the time. 

10. A shell weighing 1 lb is fired vertically upward from the earth’s surface with a 
muzzle velocity of 1000 ft/sec. The air resistance (in pounds) is numerically equal 
to 10” 4 y 2 , where v is the velocity (in feet per second). 

(a) Find the velocity of the rising shell as a function of the time. 

(b) How long will the shell rise? 

11. An object weighing 16 lb is dropped from rest on the surface of a calm lake and 
thereafter starts to sink. While its weight tends to force it downward, the 
buoyancy of the object tends to force it back upward. If this buoyancy force is one 
of 6 lb and the resistance of the water (in pounds) is numerically equal to twice the 
square of the velocity (in feet per second), find the formula for the velocity of the 
sinking object as a function of the time. 

12. An object weighing 12 lb is placed beneath the surface of a calm lake. The 
buoyancy of the object is 30 lb; because of this the object begins to rise. If 
the resistance of the water (in pounds) is numerically equal to the square of the 
velocity (in feet per second) and the object surfaces in 5 sec, find the velocity of the 
object at the instant when it reaches the surface. 
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13. A man is pushing a loaded sled across a level field of ice at the constant speed of 
10 ft/sec. When the man is halfway across the ice field, he stops pushing and lets 
the loaded sled continue on. The combined weight of the sled and its load is 80 lb; 
the air resistance (in pounds) is numerically equal to where v is the velocity of 
the sled (in feet per second); and the coefficient of friction of the runners on the ice 
is 0.04. How far will the sled continue to move after the man stops pushing? 

14. A girl on Jier sled has just slid down a hill onto a level field of ice and is starting to 
slow down. At the instant when their speed is 5 ft/sec, the girl’s father runs up and 
begins to push the sled forward, exerting a constant force of 15 lb in the direction 
of motion. The combined weight of the girl and the sled is 96 lb, the air resistance 
(in pounds) is numerically equal to one-half the velocity (in feet per second), and 
the coefficient of friction of the runners on the ice is 0.05. How fast is the sled 
moving 10 sec after the father begins pushing? 

15. A case of canned milk weighing 24 lb is released from rest at the top of a plane 
metal slide which is 30 ft long and inclined 45° to the horizontal. Air resistance (in 
pounds) is numerically equal to one-third the velocity (in feet per second) and the 
coefficent of friction is 0.4, 

(a) What is the velocity of the moving case 1 sec after it is released? 

(b) What is the velocity when the case reaches the bottom of the slide? 

16. A boy goes sledding down a long 30° slope. The combined weight of the boy and 
his sled is 72 lb and the air resistance (in pounds) is numerically equal to twice their 
velocity (in feet per second). If they started from rest and their velocity at the end 
of 5 sec is 10 ft/sec, what is the coefficient of friction of the sled runners on the 
snow? 

17. An object weighing 32 lb is released from rest 50 ft above the surface of a calm 
lake. Before the object reaches the surface of the lake, the air resistance (in pounds) 
is given by 2v , where v is the velocity (in feet per second). After the object passes 
beneath the surface, the water resistance (in pounds) is given by 6v. Further, the 
object is then buoyed up by a buoyancy force of 8 lb. Find the velocity of the 
object 2 sec after it passes beneath the surface of the lake. 

18. A rocket of mass m is fired vertically upward from the surface of the earth with 
initial velocity v — v 0 . The only force on the rocket that we consider is the 
gravitational attraction of the earth. Then, according to Newton’s law of 
gravitation, the acceleration a of the rocket is given by a = —k/x 2 , where k > 0 
is a constant of proportionality and x is the distance “upward” from the center 
of the earth along the line of motion. At time t = 0, x - R (where R is the radius 
of the earth), a - —g (where g is the acceleration due to gravity), and v - v 0 . 
Express a = dv/dt as in Equation (3.23), apply the appropriate initial data, and 
note that v satisfies the differential equation 

dv _ gR 2 
dx x 2 

Solve this differential equation, apply the appropriate initial condition, and thus 
express v as a function of x. In particular, show that the minimum value of r 0 for 
which the rocket will escape from the earth \s^^2gR. This is the so-called veloc - 
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ity of escape ; and using JR = 4000 miles, g - 32 ft /sec 2 , one finds that this is 
approximately 25,OOOmiph (or 7 mi/sec). 

19. A body of mass m is m rectilinear motion along a horizontal axis. The resultant 
force acting on the body is given by — kx , where k > 0 is a constant of 
proportionality and x is the distance along the axis from a fixed point O. The body 
has initial velocity v = v 0 when x ~ x 0 . Apply Newton’s second law in the form 
(3.23) and thus write the differential equation of motion in the form 

dv 

mv — — —kx. 
dx 

Solve the differential equation, apply the initial condition, and thus express the 
square of the velocity v as a function of the distance x. Recalling that v = dx/dt , 
show that the relatioirbetween v and x thus obtained is satisfied for all time t by 

where </> is a constant. 



3.3 RATE PROBLEMS 

In certain problems the rate at which a quantity changes is a known function of the 
amount present and/or the time, and it is desired to find the quantity itself. If x denotes 
the amount of the quantity present at time t, then dx/dt denotes the rate at which the 
quantity changes and we are at once led to a differential equation. In this section we 
consider certain problems of this type. 

A. Rate of Growth and Decay 

► Example 3.8 

The rate at which radioactive nuclei decay is proportional to the number of such nuclei 
.that are present in a given sample. Half of the original number of radioactive nuclei 
have undergone disintegration in a period of 1500 years, 

1 . What percentage of the original radioactive nuclei will remain after 4500 years? 

2. In how many years will only one-tenth of the original number remain? 

Mathematical Formulation. Let x be the amount of radioactive nuclei present 
after t years. Then dx/dt represents the rate at which the nuclei decay. Since the nuclei 
decay at a rate proportional to the amount present, we have 

~~ — Kx, (3.40) 

f/\ dt 

where K is a constant of proportional ity. The amount x is clearly positive; further, since 
x is decreasing, dx/dt < 0. Thus, from Equation (3.40), we must have K < 0. In order to 
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emphasize that x is decreasing, vve prefer to replace K by a positive constant preceded 
by a minus sign. Thus we let fc == — K > 0 and write the differential equation (3.40) in 
the form 

^ =-kx. (3.41) 

dt 

Letting x 0 denote the amount initially present, we also have the initial condition 

x(0) = x o . (3.42) 

We know that we shall need such a condition in order to determine the arbitrary 
constant that will appear in ai one-parameter family of solutions of the differential 
equation (3.41). However, we shall apparently need something else, for Equation (3.41) 
contains an unknown constant of proportionality k. This “something else” appears in 
the statement of the problem, for we are told that half of the original number 
disintegrate in 1 500 years. Thus half also remain at that time, and this at once gives the 
condition 

x(1500) = jx 0 ‘ (3.43) 

Solution. The differential equation (3.41) is clearly separable; separating variables, 
integrating, and simplifying, we have at once 



Applying the initial condition (3.42), x - x 0 when t = 0, we find that c - x 0 and hence 
we obtain 



We have not yet determined k . Thus we now apply condition (3.43), x = jx 0 when 
t = 1500, to Equation (3.44). We Ifind 

± Y _ x « “ 1 500k 

2' v 0 *“ i 



~fc>l500 _ 1 

? J — 2» 



or finally 



From this equation we could deitermine k explicitly and substitute the result into 
Equation (3.44). However, we see f rom Equation (3.44) that we actually do not need k 
itself but rather only e ~ k , which vve have just obtained in Equation (3.45). Thus we 
substitute e~ k fr om (3.45) into (3.44) to obtain 

x = x 0 (e- k y = x 0 [$y'“°°y 



x *= xS)' n ™ (3.46) 

Equation (3.46) gives tlhe number x of radioactive nuclei that are present at time r. 
Question 1 asks us what percentage of the original number will remain after 4500 years. 
We thus let t — 4500 in Equation (3.46> and find 

x — -*o(l ) 3 ~ 8 X 0* 




t 



1500 In 10 
In 2 



« 4985 (years). 



B. Population Growth 

We next consider the growth of a population (for example, human, an animal species, 
or a bacteria colony) as a function of time. Note that a population actually increases 
discontinuously by whole number amounts. However, if the population is very large, 
such individual increases in it are essentially negligible compared to the entire pop- 
ulation itself. In other words, the population increase is approximately continuous. 
We shall therefore assume that this increase is indeed continuous and in fact that the 
population is a continuous and differentiable fucntion of time. 

Given a population, we let x be the number of individuals in it at time t. If we assume 
that the rate of change of the population is proportional to the number of individuals 
in it at any time, we are led to the differential equation 

~ = kx, (3.47) 

dt 

where k is a constant of proportionality. The population x is positive and is increasing 
and hence dx/dx > 0. Therefore, from (3.47), we must have k > 0. Now suppose that at 
time f 0 the population is x 0 . Then, in addition to the differential equation (3.47), we 
have the initial condition 

*('o) = *o- (3.48) 

The differential equation (3.47) is separable. Separating variables, integrating, and 
simplifying, we obtain 

x = ce k1 . 

Applying the initial condition (3.48), x = x 0 at f = f 0 , to this, we have x 0 = ce kt0 . From 
this we at once find c = x 0 e' kt ° and hence obtain the unique solution 

x - x 0 e ku ~ to) (3.49) 

of the differential equation (3.471, which satisfies the initial condition (3.48). 
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From (3.49) we see that a population governed by the differential equation (3.47) with 
k > 0 and initial condition (3.48) is one that increases exponentially with time. This law 
of population growth is called the Malthusian law . We should now inquire whether or 
not there are cases in which such a model for population growth is indeed realistic. In 
answer to this, it can be shown that this model, with a suitable value of k , is remarkably 
accurate in the case of the human population of the earth during the last several 
decades (see Problem 8(b)). U is also known to be outstandingly accurate for certain 
mammalian species, with suitable fc, under certain realizable conditions and for certain 
time periods. On the other hand, turning back to the case of the human population of 
the earth, it can be shown that the Malthusian law turns out to be quite unreasonable 
when applied to the distant future (see Problem 8(e)). It is also completely unrealistic for 
other populations (for example, bacteria colonies) when applied over sufficiently long 
periods of time. The reason for this is not hard to see. For, according to (3.49), a 
population modeled by this law always increases and indeed does so at an ever 
increasing rate; whereas observation shows that a given population simply does not 
grow indefinitely. 

Population growth is represented more realistically in many cases by assuming that 
the number of individuals x in the population at time t is described by a differential 
equation of the form 

d.x 

— = kx — Ax 2 , (3.50) 

at 

where k > 0 and A > 0 are constants. The additional term — Ax 2 is the result of some 
cause that tends to limit the ultimate growth of the population. For example, such a 
cause could be insufficient living space or food supply, when the population becomes 
sufficiently large. Concerning the choice of — Ax 2 for the term representing the effect of 
the cause, one can argue as follows : Assuming the, cause affects the entire population of 
x members, then the effect on any one individual is proportional to x. Thus the effect on 
all x individuals in the population would be proportional to x ■ x = x 2 . 

We thus assume that a population is described by a differential equation of the form 
(3.50), with constants k > 0 and A > 0, and an initial condition of the form (3.48). In 
most such cases, it turns out that the constant A is very small compared to the constant 
k. Thus for sufficiently small x, the term kx predominates, and so the population grows 
very rapidly for a time. However, when x becomes sufficiently large, the term - Ax 2 is of 
comparatively greater influence, and the result of this is a decrease in the rapid growth 
rate. We note that the differential equation (3.50) is both a separable equation and a 
Bernoulli equation. The law of population growth so described is called the logistic law 
of growth. We now consider a specific example of this type of growth. 



► Example 3 »9 

The population x of a certain.city satisfies the logistic law 

dx _ 1 1 2 

Tt~ Too x_ (Toj rX 



(3.51) 



where time t is measured in years. Given that the population of this city is 100,000 in 
1980, determine the population as a function of time for t > 1980. In particular, answer 
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the following questions: 

(a) What will be the population in 2000? 

(b) In what year does the 1980 population double? 

(c) Assuming the differential equation (3.51) applies for all t > 1980, how large will 
the population ultimately be? 



Solution. We must solve the separable differential equation (3.51) subject to the 
initial solution 



x(1980) = 100,000. 

Separating variables in (3.51), we obtain 

dx 



and hence 



(10)~ 2 x — (10)~ 8 x 2 
dx 






dt. 



(10)- 2 x[l -(10 T 6 x] 

Using partial fractions, this becomes 

,Jl + / L.* 

1_X 1 - (10) 6 xJ 

Integrating, assuming 0 < x < 10 6 , we obtain 

100{ln x - ln[l - (10) 6 x]} ~t + Cl 



and hence 



Thus we find 






1 - (10)' 6 x ce 

Solving this for x, we finally obtain 



_ retnoo 



x — ' 



ce 



,*/ioo 



1 + (10)^ 6 ce J/i00, 
Now applying the initial condition (3.52) to this, we have 



(10) 5 



ce 



from which we obtain 



!■ 



c ~ 



1 +(I0)- 6 ce 19 *’ 

(10) 5 (10) 6 



i ,9 8 [l -(10) 5 (10)- 6 ] 9e 



19. e • 



(3-52) 



(3.53) 



Substituting this value for c back into (3.53) and simplifying, we obtain the solution in 
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the form 



x =5 



(10) 6 



1 + 9 e l9.8-r/l00 • 



(3.54) 



This gives the population x as a function of time for t > 1980. 

We now consider the questions (a), (b), and (c) of the problem. Question (a) asks for 
the population in the year 2000. Thus we let t = 2000 in (3.54) and obtain 



x — 



(10) 6 



1 + 9<T°‘ 2 



119,495. 



Question (b) asks for the year in which the population doubles. Thus we let 
x = 200,000 = 2(10) 5 in (3.54) and solve for t. We have 



2 ( 10) 5 » 



m 6 



from which 



and hence 



1 + 9^19.8-1/100’ 

,,19.8-r/iOO __ 4 
t 9 ) 

f « 2061. 



Question (c) asks how large the population will ultimately be, assuming the differential 
equation (3.51) applies for all t > 1980. To answer this, we evaluate limx as t -> co using 
the solution (3.54) of (3.51). We find 



.. r ( 10) 6 



= ( 10) 6 = 1 , 000 , 000 . 



C. Mixture Problems 

We now consider rate problems involving mixtures. A substance S is allowed to flow 
into a certain mixture in a container at a certain rate, and the mixture is kept uniform by 
stirring. Further, in one such situation, this uniform mixture simultaneously flows out 
of the container at another (generally different) rate; in another situation this may not 
be the case. In either case we seek to determine the quantity of the substance S present 
in the mixture at time t. 

Letting x denote the amount of S present at time t , the derivative dx/dt denotes the 
rate of change of x with respect to t. If IN denotes the rate at which S enters the mixture 
and OUT the rate at which it leaves, we have at once the basic equation 

^ = IN — OUT (3.55) 

dt 

from which to determine the amount x of S at time t. We now consider examples. 

► Example 3.10 

A tank initially contains 50 gal of pure water. Starting at time t = 0 a brine containing 
2 lb of dissolved salt per gallon flows into the tank at the rate of 3 gal/min. The mixture 
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is kept uniform by stirring and the well-stirred mixture simultaneously flows out of the 
tank at the same rate. 

1. How much salt is in the tank at any time t > 0? 

2. How much salt is present at the end of 25 min? 

3. How much salt is present after a long time? 

Mathematical Formulation. Let x denote the amount of salt in the tank at time t. 
We apply the basic equation (3.55), 

^ = IN — OUT. 

dt 

The brine flows in at the rate of 3 gal/min, and each gallon contains 2 lb of salt. Thus 
IN = (2 lb/gal)(3 gal/min) = 6 Ib/min. 

Since the rate of outflow equals the rate of inflow, the tank contains 50 gal of the 
mixture at any time r. This 50 gal contains x lb of salt at time f, and so the concentration 
of salt at time t is jqX lb/gal. Thus, since the mixture flows out at the rate of 3 gal/min, we 
have 



- (55 ib/gai ) 



(3 gal/min) = ^ lb/min. 



Thus the differential equation for x as a function of t is 



Since initially there was no salt in the tank, we also have the initial condition 

x(0) = 0. (3.57) 

Solution. Equation (3.56) is both linear and separable. Separating variables, we 
have 

dx _ 3 
100 -x = 50 dL 

Integrating and simplifying, we obtain 

x ~ 100 + ce‘ 3f/5 °. 

Applying the condition (3.57), x = 0 at t — 0, we find that c = — 100. Thus we have 

x = 100(1 -<r 3 " 50 ). (3.58) 

This is the answer to question 1. As for question 2, at the end of 25 min, t — 25, and 
Equation (3.58) gives 

x(25)= 100(1 5 )* 78(ib). 

Question 3 essentially asks us how much salt is present as t -> 00 . To answer this we let 
t -► 00 in Equation (3.58) and observe that x -* 100. 
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► Example 3.1 1 

A large tank initially contains 50 gal of brine in which there is dissolved 10 lb of salt. 
Brine containing 2 lb of dissolved salt per gallon flows into the tank at the rate of 
5 gal/min. The mixture is kept uniform by stirring, and the stirred mixture simulta- 
neously flows out at the slower rate of 3 gal/min. How much salt is in the tank at any 
time t > 0? 

Mathematical Formulation. Let x = the amount of salt at time t . Again we shall 
use Equation (3.55): 

^ = IN — OUT. 

dt 

Proceeding as in Example 3.10, 



IN = (2 lb/gal)(5 gal/min) = 10 Ib/min; 



Iso, once again 



OUT = (C lb/gal)(3 gal/min), 

where C lb/gal denotes the concentration. But here, since the rate of outflow is different 
from that of inflow, the concentration is not quite so simple. At time t = 0, the tank 
contains 50 gal of brine. Since brine flows in at the rate of 5 gal/min but flows out at the 
slower rate of 3 gal/min, there is a net gain of 5 — 3 = 2 gal/min of brine in the tank. 
Thus at the end of t minutes the amount of brine in the tank is 

50 + 2 1 gal. 

Hence the concentration at time t minutes is 



50 + It 



lb/gal. 



and so 



OUT - — — Ib/min. 

50 + It 

Thus the differential equation becomes 

^=10 *L_ 

dt 50 + 2f 

Since there was initially 10 lb of salt in the tank, we have the initial condition 

x(0) = 10. 



Solution. The differential equation (3.59) is not separable but it is linear. Putting it 
in standard form, 

dx 3 

— + — x - 10, 

dt 2 1 50 
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we find the integrating factor 



^}dho‘' , )‘ (2,+SO) ’” 

Multiplying through by this, we have 

(2f + 50) 3/2 ~ + 3(2 1 + 50) 1/2 x = I0(2f + 50) 3/2 
at 

or 



Thus 



^ [(2 1 + 50) 3/2 x] = 10(2r + 50) 3/2 . 
dt 

(2 1 + 50) 3/2 a: = 2(2f + 50) 5/2 + c 



or 



x = 4(f + 25) + • 



(2t + 50) 3/2 ‘ 

Applying condition (3.60), x = 10 at r — 0, we find 

c 



10 ^ 100 + 



(50) 3/2 



or 



c = — (90)(50) 3/2 = -22,500^2. 
Thus the amount of salt at any time t > 0 is given by 

22,500^/2 
(2 1 + 50) 3/2 ' 



x = 4( + 100 



Exercises 

t. Assume that the rate at which radioactive nuclei decay is proportional to the 
number of such nuclei that are present in a given sample. In a certain sample 10% 
of the original number of radioactive nuclei have undergone disintegration in a 
period of 100 years. 

(a) What percentage of the original radioactive nuclei will remain after 1000 
years? 

(b) In how many years will only one+ourth of the original number remain? 

2. A certain chemical is converted into another chemical by a chemical reaction. The 
rate at which the first chemical is converted is proportional to the amount of this 
chemical present at any instant. Ten percent of the original amount of the first 
chemical has been converted in 5 min. 

(a) What percent of the first chemical will have been converted in 20 min? 

(b) In how many minutes will 60% of the first chemical have been converted? 
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3. A chemical reaction converts a certain chemical into another chemical, and the 
rate at which the first chemical is converted is proportional to the amount of this 
chemical present at any time. At the end of one hour, 50 gm of the first chemical 
remain; while at the end of three hours, only 25 gm remain. 

(a) How many grams of the first chemical were present initially? 

(b) How many grams of the first chemical will remain at the end of five hours? 

(c) In how many hours will only 2 gm of the first chemical remain? 

4. A chemical reaction converts a certain chemical into another chemical, and the 
rate at which the first chemical is converted is proportional to the amount of this 
chemical present at any time. At the end of one hour, two-thirds kg of the first 
chemical remains, while at the end of four hours, only one-third kg remains. 

(a) What fraction of the first chemical remains at the end of seven hours? 

.(b) When will only one-tenth of the first chemical remain? 

5. Assume that the population of a certain city increases at a ra te proportional to the 
number of inhabitants at any time. If the population doubles in 40 years, in how 
many years will it triple? 

6. The population of the city of Bingville increases at a rate proportional to the 
number of its inhabitants present at any time L If the population of Bingville was 
30,000 in 1970 and 35,000 in 1980, what will be the population of Bingville in 
1990? 

7. In a certain bacteria culture the rate of increase in the number of bacteria is 
proportional to the number present. 

(a) If the number triples in 5 hr, how many will be present in 10 hr? 

(b) When will the number present be 10 times the number initially present? 

8. Assume that the rate of change of the human population of the earth is pro- 
portional to the number of people on earth at any time, and suppose that this 
population is increasing at the rate of 2% per year. The 1979 World Almanac gives 
the 1978 world population estimate as 4,219 million; assume this figure is in fact 
correct. 

fa) Using this data, express the human population of the earth as a function of 
time. 

(b) According to the formula of part (a), what was the population of the earth in 
1950? The 1979 World Almanac gives the 1950 world population estimate as 
2,510 million. Assuming this estimate is very nearly correct, comment on the 
accuracy of the formula of part (a) in checking such past populations. 

(c) According to the formula of part (a), what will be the population of the earth 
in 2000? Does this seem reasonable? 

(d) According to the formula of part (a), what was the population of the earth in 
1900? The 1979 World Almanac gives the 1900 world population estimate as 
1,600 million. Assuming this estimate is very nearly correct, comment on the 
accuracy of the formula of part (a) in checking such past populations. 

(e) According to the formula of part (a), what will be the population of the earth 
in 2100? Does this seem reasonable? 
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9. The human population of a certain island satisfies the logistic law (3.50) with 
k = 0.03, A = 3(10)' 8 , and time t measured in years. 

(a) If the population in 1980 is 200,000, find a formula for the population in 
future years. 

(b) According to the formula of part (a), what will be the population in the year 
2000 ? 

(c) What is the limiting value of the population as t -> oo? 

10. This is a general problem about the logistic law of growth. A population satisfies 
the logistic law (3.50) and has x 0 members at time t 0 * 

(a) Solve the differential equation (3.50) and thus express the population x as a 
function of t. 

(b) Show that as t oo, the population x approaches the limiting value /c/A. 

(c) Show that dx/dt is increasing if x < kjll and decreasing if x > k(2L 

(d) Graph x as a function of t for t > t 0 . 

(e) Interpret the results of parts (b), (c), and (d). 

11. The human population of a certain small island would satisfy the logistic law 
(3.50), with k = 4^0, A « (10p 8 , and t measured in years, provided the annual 
emigration from the island is neglected. However, the fact is that every year 100 
people become disenchanted with island life and move from the island to the 
mainland. Modify the logistic differential equation (3.50) with the given k and A so 
as to include the stated annual emigration. Assuming that the population in 1980 
is 20,000, solve the resulting initial- value problem and thus find the pop- 
ulation of the island as a function of time. 

12. Under natural circumstances the population of mice on a certain island would 
increase at a rate proportional to the number of mice present at any time, 
provided the island had no cats. There were no cats on the island from the 
beginning of* 1970 to the beginning of 1980, and during this time the mouse 
population doubled, reaching an all-time high of 100,000 at the beginning of 
1980. At this time the people of the island, alarmed by the increasing number of 
mice, imported a number of cats to kill the mice. If the indicated natural rate of 
increase of mice was thereafter offset by the work of the cats, who killed 1000 mice 
a month, how many mice remained at the beginning of 1981? 

1 3. An amount of invested money is said to draw interest compounded continuously if 
the amount of money increases at a rate proportional to the amount present. 
Suppose SI 000 is invested and draws interest compounded continuously, where 
the annual interest rate is 6%. 

(a) How much money will be present 10 years after the original amount was 
invested? 

(b) How long will it take the original amount of money to double? 

14. Suppose a certain amount of money is invested and draws interest compounded 
continuously. 

(a) If the original amount doubles in two years, then what is the annual interest 
rate? 



APPLICATIONS OF FIRST-ORDER EQUATIONS 

(b) If the original amount increases 50% in six months, then how long will it take 
the original amount to double? 

15. A tank initially contains 100 gal of brine in which there is dissolved 20 lb of salt. 
Starting at time t = 0, brine containing 3 lb of dissolved salt per gallon flows into 
the tank at the rate of 4 gal/min. The mixture is kept uniform by stirring and the 
well-stirred mixture simultaneously flows out of the tank at the same rate. 

(a) How much salt is in the tank at the end of 10 min? 

(b) When is there 160 lb of salt in the tank? 

16. A large tank initially contains 100 gal of brine in which 10 lb of salt is dissolved. 
Starting at t = 0, pure water flows into the tank at the rate of 5 gal/min. The 
mixture is kept uniform by stirring and the well-stirred mixture simultaneously 
flows out at the slower rate of 2 gal/min. 

(a) How much salt is in the tank at the end of 15 min and what is the 
concentration at that time? 

(b) If the capacity of the tank is 250 gal, what is the concentration at the instant 
the tank overflows? 

17. A tank initially contains 100 gal of pure water. Starting at t = 0, a brine 
containing 4 lb of salt per gallon flows into the tank at the rate of 5 gal/min. The 
mixture is kept uniform by stirring and the well-stirred mixture flows out at the 
slower rate of 3 gal/min. 

(a) How much salt is in the tank at the end of 20 min? 

(b) When is there 50 lb of silt in the tank? 

IS. A large tank initially contains 200 gal of brine in which 15 lb of salt is dissolved. 
Starting at t = 0, brine containing 4 lb of salt per gallon flows into the tank at the 
rate of 3.5 gal/min. The mixture is kept uniform by stirring and the well-stirred 
mixture leaves the tank at the rate of 4 gal/min. 

(a) How much salt is in th^ tank at the end of one hour? 

(b) How much salt is in the tank when the tank contains only 50 gal of brine? 

19. A 500 liter tank initially contains 300 liters of fluid in which there is dissolved 
50 gm of a certain chemical. Fluid containing 30 gm per liter of the dissolved 
chemical flows into the tank at the rate of 4 liters/min. The mixture is kept 
uniform by stirring, and the stirred mixture simultaneously flows out at the rate of 
2.5 liters/min. How much of the chemical is in the tank at the instant it overflows? 

20. A 200 liter tank is initially full of fluid in which there is dissolved 40 gm of .a 
certain chemical. Fluid containing 50 gm per liter of this chemical flows into the 
tank at the rate of 5 liters/min. The mixture is kept uniform by stirring, and the 
stirred mixture simultaneously flows out at the rate of 7 liters/min. How much of 
the chemical is in the tank when it is only half full? 

21. The air in a room whose volume is 10,000 cu ft tests 0.15% carbon dioxide. 
Starting at t = 0, outside air testing 0.05% carbon dioxide is admitted at the rate 
of 5000 cu ft/min. 

(a) What is the percentage of carbon dioxide in the air in the room after 3 nun? 

(b) When does the air in the room test 0.1% carbon dioxide? 
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22. The air in a room 50 ft by 20 ft by 8 ft tests 0.2% carbon dioxide. Starting at t = 0, 
outside air testing 0.05% carbon dioxide is admitted to the room. How many 
cubic feet of this outside air must be admitted per minute in order that the air in 
the room test 0.1% at the end of 30 min? 

23. Newton’s law of cooling states that the rate at which a body cools is proportional 
to the difference between the temperature of the body and that of the medium in 
which it is situated. A body of temperature 80 °F is placed at time t = 0 in a 
medium the temperature of which is maintained at 50 °F. At the end of 5 min, the 
body has cooled to a temperature of 70 °F. 

(a) What is the temperature of the body at the end of 10 min? 

(b) When will the temperature of the body be 60 °F? 

24. A body cools from 60 °C to 50 °C in 15 min in air which is maintained to 30 °C. 
How long will it take this body to cool from 100 °C to 80 °C in air that is main- 
tained at 50 °C? Assume Newton's law of cooling (Exercise 23). 

25. The rate at which a certain substance dissolves in water is proportional at the 
product of the amount undissolved and the difference c x — c 2 > where c x is the 
concentration in the saturated solution and c 2 is the concentration in the actual 
solution. If saturated, 50 gm of water would dissolve 20 gm of the substance. If 
10 gm of the substance is placed in 50 gm of water and half of the substance is 
then dissolved in 90 min, how much will be dissolved in 3 hr? 






CHAPTER FOUR 

Explicit Methods of Solving Higher-Order 
Linear Differential Equations 



The subject of ordinary linear differential equations is one of great theoretical and 
practical importance. Theoretically, the subject is one of simplicity and elegance. 
Practically, linear differential equations originate in a variety of applications to science 
and engineering. Fortunately many of the linear differential equations that thus occur 
are of a special type, linear with constant coefficients, for which explicit methods of 
solution are available. The main purpose of this chapter is to study certain of these 
methods. First, however, we need to consider certain basic theorems that will be used 
throughout the chapter. These theorems are stated and illustrated in Section 4.1, but 
proofs are omitted in this introductory section. By far the most important case is that of 
the second-order linear differential equation, and we shall explicitly consider and 
illustrate this case for each important concept and result presented. In the final section 
of the chapter we return to this fundamental theory and present theorems and proofs in 
this important special case. Proofs in the general case are given in Chapter 11. 



4.1 BASIC THEORY OF LINEAR DIFFERENTIAL EQUATIONS 
A. Definition and Basic Existence Theorem 
DEFINITION 

A linear ordinary differential equation of order n in the dependent variable y and the 
independent variable x is an equation that is inf or can be expressed in , the form 

• -/V !<x) ^ '+ , a,(x)y = /(*), (4.1) 
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where a 0 is not identically zero . We shall assume that a 0y a i9 ...,a„ and F are continuous 
real functions on a real interval a < x < b and that a 0 (x) ^ 0 for any x on a <; x < 6. The 
right-hand member F(x) is called the nonhomogeneous term. If F is identically zero , 
Equation {4.1) reduces to 

d”y d nl y dy 

°o(x) + ”’ + «*- 1 w ^ + a "My = 0 ( 4 - 2 > 

and is then called homogeneous. 

For n =: 2, Equation (4.1) reduces to the second-order nonhomogeneous linear 
differential equation 

d 2 v dv 

a o(*) + a 2 (x)y = F(x) (4.3) 

and (4.2) reduces to the corresponding second-order homogeneous equation 

a o( x ) + a ^% + a 2( x )y = a < 44 ) 

Here we assume that a 09 a l9 a l9 and F are continuous real functions on a real interval 
a<x < b and that a Q {x) ^ 0 for any x on a < x ^ b. 



► Example 4.1 



The equation 



d 2 y * dy 



+ 3x — + x 3 y = e x 



is a linear ordinary differential equation of the second order. 



► Example 4.2 



The equation 



d 3 y d 2 y 2 dy 

-r-r + x — f + 3x 2 5y = sin x 

dx 3 dx 2 dx 



is a linear ordinary differential equation of the third order. 

We now state the basic existence theorem for initial-value problems associated with 
and nth-order linear ordinary differential equation: 

THEOREM 4.1 
Hypothesis 

I . Consider the nth-order linear differential equation 



<ioto xt + fliM -r^T + — + A. - lto 37 + a-toy - Ffr). 



f 
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where a 0 ,a it ... t a n and F are continuous real functions on a real interval a < x < b and 
u 0 (x) ^ 0 for any x on a < x < b. 

2. Let x 0 be any point of the interval a <, x $ h, and let c 0y c l ,...,c n _ l ben arbitrary 
real constants . 



Conclusion. There exists a unique solution f of (4.1) such that 
fM = CqJ'(x o) = c„...,/ ( "'°(x 0 ) = c R , u 
and this solution is defined over the entire interval a < x < 6. 

Suppose that we are considering an nth-order linear differential equation (4.1), the 
coefficients and nonhomogeneous term of which all possess the continuity require- 
ments set forth in Hypothesis 1 of Theorem 4.1 on a certain interval of the x axis. Then, 
given any point x 0 of this interval and any n real numbers c 0 ,Cj , . . . , c n _ l9 the theorem 
assures us that there is precisely one solution of the differential equation that assumes 
the value c Q at x = x 0 and whose kth derivative assumes the value c k for each k - 
1 at x - x 0 . Further, the theorem asserts that this unique solution is 
defined for all x in the above-mentioned interval. 

For the second-order linear differential equation, 

a°(x) + tfi(x) ~ + a 2 (x)y = F(x), (4.3) 

the requirements of Hypothesis l of Theorem 4.1 are that a^ y a^a 2> and F be 
continuous on a real interval a< x < b and that a 0 (x) £ 0 for any x on this interval. 
Then, if x 0 is any point of the interval a < x < b and c 0 and Cj are any two real numbers, 
the theorem assures us that there is precisely one solution / of the second-order 
differential equation (4.3) which assumes the vqlue c Q at x = x 0 and whose first 
derivative assumes the value c t at x — x 0 : 

f(x 0 ) = c 0y f'(x 0 ) -c x , (4.5) 

Moreover, the theorem asserts that this unique solution / of Equation (4.3) which 
satisfies conditions (4.5) is defined for all x on the interval a < x < b. 



► Example 4.3 



Consider the initial-value problem 



dx 2 + dx 



+ x*y — e x 



y( 1 ) = 2 , 

/(!)= “5. 

The coefficients 1, 3x, and x 3 , as well as the nonhomogeneous term e x > in this second- 
order differential equation are all continuous for all values of x, — oo < x < oo. The 
point x 0 here is the point 1, which certainly belongs to this interval; and the real 
numbers c 0 and are 2 and —5, respectively. Thus Theorem 4.1 assures us that a 
solution of the given problem exists, is unique, and is defined for all x, — oo < x < oo. 
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► Example '4.4 

Consider the initial- value problem 

2 fl + x d -l 

Jx> Jx 1 



„ , dy 

+ Sx'' 5y — sin x, 



y(4) = 3, 

/( 4) = 5, 

y”( 4 ) = -l 

Here we have a third-order problem. The coefficients 2,x, 3x 2 , and —5, as well as the 
nonhomogeneous term sin x, are all continuous for all x, — oo < x < oo. The point 
x 0 = 4 certainly belongs to this interval; the real numbers c 0 ,c 1 , and c 2 in this prob- 
lem are 3, 5, and — respectively. Theorem 4.1 assures us that this problem also has 
a unique solution which is defined for all x, — oo < x < oo. 

A useful corollary to Theorem 4.1 is the following: 



COROLLARY 

Hypothesis. Let f be a solution of the nth-order homogeneous linear differential 
equation 

W + a ' (x) lx^ + • ■ ■ ■ ■ + 1 (*> % + a « <*>* = 0 < 4 - 2 ) 



such that 



fix 0 ) = 0,f’{x o ) - 0,...,/<"- l >(x 0 ) = 0, 



where x 0 is a point of the interval a < x < b in which the coefficients a 0y a l> ... 9 a„ are 
all continuous and u 0 (x) ^ 0. 

Conclusion. Then f(x) — 0 for all x on a < x < b. 

Let us suppose that we are considering a homogeneous equation of the form (4.2), all 
the coefficients of which are continuous on a certain interval of the x axis. Suppose 
further that we have a solution / of this equation which is such that / and its first n - 1 
derivatives all equal zero at a point x 0 of this interval. Then this corollary states that 
this solution is the “trivial” solution / such that f(x) = 0 for all x on the above- 
mentioned interval. 



► Example 4.5 

The unique solution / of the third-order homogeneous equation 

d 3 y „ d 2 y A dy - 

^ + 2 ^ + 4x ^ + x ^ 0 ’ 
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which is such that 

/(2) = f\2) = f; 2) = 0, 

is the trivial solution / such that f(x) = 0 for all x. 

B. The Homogeneous Equation 

We now consider the fundamental results concerning the homogeneous equation 

d n v d n ~ l v dv 

M*) dx!" + dx 1 ^ + " +a "~ + a *M y = °- ( 4 - 2 > 

We first state the following basic theorem: 



THEOREM 4.2 BASIC THEOREM ON LINEAR HOMOGENEOUS 
DIFFERENTIAL EQUATIONS 

Hypothesis. Let fi>fi 9 ..~ 9 f m be any m solutions of the homogeneous linear differ- 
ential equation (4.2). 

Conclusion. Then c t f { + c 2 / 2 + h c m f m is also a solution of (4.2), where 

c l ,c 2 ,... J c m are m arbitrary constants . 

Theorem 4.2 states that if m known solutions of (4.2) are each multiplied by an 
arbitrary constant and the resulting products are then added together, the resulting 
sum is also a solution of (4.2). We may put this theorem in a very simple form by means 
of the concept of linear combination, which we now introduce. 



DEFINITION 

If are m given functions , and c x ,c 2 ,,.. t c m are m constants , then the ex- 

pression 

c,/i + c 2 / 2 + + 

is called a linear combination of f x ,f 2 >- * ■ ♦/«• 

In terms of this concept. Theorem 4.2 may be stated as follows: 



THEOREM 4.2 (RESTATED) 

Any linear combination of solutions of the homogeneous linear differential equation (4.2) 
is also a solution of (4.2). 

In particular, any linear combination 

c lfl + c 2 f% + ’ ' ' \ c mfm 



4.1 
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of m solutions ...,/ w of the second- order homogeneous linear differential 
equation 

*o(*) lx 1 + ° 1 ' W %. + ° 2 * 0 ( 4 - 4 ) 



is also a solution of (4.4). 

► Example 4.6 

The student will readily verify that sin x and cos x are solutions of 

d 2 y 

lx* + y = °- 



Theorem 4.2 states that the linear combination c x sin x + c 2 cos x is also a solution for 
any constants c x and c 2 . For example, the particular linear combination 



5 sin x + 6 cos x 



is a solution. 



► Example 4.7 

The student may verify that e x , e~ x y and e 2x are solutions of 
d 2 y ~ d 2 y dy ■ 

Theorem 4.2 states that the linear combination c v e x + c 2 e~ x + c 3 e 2x is also a solution 
for any constants c u c 2l and c 3 . For example, the particular linear combination 

2e* — 4- \e 2x 

is a solution. 

We now consider what constitutes the so-called general solution of (4.2). To 
understand this we first introduce the concepts of linear dependence and linear 
independence . ... . .. 



DEFINITION 



The n functions /i,/ 2 , ■**»/* are called linearly dependent on a < x < b if there exist 
constants C|,c 2 > - *»c n , not all zero, such that 

Cl fi(x) + c 2 f 2 (x) + ••• + cj„(x) = 0 
for all x such that a < x < b, 

In particular , two functions f x and f z are linearly dependent on a < x £ b if there exist 
constants c u c 2i not both zero, such that 



c't/iW + c 2 / 2 (x) = 0 



for all x such that a < x < b. 
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► Example 4,8 

We observe that x and 2x are linearly dependent on the interval 0 3 x <, 1. For there 
exist constants c x and c 2 , not both zero , such that 

c x x 4 c 2 ( 2x) = 0 

for all x on the interval 0 < x ^ 1. For example, let c Y =* 2, c 2 - — t. 

► Example 4.9 

We observe that sin x, 3 sin x, and —sin x are linearly dependent on the interval 
— 1 ^ x < 2. For there exist constants c 1 ,c 2 ,c 3 , not ail zero , such that 

' \ 

c l sin x 4- c 2 (3 sin x) 4- c 3 (-sin x) = 0 
for all x on the interval — 1 < x < 2. For example, let c y = l,c 2 — 1, c 3 = 4. 



DEFINITION 

The n functions f i9 f 2 ,... f f m are called linearly independent on the interval a < x < b if 
they are not linearly dependent there. That is, the functions / u / 2 ,...,/„ arc linearly 
independent on a <, x <>b if the relation 

cJx(x) 4 c 2 / 2 (x) 4 — 4 c„/„(x) = 0 

/or all x such that a <x < b implies that 

Ci =c 2 = * — — 0. 

//I otAcr words, t/ie on/y /inear combination of f u f 2i ...,f n that is identically zero on 
a< x < bis the trivial linear combination 

°7i+0‘/ 2 + “' + o-/„. 

In particular , two /unctions / x and / 2 ore linearly independent on a < x < b if the 
relation 

Ci/i(x) + c 2 / 2 (x) = 0 
/or all x on a < x < b implies that 

c i - c 2 - 0. 



► Example 4.10 

We assert that x and x 2 are linearly independent on 0 < x < 1, since c,x 4 c 2 x 2 — 0 
for all x on 0 < x < 1 implies that both c t = 0 and c 2 - 0. We may verify this in 
the following way. We differentiate both sides of qx + c 2 x 2 - 0 to obtain c t + 
2c 2 x = 0, which must also hold for all x on 0 ^ x < 1. Then from this we also have 
c t x + 2c 2 x 2 = 0 for all such x. Thus we have both 

c,x 4 c 2 x 2 - 0 and c t x 4 2c 2 x 2 = 0 



(4.6) 
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for all x on 0 £ x < 1. Subtracting the first from the second gives c 2 x 2 — 0 for all x on 
0 £ x < 1, which at once implies c 2 = 0. Then either of (4.6) show similarly that c x = 0. 

The next theorem is concerned with the existence of sets of linearly independent 
solutions of an nth-order homogeneous linear differential equation and with the 
significance of such linearly independent sets. 

THEOREM 4.3 

The nth-order homogeneous linear differential equation 

«oW ~ + «, M ~T + ”■+«.- iW % + aMy = 0 (4.2) 

always possesses n solutions that are linearly independent . Further , if are n 

linearly independent solutions of (4.2), then every solution f of (4.2) can be expressed as a 
linear combination 



C \f\ + C lfl + **' + Cnfn 

of these n linearly independent solutions by proper choice of the constants 

Given an nth-order homogeneous linear differential equation, this theorem assures 
us first that a set of n linearly independent solutions actually exists. The existence of 
such a linearly independent set assured, the theorem goes on to tell us that any solution 
whatsoever of (4.2) can be written as a linear combination of such a linearly inde- 
pendent set of n solutions by suitable choice of the constants c l ,c 2 ,* . .,c„. 

For the second-order homogeneous linear differential equation 

d 2 y dv 

fl ° W + ^ lx + ° 2 ^ y " °* ( 4 4 ) 

Theorem 4.3 first assures us that a set of two linearly independent solutions exists. The 
existence of such a linearly independent set assured, let /, and / 2 be a set of two linearly 
independent solutions. Then if / is any solution of (4.4), the theorem also assures us 
that / can be expressed as a linear combination c l f l + c 2 f 2 of the two linearly 
independent solutions f x and f 2 by proper choice of the constants c x and c 2 . 



► Example 4.1 1 

We have observed that sin x and cos x are solutions of 

d 2 y 

d + >-° (4 - 7 > 

for all x, — oo < x < oo. Further, one can show that these two solutions are linearly 
independent. Now suppose / is any solution of (4.7). Then by Theorem 4.3 / can be 
expressed as a certain linear combination c x sin x + c 2 cos x of the two linearly 
independent solutions sin x and cos x by proper choice of c x and c 2 . That is, there exist 
two particular constants c x and c 2 such that 

f(x) = c x sin x + c 2 cos x 



(4.8) 



110 EXPLICIT METHODS OF SOLVING HIGHER-ORDER LINEAR DIFFERENTIAL EQUATIONS 



for all x, — oo < x < oo. For example, one can easily verify that/(x) = sin(x 4- n/6) is a 
solution of Equation (4.7). Since 

. ( n\ n m n J 3 . 1 

sin x + - - sin x cos - + cos x sin — - ~r- sin x + - cos x, 

\ 6 / 6 6 2 2 

we see that the solution sin(x + 7r/6) can be expressed as the linear combination 

yfi . , 

s,n x 4* 2COS X 

2 2 

of the two linearly independent solutions sin x and cos x, Note that this is of the form in 
the right member of (4.8) with c x = <^3/2 and c 2 = 2 * 

Now let f x ,f 2 be a set of n linearly independent solutions of (4.2). Then by 
Theorem 4.2 we know that the linear combination 

c Jt + c 2 f 2 +•' + cj„, ( 4 -9) 

where c l ,c 2> . ..,c„ are n arbitrary constants, is also a solution of (4.2), On the other 
hand, by Theorem 4.3 we know that if / is any solution of (42), then it can be expressed 
as a linear combination (4.9) of the n linearly independent solutions /i>/ 2 , . . ■,/„> by a 
suitable choice of the constants c 1 ,c 2 ,..*,c„. Thus a linear combination (4.9) of the 

n linearly independent solutions f ly f 2 /„ in which Cj,c 2 ,. . .,c„ are arbitrary 

constants must include all solutions of (4.2). For this reason, we refer to a set of n 
linearly independent solutions of (4.2) as a “fundamental set” of (4.2) and call a 
“general” linear combination of n linearly independent solutions a “general solution” 
of (4.2), in accordance with the following definition: 



DEFINITION 

are n linearly independent solutions of the nth-order homogeneous linear 
differential equation 

«o(*) + a t (x) + •" + «»- lW ^ + a J x )y = 0 (4.2) 

on a <: x < b, then the set ■ ■>/„ is called a fundamental set of solutions of (4.2) 

and the function f defined by 

f(x) = c i f i (*) + Cifi (*) + ■•■+ cj„(x) y a<x<b , 

where c ly c 2t ... t c n are arbitrary constants , is called a general solution of (4.2) on 
a < x < b. 

Therefore, if we can find n linearly independent solutions of (4.2), we can at once 
write the general solution of (4.2) as a general linear combination of these n solutions. 
For the second- order homogeneous linear differential equation 

d 2 v dy 

a °Wd^ + a iMfa +a2 M y=:() ' ( 4 . 4 ) 

'■indamental set consists of two linearly independent solutions. If /, and f 2 are a 
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fundamental set of (4.4) on a < x < b, then a general solution of (4.4) on a < x < b is 
defined by 

Cifi (x) 4- c 2 / 2 (x), a<x<b , 

where c x and c 2 are arbitrary constants. 

► Example 4.12 

We have observed that sin x and cos x are solutions of 

d 2 y 

j£ +y =° 

for all x, —oo < x < oo. Further, one can show that these two solutions are linearly 
independent. Thus, they constitute a fundamental set of solutions of the given 
differential equation, and its general solution may be expressed as the linear 
combination 

c t sin x + c 2 cos x, 

where and c 2 are arbitrary constants. We write this as y = c r sin x 4- c 2 cos x. 



► Example 4.1 3 

The solutions e x , e~* t and e 2x of 

may be shown to be linearly independent for all x, - go < x < oo. Thus, e\e x , and e 2x 
constitute a fundamental set of the given differential equation, and its general solution 
may be expressed as the linear combination 

c l e x + c 2 e~ x 4- c z e 2x t 

where c u c ls and c 3 are arbitrary constants. We write this as 

y = c x e x + c 2 e~ x 4- c 3 e 2x . 

The next theorem gives a simple criterion for determining whether or not n solutions 
of (4.2) are linearly independent. We first introduce another concept. 



DEFINITION 



Let fu f 2 y . ,/„ be n real functions each of which has an (n — I)st derivative on a real 
interval a < x < b. The determinant 






h fz ... /. 

f\ fz -. f'n 

f(r>- 1) fi ft-i) 
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in which primes denote derivatives , is called the Wronskian of these n functions. We 
observe that W(f l7 f 2i . . . 9 f„) is itself a real function defined on a < x <b. Its value at x 
is denoted by W(f u f 2 ,., . ,/„)(*) or by W[f y (x)J x (x )+. . . ,/ n M3. 



THEOREM 4.4 

The n solutions f u f 2 of the nth-order homogeneous linear differential equation 
(4.2) are linearly independent on a < x < b if and only if the Wronskian of f u f l9 , . . >/ n 
is different from zero for some x on the interval a < x < b. 

We have further: 



THEOREM 4.5 

The Wronskian of n solutions --,/„<>/ (4 .2) is either identically zero on a ^ x < b 

or else is never zero on a < x < b. 



Thus if we can find n solutions of (4.2), we can apply the Theorems 4.4 and 45 to 
determine whether or not they are linearly independent. If they are linearly inde- 
pendent, then we can form the general solution as a linear combination of these n 
linearly independent solutions. 

In the case of the general second- order homogeneous linear differential equation 
■3°W M ^ + a 2 {x)y = 0, (4.4) 



the Wronskian of two solutions f x and f 2 is the second-order determinant 






/. 

A 



h 

fi 






By Theorem 4.4, two solutions f x and f 2 of (4.4) are linearly independent on a < x < b 
if and only if their Wronskian is different from zero for some x on a < x < b; and by 
Theorem 45, this Wronskian is either always zero or never zero on a < x < b. Thus if 
W\_f W* / 2 W] ^0 on a<x<b, solutions A and f 2 of (44) are linearly independent 
on a < x < b and the general solution of (4.4) can be written as the linear combination 



Cl/lW + c 2 / 2 (*). 



where c t and c 2 are arbitrary constants. 



► Example 4.14 



We apply Theorem 4.4 to show that the solutions sin x and cos x of 
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are linearly independent. We find that 



lV(sin x, cos x) = 



sin x 
cos x 



cos x 
— sin x 



— —sin 2 x — cos 2 x 



-1 *0 



for all real x. Thus* since W (sin x, cos x) ^ 0 for all real x, we conclude that sin x and 
cos x are indeed linearly independent solutions of the given differential equation on 
every real interval. 



► Example 4.15 

The solutions e x 7 e~ x , and e 2x of 

_ 2 — ^ - — + 2y = 0 

dx 3 dx 2 dx + } 

are linearly independent on every real interval* for 





e x 


e~ x 


e 2x 




1 


1 


1 


W{e x ,e- X ,e 2x ) = 


e x 


—e~ x 


2e 2x 


— e 2x 


1 


-1 


2 




e x 


e’ x 


4e 2x 




1 


1 


4 



= —6e 2x / 0 



for all real x. 



Exercises 

1. Theorem 4.1 applies to one of the following problems but not to the other. 
Determine to which of the problems the theorem applies and state precisely the 
conclusion which can be drawn in this case. Explain why the theorem does not 
apply to the remaining problem. 

(a) & +5 £ + 6>, = e ‘’ > ,(0) = 5 ’ /(°) = 7 - 

(b) S +5 £ +6) ' =ex ’ y{0)=5> ^ ,)=7 - 

2. Answer orally: What is the solution of the following initial-value problem? Why? 

0 + x£ + x^ = °, y(l) = 0, /(l)-Q. 

3. Prove Theorem 4.2 for the case m = n = 2. That is, prove that if /Jx) and / 2 (x) 
are two solutions of 

d 2 v dv 

a o(x)^ + a Mj; + a 2 (x)y = 0 y 

then c x f x (x) + c 2 / 2 (x) is also a solution of this equation, where c x and c 2 are 
arbitrary constants. 
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! 



4. Consider the differential equation 



d 2 y 4 dy 
dx 2 dx 



+ 3y = 0. 



(A) 



(a) Show that each of the functions e* and e 3x is a solution of differential 
equation (A) on the interval a < x <; fc, where a and b are arbitrary real 
numbers such that a < b. 

(b) What theorem enables us to conclude at once that each of the functions 

5e * + 2e 3x y 6e x — 4e 3x , and — le* + 5e 3x 

is also a solution of differential equation (A) on a <> x < 6? 

(c) Each of the functions 

3e x , -4e x , 5e\ and 6e x 

is also a solution of differential equation (A) on a <, x < b. Why? 

5. Again consider the differential equation (A) of Exercise 4. 

(a) Use the definition of linear dependence to show that the four functions of 
part (c) of Exercise 4 are linearly dependent on a < x < b. 

(b) Use Theorem 4.4 to show that each pair of the four solutions of differential 
equation (A) listed in part (c) of Exercise 4 are linearly dependent on 
a £ x < b. 



6. Again consider the differential equation (A) of Exercise 4. 

(a) Use the definition of linear independence to show that the two functions e x 
and e 3x are linearly independent on a < x < b. 

(b) Use Theorem 4.4 to show that the two solutions e x and e 3x of differential 
equation (A) are linearly independent on a < x < b. 

7. Consider the differential equation 



dx 2 dx 



+ 6y = 0. 



(a) Show that e 2x and e 2x are linearly independent solutions of this equation on 
the interval — oo < x < oo. 



(b) Write the general solution of the given equation, 

(c) Find the solution that satisfies the conditions y(0) — 2, y'(0) = 3. Explain 
why this solution is unique. Over what interval is it defined? 



8. Consider the differential equation 




d^ 

dx 2 




+ y « 0. 



(a) Show that e x and xe x are linearly independent solutions of this equation on 
the interval — oo < x < oo. 



(b) Write the general solution of the given equation. 

(c) Find the solution that satisfies the condition y(0) = 1 , y '(0) = 4. Explain why 
this solution is unique. Over what interval is it defined? 
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9. Consider the differential equation 



*’£;-2*£ + 2r-ft 



dx 2 



dx 



(a) Show that x and x 2 are linearly independent solutions of this equation on 
the interval 0 < x < oo. 

(b) Write the general solution of the given equation. 

(c) Find the solution that satisfies the conditions y(l) = 3, y'(l) = 1 Explain 
why this solution is unique. Over what interval is this solution defined? 

10. Consider the differential equation 



x 



2 



d 2 y dy 

d? + x Tx 



- Ay - 0. 



(a) Show that x 2 and 1/x 2 are linearly independent solutions of this equation on 
the interval 0 < x < oo. 

(b) Write the general solution of the given equation. 

(c) Find the solution that satisfies the conditions y( 2) = 3, y'(2) — — 1. Explain 
why this solution is unique. Over what interval is this solution defined? 



1 1 . Consider the differential equation 



(a) 

(b) 

(c) 

(d) 



d 2 y 
dx 2 



dy 

5/ + 4y = 0. 

dx 



Show that each of the functions e x , e 4x , and 2e x — 3e 4x is a solution of this 
equation on the interval - oo < x < oo. 

Show that the solutions e x and e 4x are linearly independent on - oo < 

x < oo. 

Show that the solutions e x and 2e x — 3e 4x are also linearly independent on 



-OO < X < CO. 

Are the solutions e 4x and 2e x - 3e 4x still another pair of linearly indepen 
dent solutions on — oo < x < col Justify your answer. 



12. Given that e~ x , e 3x , and e 4x are all solutions of 



d 3 y 
dx 3 



, d l y dy 

6 dS + 5 i; + 



I2y - 0, 



show that they are linearly independent on the interval — oo < x < oo and write 
the general solution. 

13. Given that x, x 2 , and x 4 are all solutions of 



, d 3 y 
dx c 3 



— 4x : 



d 2 y 

dx 2 



dy 

+ 8x 8y = 0, 

dx 



show that they are linearly independent on the interval 0 < x < oo and write the 
general solution. 



L 
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C. Reduction of Order 

In Section 4.2 we shall begin to study methods for obtaining explicit solutions of 
higher-order linear differential equations. There and in later sections we shall find that 
the following theorem on reduction of order is often quite useful. 

THEOREM 4.6 

Hypothesis. Let f be a nontrivial solution of the nth-order homogeneous linear dif- 
ferential equation 

d n v d nl y dy * 

a ° ix) d£ + a ' {x) d^ + ' + a *~ l(x) Jx + a ” (x)y = a (4 ‘ 2) 

Conclusion. The transformation y - f(x) v reduces Equation (4.2) to an (« - \)st- 
order homogeneous linear differential equation in the dependent variable w =* dv/dx. 

This theorem states that if one nonzero solution of the nth-order homogeneous lin- 
ear differential equation (4.2) is known, then by making the appropriate transforma- 
tion we may reduce the given equation to another homogeneous linear equation that is 
one order lower than the original. Since this theorem will be most useful for us in 
connection with second-order homogeneous linear equations (the case where n = 2), 
we shall now investigate the second-order case in detail. Suppose / is a known 
nontrivial solution of the second-order homogeneous linear equation 

°°W^T + o, W^ + fl 2 (x)y = 0. (4.10) 

Let us make the transformation 

y = /(x)y, (4.11) 

where / is the known solution of (4.10) and v is a function of x that will be determined. 
Then, differentiating, we obtain 

§-/(->£ + /», (4.12) 

d 2 y d 2 v dv 

d? = m d?‘ + 2f ' (x) Tx + r(x)v - (413) 

Substituting (4.11), (4.12), and (4.13) into (4.10), we obtain 

flo (*)|V(x) + 2 f'(x) ^ + f"(x)v J + a, (x)|V(x) ^ + f'{x)v J + a 2 {x)f(x)v = 0 

or 

"o.W/W^r + [2a 0 (*)/'(*) + a,(x)/(x)] 

+ [«o( x )/"M + a i( x )f‘( x ) + = 0. 
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Since / is a solution of (4.10), the coefficient of t; is zero, and so the last equation reduces 
to 



tfo (*)/(*) + [2 a 0 (x)f'(x) + <Mx)/(x)] — = 0. 

Letting w = dv/dx , this becomes 

dw 

floW/W ^ + [2a 0 ( x )/( x ) + *i (*)/(*)>' = 0. (4.14) 

This is a ^irsf-order homogeneous linear differential equation in the dependent variable 
w'. The equation is separable; thus assuming /(x) / 0 and a 0 (x) # 0, we may write 



= _r 2 /M + fiw] dx 

L m oowj 



Thus integrating, we obtain 



In 1 w| = — ln[/(x)] 2 - -^r-rdx + In |c| 






[/(- V )] 2 

This is the general solution of Equation (4.14); choosing the particular solution for 
which c = 1, recalling that dv/dx — vv, and integrating again, we now obtain 



Finally, from (4,1 1), we obtain 



,exp 


[-J 


;a«*l 

fl oM J 


[/(X)] 2 



y = A*) 



hiih 



C/(X)] 2 



Th^ function defined in the right member of (4. 1 5), which we shall henceforth denote by 
g , is actually a solution of the original second-order equation (4. 10). Furthermore, this 
new solution g and the original known solution / are linearly independent, since 

wlf(Mx \ = /M 9(x) _ f(x) f(x)v 

J ’ fix) g'(x) ~ f’(x) f(x)v' + f'(x)v 



.CK«))V.exp[-{^^]^0. 



Thus the linear combination 



Cj+Crf 



is the general solution of Equation (4,10). We now summarize this discussion in the 
following theorem. 
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THEOREM 4.7 

Hypothesis. Let f be a nontrivial solution of the second-order homogeneous 
linear differential equation 



d 2 v -r 

a o{x)-^ + a 1 (x)— + a 2 (x)y = 0. 



)d T 



(4.10) 



Conclusion 1, The transformation y — f(x) v reduces Equation (4 AO) to the first- 
order homogeneous linear differential equation 

a 0 M/M + [2 a 0 M/'M + a , M/M] w = 0 (4.14) 

in the dependent variable w, where w *= dv/dx. 



Conclusion 2. The particular solution 




of Equation (4 A 4) gives rise to the function v y where 




The function g defined by g(x) = f(x)v(x) is then a solution of the second-order equation 
(4 AO). 



Conclusion 3. The original known solution f and the “new” solution g are linearly 
independent solutions of (4 AO), and hence the general solution of (4. 10) may be expressed 
as the linear combination 



Cif+c 2 g: 



Let us emphasize the utility of this theorem and at the same time clearly recognize its 
limitations. Certainly its utility is by now obvious. It tells us that if one solution of the 
second-order equation (4.10) is known, then we can reduce the order to obtain a linearly 
independent solution and thereby obtain the general solution of (4.10). But the 
limitations of the theorem are equally obvious. One solution of Equation (4.10) must 
already be known to us in order to apply the theorem. How does one “already know” a 
solution? In general one does not. In some cases the form of the equation itself or 
related physical considerations suggest that there may be a solution of a certain special 
form: for example, an exponential solution or a linear solution. However, such cases are 
not too common and if no solution at all can be so ascertained, then the theorem will 
not aid us. 

We now illustrate the method of reduction of order by means of the following 
example. 
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► Example 4.16 

Given that y = x is a solution of 



J 2 y ~<iy 



{x2 + l) l^~ 2 x £ + 2>, = 0 ’ < 416 ) 

find a linearly independent solution by reducing the order. 

Solution. First observe that y = x does satisfy Equation (4.16). Then let 

y = xv. 

Then 

dy dv . d 2 y d 2 v . 

+ " j?-*3? + 2 S' 

Substituting the expressions for y, dy/dx, and d 2 y/dx 2 into Equation (4.16), we obtain 
/ 2 n/ dtA f dv \ 

<r + 1 VE I+ (*E + 7 + 



or 



. , <? 2 i? * <iu 

x,x +1) 5? + 2 &- a 

Letting w - dv/dx we obtain the yirst-order homogeneous linear equation 

x(x 2 + 1)^ + 2w - 0. 
dx 

Treating this as a separable equation, we obtain 

dw 2 dx 

w 



or 



x(x 2 + 1) 
2x 



dw (2 2x \ . 

— — I 1 — 2 7 f ) 

w \ x x* + l } 



Integrating, we obtain the general solution 

c(x 2 + 1) 



w = - 

Choosing c = 1, we recall that dv/dx — w and integrate to obtain the function v given 
by 

u(x) = x — — . 

X 

Now forming g = fv y where / (x) denotes the known solution x, we obtain the function g 
defined by 



9(x) 



-*H)- 
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By Theorem 4.7 we know that this is the desired linearly independent solution. The 
general solution of Equation' (4. 16) may thus be expressed as the linear combination 
c l x *f c 2 (x 2 — 1) of the linearly independent solutions / and g. We thus write the 
general solution of Equation (4.16) as 

y - CjX + c 2 {x 2 - 1). 



D. The Nonhomogeneous Equation 



We now return briefly to the nonhomogeneous equation 

a oW ^ (*) ^ + ' - ' + - , (x) ~ + <U.x), = F(x). (4.1) 

The basic theorem dealing with this equation is the following. 



THEOREM 4.8 
Hypothesis 

(7) Let v be any solution of the given ( nonhomogeneous ) nth-order linear differential 
equation (4.7). (2) Let u be any solution of the corresponding homogeneous equation 

flo(x) + «|W + ' • • + a n- l(x) + «„(*)>’ = o. (4.2) 

Conclusion. Then u + v is also a solution of the given (nonhomogeneous) equation 
(4.1). 



► Example 4.17 

Observe that y = x is a solution of the nonhomogeneous equation 

d 2 y 



2 +y = x. 



dx 

and that y = sin x is a solution of the corresponding homogeneous equation 

d 2 y 



Er + r-o. 

Then by Theorem 4.8 the sum 

sin x + x 

is also a solution of the given nonhomogeneous equation 

d 2 y 

±e +y=x - 

The student should check that this is indeed true. 



